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y  16.  AbstruCt 

This  report  documents  the  results  of  the  Federal  Aviation  Administration  (FAA) /  Environmental  Protection  Agency 
(EPA)  air  quality  study  which  has  been  conducted  to  assess  the  Impact  of  aircraft  emissions  of  carbon  monoxide 
(CO)  ,  hydrocarbons  (HC)  ,  and  oxides  of  nitrogen  (NO*)  in  the  vicinity  of  airports.  This  assessment  includes  the 
results  of  recent  modeling  and  monitoring  efforts  at  Washington  National  (1)CA)  ,  Los  Angeles  International  (LAX), 
Dulles  International  (1AD) ,  and  Lakeland,  Florida  airports  and  an  updated  modeling  of  aircra! t  generated  pollution 
at  LAX,  John  F.  Kennedy  (JFK)  and  Chicago  O'Hare  (ORD)  airports.  The  Airport  Vicinity  Air  Po.lution  (AVAP)  model 
which  was  designed  for  use  at  civil  airports  was  used  in  this  assessment.  In  addition  th.  results  of  the 
application  of  the  military  version  of  the  AVAP  model  the  Air  Quality  Assessment  Model  (AQAM) ,  are  summarized  ^ 

Toth  the  results  of  the  pollution  monitoring  analyses  in  Volume  I  and  the  modeling  studies  in  Volume  II  suggest 

o  Maximum  hourly  average  CO  concentrations  from  aircraft  are  unlikely  to  exceed  5  parts  per  million  (ppm) 
in  areas  of  public  exposure  and  are  thus  small  In  comparison  to  the  National  Ambient  Air  Quality  Standard  of  35  ppm. 

o  Maximum  hourly  HC  concentrations  from  aircraft  can  exceed  0.25  ppm  over  an  area  several  times  the  size  of 
the  airport. 

o  While  annual  average  NO2  concentrations  from  aircraft  are  estimated  to  contribute  only  10  to  20  percent  ot 
the  NAAQS  limit  level,  these  concentrations,  when  averaged  over  a  one  hour  time  periodare  estimated  to  produce  con¬ 
centrations  as  high  as  0.5  ppm  if  one  assumes  that  all  engine  produced  NO  Is  converted  to  NO2  by  the  time  these 
emissions  reach  places  of  public  exposure.  This  value  is  at  the  upper  end  of  the  concentration  range  being  considered 
for  the  short  term  N02  standard  presently  under  review  and  can  not  be  ignored.  However  a  decision  on  the  Influence 
of  aircraft  emissions  on  a  short  term  standard  must  await  the  release  of  this  standard  and  the  result*  of  research 
under  way  to  determine  the  rate  at  which  engine  produced  NO  converts  to  N02  in  the  ambient  air. 

This  report,  IMPACT  OF  AIRCRAFT  EMISSIONS  ON  AIR  QUALITY  IN  THE  VICINITY  OF  AIRPORTS,  la  printed  in 
two  volumes  with  the  following  subtitles;  VOLUME  I:  Recent  Airport  Measurement  Programs.  Data 
Analyses,  and  Sub-Model  Development;  and  VO^bME  XI:  An  Updated  Model  Assessment  of  Aircraft 
Generated  Air  Pollution  at  LAX.  JFK,  and  ORD . 
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PREFACE 


This  document  has  been  prepared  by  the  Energy  and  Environmental  Systems 
Division  of  Argonne  National  Laboratory  at  the  request  of  the  Federal  Aviation 
Administration.  The  report  attempts  to  realistically  simulate  the  air 
quality  impact  of  aircraft  in  and  around  the  airport  property  under  adverse 
dispersion  conditions.  No  attempt  has  been  made  to  include  the  effect  of 
non-aircraft  sources. 
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1  INTRODUCTION  AND  SUMMARY 


1 . 1  INTRODUCTION 

The  1970  ammendments  to  the  Clean  Air  Act  (U.S.  Congress,  1970) 
directed  the  Environmental  Protection  Agency  (EPA)  to  establish  emission 
standards  for  aircraft  and  aircraft  engines  if  such  emissions  are  judged  to 
cause  or  are  likely  to  cause  or  contribute  to  air  pollution  which  endangers 
the  public  health.  The  1970  amendments  also  directed  the  EPA  to  conduct  a 
study  of  the  extent  to  which  aircraft  emissions  affect  air  quality  in  Air 
Quality  Control  Regions  (AQCRs)  throughout  the  United  States.  Based  upon 
information  available  in  the  early  1970's  the  EPA  judged  that  aircraft  and 
aiports  were  then  or  were  projected  to  be  significant  sources  of  emissions  of 
carbon  monoxide  (CO),  hydrocarbons  (HC),  nitrogen  oxides  (N0X) ,  and  smoke  in 
some  of  the  AQCRs  where  the  National  Ambient  Air  Quality  Standards  (NAAQS) 
were  being  violated.  Therefore,  engine  emission  standards  were  promulgated  on 
July  17,  1973  for  Commercial  and  General  Aviation  aircraft.  The  basis  for 
these  emission  standards  was  the  air  quality  and  technology  assessments  made 
during  the  early  1970' s.  Since  that  time  major  advancements  have  been  made  in 
the  techniques  for  monitoring  and  modeling  aircraft  emissions  and  in  the 
control  of  aircraft  engine  emissions,  though  this  latter  topic  is  beyond  the 
scope  of  this  report. 

On  March  24,  1978,  a  Notice  of  Proposed  Rule  Making  (NPRM)  was  published 
in  the  Federal  Register  to  announce  the  intention  of  the  EPA  to  amend  the  1973 
engine  emission  regulation.  Included  in  the  NPRM  was  the  establishment  of  an 
FAA/EPA  air  quality  study  to  relate  aircraft  emissions  to  ambient  air  quality. 

This  report  presents  the  results  of  this  FAA/EPA  air  quality  study  to 
assess  the  impact  of  aircraft  emissions  of  CO,  HC,  and  N0X  at  populated 
locations  in  the  vicinity  of  airports.  This  assessment  includes  the  results 
of  the  recent  modeling  and  monitoring  efforts  at  Washington  National  (DCA) , 

Los  Angeles  International  (LAX),  Dulles  International  (IAD),  and  Lakeland 
Florida  airports  and  a  modeling  update  of  aircraft  generated  pollution  at  Los 
Angeles  International  (LAX),  John  F.  Kennedy  (JFK),  and  Chicago  O'Hare  (ORD) 
airports,  as  pollution  at  these  three  airports  was  cited  in  the  1973  engine 
emission  standards. 

ihe  study  methodology  is  indicated  in  Figure  1.1.  Measured  air  quality 
concentrations  associated  with  the  operation  of  commercial  aircraft  were  both 
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extrapolated  to  worst  case  conditions,  as  implied  by  the  NAAQS,  and  also  used  to 
verify  oub  models.  Features  of  these  sub  models  were  subsequently  incorporated 
into  a  major  airport  model  which  was  then  used  to  reassess  the  impacc  of 
aircraft  emissions  on  air  quality  at  JFK,  OPD,  and  LAX.  In  addition,  the 
results  of  the  modeling  and  monitoring  of  emissions  from  general  aviation 
airports  and  the  results  of  an  airport  model  validation  effort  conducted  at 
Williams  Air  Force  Base  are  reported. 

This  report  first  presents  a  brief  discussion  of  the  history  of  air 
pollution  measurements  and  the  air  quality  models  used  in  the  EPA's  most 
recent  assessment  of  aircraft  impacts  as  a  basis  for  the  NPRM .  This  discussion 
highlights  the  discrepancies  between  studies  which  have  led  to  the  present 
uncertainties  about  overall  impacts  of  aircraft  operations  and  the  need  for 
additional  engine  emission  standards.  Subsequent  sections  of  the  report 
discuss  the  recent  monitoring  and  modeling  studies  which  have  been  identified 
above,  and  the  relationship  of  their  results  to  these  issues.  It  will  be  seen 
that  one  of  the  most  difficult  questions  addressed,  but  not  yet  completely 
answered  by  these  studies,  is  how  significant  the  effect  of  enhanced  initial 
dispersion  of  aircraft  plumes  is  upon  concentrations  of  pollutant  reaching 
populated  or  public  areas  near  an  airpor*  .  Each  of  the  recent  studies  sheds 
light  on  a  facet  of  this  problem,  and  leads  to  the  development  or  refinements 
of  mode  specific,  sub-model  components  of  the  overall  airport  model  via  the 
methodology  indicated  in  Figure  1.1.  Review  of  the  previous  modeling  efforts, 
as  well  as  the  resuits  presented  in  this  report,  shows  that  the  significance 
of  these  improvements  is  sometimes  a  matter  of  perspective.  That  is,  if 
concentrations  are  predicted  to  be  well  below  applicable  standards  by  conserva¬ 
tive  screening  models,  identification  of  additional  dilution  factors  is 
relatively  unimportant.  An  underestimate  of  the  source  inventory,  on  the 
other  hand,  can  lead  to  associated  distrust  of  the  capabilities  of  a  pollutant 
dispersion  model .  Much  of  this  volume  and  the  supporting  Airport  Vicinity  Air 
Pollution  (AVAP)  model  calculations  included  in  Volume  II  are,  therefore, 
concerned  with  identifying  worst  case  meteorological  and  aircraft  operational 
conditions  leading  to  the  highest  emission  rates  and  concentrations  of  problem 
pollutants.  In  these  cases,  the  accuracy  and  precision  of  modeling  results 
supported  by  measurement  data,  can  become  a  critical  issue  in  defining  the 
needs  for  additional  standards. 
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Volume  II  contains  the  results  of  a  detailed  application  of  the  updated 
AVAP  model  to  LAX,  JFK  and  ORD  airports.  The  model  incorporates  the  most 
recent  submodels  and  sub-model  parameters  obtainable  from  the  airport  monitor¬ 
ing  programs  described  herein,  cs  well  as  the  most  recent  updates  to  the 
theory  of  dispersion  of  atmospheric  pollutants.  The  principal  findings  of 
Volume  li  have  also  been  incorporated  into  this  Summary. 

1.2  OBJECTIVES 

During  the  past  several  years  a  number  of  government  agencies,  includ¬ 
ing  the  EPA,  FAA,  and  USAF  have  been  engaged  in  a  comprehensive  program  to 
assess  the  effects  of  aircraft  emissions  upon  air  quality.  While  the  motiva¬ 
tion  to  evaluate  such  impacts  originated  with  the  1970  Amendments  to  the  Clean 
Air  Act,  the  March  24,  1978  Notice  of  Proposed  Rule  Making  (NPRM) ,  announcing 
the  EPA's  intention  to  modify  the  i973  engine  emission  standards,  provided  a 
clear  mandate  to 

•  resolve  the  ambiguities  of  previous  monitoring  and  model¬ 
ing  efforts, 

•  update  airport  dispersion  modeling  assessments  to  icflect 
recent  modeling  improvements,  and 

•  measure  pollutant  levels  near  aircraft  in  a  manner  that 
would  clearly  determine  aircraft  emissions  impact, 

so  that  realistic  engine  emissions  standards  could  be  established  on  the  basis 
of  the  best  available  information. 

While  the  objective  as  stated  above  suggests  a  program  of  all-encompass¬ 
ing  scope,  it  is  useful  at  the  outset  to  consider  the  limitations  of  this 
endeavor.  The  project  concerns  itself  solely  with  the  ground- level ,  air 
quality  impacts  of  aircraft  exhaust 

•  on  or  in  the  near  vicinity  of  airports, 

•  in  areas  of  possible  public  exposure,  and 

•  relative  to  existing  or  potential  National 
Ambient  Air  Quality  Standards  (NAAQS) . 

Thus,  for  example,  this  report  is  not  concerned  with  stratospheric  impacts, 
the  hydrocarbon  odor  nuisance  problem  at  airports,  the  combined  effects  of 
pollution  from  aircraft,  access  vehicles,  and  service  vehicles,  or  the  level 
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of  pollution  inside  the  passenger  te...iinal,  despite  the  fact  that  such  consider¬ 
ation0  may  have  explicit  or  implicit  effects  on  the  determination  of  adequate 
aircraft  emission  standards. 

1.3  APPROACH 

The  principal  strategy  was  to  assess  the  air  quality  impacts  of  aircraft 
exhaust  through  monitoring  of  aircraft  pollution  impacts  within  0.5  km  of  the 
aircraft.  Tills  served  the  dual  purpose  of  supplying  actual  men  ured  impacts 
with  which  one  could  infer  average  and  worst  case+  pollutant  concentrations 
and  of  providing  a  research  grade  data  base  with  which  one  could  investigate 
and  parameterize  the  aircraft  plane  dispersion  physics  in  or.ier  to  improve  the 
predictive  accuracy  of  sub-models  within  an  airport  model,  and  hence,  ultimate¬ 
ly  improve  the  predictive  power  of  airport  air  quality  assessment  models. 

With  one  such  improved  model,  AVAP,  -t  was  then  possible  to  simulate  worst 
case  pollutant  conditions  at  major  S.  airports:  an  objective  that  would 
have  been  unacceptably  expensive  to  attain  solely  through  ambient  air  monitor¬ 
ing  programs. 

The  monitoring  programs,  that  provided  the  basis  for  pursuit  of  the 
aircraft  exhaust  impact  assessment  objective  through  the  above-described 
strategy,  are  summarized  in  Table  1.1.  These  experiments,  described  in  this 
report  and  in  other  indicated  documentation,  share  two  important  character¬ 
istics  that  set  them  apart  from  previous  monitoring  programs.  First,  they 
were  designed  to  focus  in  on  specific  aircraft  modes  of  operation.  Hie 
orientation  of  receptors  at  Washington  National  (DCA),  seen  in  Figure  1.2, 
provides  an  example  of  such  3  modal  locus.  Under  winds  from  the  NW~N  direc¬ 
tions,  the  pollution  cloud  from  queueing  aircraft  is  transported  across  the 
network  of  monitors  while  under  winds  from  the  NNF.-SK  directions,  the  plumes 
created  by  the  high-thrust  takeoff  mode  are  sampled  by  the  same  monitors. 

Second,  these  experiments  achieve  a  separation  or  aircraft  pollution  from 
other  source  related  pollution  via  a  mul t i-s t at i on  receptor  array  either 
operating  in  a  low-background  environment,  such  as  at  Williams  AFB,  Arizona 
(Yamartino  et  al.,  1980)  or  else  campling  at  a  sufficiently  high  rate  to 

+l,worst  case"  is  generally  taken  here  to  indicate  the  highest  hourly  average 
concentration  per  annum  or  the  meteorological  and  aircraft  operations  condi¬ 
tions  leading  to  such  concentrations. 


Table  1.1.  Summary  of  Recent  Airport  Air  Quality  Experiment 


7 


MONITORING  LOCATIONS  AT  DCA 
FOCUSED  ON  QUEUEING  AND  TAKEOFF  MODES 


NORTH 


Fig.  1.2.  Monitoring  Locations  at  DCA.  With  the  exception 
of  site  2,  positioned  to  monitor  ambient  back¬ 
ground  concentrations,  the  monitors  were  located 
to  measure  the  impacts,  as  well  as  characterize 
the  advection  and  dispersion,  of  pollution  origi¬ 
nating  from  the  queuing  mode  under  winds  from  the 
NNW  thru  NNE  directions  and  from  the  takeoff  mode 
under  winds  from  the  N  thru  E  directions. 


RAPID  SAMPLING  ENABLES  CONFIDENT  SEPARATION 
OF  AIRCRAFT  PRODUCED  AND  BACKGROUND  SIGNALS 
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facilitate  such  a  separation.  The  time  history  of  N0X  pollution  seen  at 
two  of  the  DCA  monitors  during  a  one  hour  period  (Figure  1.3)  illustrates 
the  analysis  possibilities  created  through  thxs  high  sampling  rate  approach. 
Takeoff  plume  pollution,  easily  identified  by  the  N0X  rise  above  background 
levels,  are  seen  to  disperse  (i.e.  reduce  in  peak  concentration  while  spreading 
out  in  time)  as  they  move  from  the  nearby  station  1  to  the  more  distant  station 
3.  The  actual  impact  of  these  plumes  may  be  estimated  via  a  background  sub¬ 
traction  (i.e.  as  seen  jv  the  dotted  line)  or  through  estimation  of  the  summa¬ 
tion  of  the  areas  un  •  ;r  the  major,  aircraft  related  pulses  (i.e.  as  indicated 
by  the  shaded  areas)  via  a  technique  requiring  a  threshold  concentration  above 
background  to  De  reac.  ed  before  inclusion  of  a  particular  pulse.* 

In  addition  to  the  air  monitoring  programs  considered  above,  improve¬ 
ments  to  the  AVAP  model  and  its  predictive  power  resulted  from  a  number  of 
efforts,  detailed  in  Volume  II,  that  include: 

•  use  of  updated  engine  emission  factors 

•  code  modification  to  permit  direct  input  of  the  hourly 
number  of  departures  (previously  computed  on  the  basis  of 
arrivals 

•  use  of  airport  specific  t imes-in-mode  for  idle,  taxi, 
and  queueing  that  were  observed  during  a  three  airport 
field  study.  These  times  often  differed  substantially 
from  airport  averaged,  AP-42  LTO  cycle  times. 

One  of  the  more  important  findings  of  this  AVAP  modification  and  multi-airport 
analysis  was  that  for  worst  case,  "hot-spot"  pollution  assessments  one  need 
only  consider  the  pollution  from  a  single  runway  complex  and  the  queueing  area 
adjacent  to  it.+  These  runway  complexes  are  generally  sufficiently  well 
separated  from  one  another  that  superposition  of  pollutants  at  the  "hot-spot" 
from  other  complexes  is  negligible  to  first  order.  This  observation  should 
eliminate  many  needlessly  repetitive  analyses  of  similarly  configured  airports. 
It  should  further  be  noted  that  the  USAF/USN  sponsored,  validation  and  accuracy 
evaluation  program  for  the  Air  Quality  Assessment  Model  (AQAM) ,  provided 
accuracy  measures  for  AQAM  and  its  nearly  identical  civilian  counterpart  AVAP, 
that  were  previously  only  conjecture. 

ities  will  subsequently  be  referred  to  as  "background  subtracted" 
integrated"  concentrations. 

teristic  subsequently  referred  to  as  "runway  complex  factoriza- 
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1.4  RESULTS 

Rather  than  to  proceed  from  experiment  to  experiment,  as  is  done  in 
the  text  of  the  report,  let  us  consider  the  various  engine  emitted  pollutants 
and  the  information  regarding  those  emissions'  impacts  as  determined  by  these 
moni.toring  and  modeling  exercises. 

1.4.1  Carbon  Monoxide 

Experimental  and  modeling  efforts  of  the  early  70' s  indicated  that 
violations  of  the  NAAQS  one-hour  CO  standard  of  35  ppm  were  indeed  possible  at 
airports.  Measured  and  modeled  peak  hourly  levels  of  46  ppm  and  24  ppm 
respectively  at  LAX  (Platt  et  al . ,  1971),  for  example,  suggested  that  aircraft 
emissions  were  a  serious  problem  relative  to  the  one-hour  standard;  however,  a 
number  of  factors  contributed  to  this  misleading  implication  of  aircraft 
inc lud ing : 

•  aircraft  source  characteristics  (i.e.,  initial  plume 
volume  and  rise)  were  not  understood  and  thus  not  modeled. 

Modeling  of  aircraft  emissions  neglecting  initial  dispersion 
can  lead  to  arbitrarily  high  concentrations  depending  on 
source  location. 

•  background  concentrations  were  often  not  measured,  making  it 
difficult  to  isolate  the  aircraft  or  airport  contributed  con¬ 
centrations  from  those  of  the  surrounding  region. 

•  building  wake  effects  can  greatly  magnify  the  impact 

of  the  multitude  of  CO  sources  around  the  terminal.  The 
modeling  of  such  enhancing  effects  was  (and  still  generally 
is  because  of  the  complexity)  ignored. 

•  emissions  from  other  sources,  particularly  service  and 
access  vehicles  may  dominate  aircraft  sources  in  the 
vicinity  of  the  terminal  where  the  highest  concentrations 
were  observed . 

•  peak  observed  concentrations  were  underpredicted  by  a 
factor  of  2-3. 

This  latter  consideration  of  2-3-fold  model  underprediction,  coupled  with  the 
other  uncertainties,  particularly  the  unknown  and  unmodeled  characteristics  of 
the  aircraft  plume,  certainly  invited  the  speculation  that  aircraft  were 
responsible  for  the  modeling  deficit  and  thus  were  the  principal  source. 

The  CO  monitoring  experiment  along  the  main  taxiway  at  Dulles  Inter¬ 
national  (Smith  et  al,  1977)  was  the  first  to  isolate  the  impact  of  aircraft 
emissions  alone  and  indicated  that  the  initial  turbulent  mixing  caused  by  the 
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plume  plus  the  subsequent  plume  rise  and  accompanying  enhanced  vertical 
dispersion  lead  to  rather  low  CO  concentrations  near  the  aircraft  (excepting, 
of  course,  directly  behind  the  engines) .  While  peak  instantaneous  CO  levels 
reached  10  ppm  at  the  first  monitor,  located  only  65m  from  the  taxiway 
centerline,  the  maximum  impact  of  a  single  aircraft  to  the  hourly  average 
concentration  remained  below  0.06  ppm  at  this  distance.  Thus,  even  adjacent 
to  a  busy  taxiway,  it  would  be  unlikely  for  hourly  average  CO  concentrations 
to  exceed  several  ppm. 

Though  substantially  higher  concentrations  were  anticipated  downwind 
cf  a  group  of  queueing  aircraft,  the  recent  5-day  experiment  at  Los  Angeles 
(Section  6)  failed  to  indicate  peak  hourly  CO  in  excess  of  3  ppm  at  a  distance 
of  200-300  meters  from  the  closest  of  the  queueing  aircraft. 

Perhaps  the  most  convincing  evidence  concerning  peak  hourly  CO  from 
taxiing  and  queueing  aircraft  comes  from  consideration  of  the  cumulative 
frequency  distributions  (CFDs)  of  CO  levels  observed  during  the  month-long  DCA 
experiment.  Keeping  in  mind  the  large  error  which  is  possible  when  estimating 
the  worst  case  once  per  year  on  the  basis  of  a  one  month  data  sample,  Figure 
1.4  indicates  that,  under  the  assumption  of  log-normality  of  the  hourly  CO 
frequency  distribution,  maximum  hourly  per  annum  CO  levels  of  10  ppm  may  be 
expected  300m  downwind  of  the  aircraft  queue,  with  a  probable  maximum  of  ^5  ppm 
due  to  aircraft  alone.  Similar  analyses  applied  to  data  from  station  1, 
located  only  100-150m  from  the  queueing  aircraft,  suggests  an  aircraft  attribu¬ 
table  maximum  hourly  level  of  10  ppm;  however,  it  is  unlikely  that  the 
general  public  would  be  sufficiently  close  to  queueing  aircraft  to  receive 
that  CO  dose. 

Finally,  we  note  that  the  AVAP  model,  adjusted  to  include  the  plume 
dynamical  effects  measured  during  the  Dulles  experiments,  predicts  "worst 
case"  CO  concentrations  of  7  ppm  100-300  meters  downwind  of  queueing  aircraft 
at  LAX.  Figure  1.5  shows  the  isopleths  resulting  from  the  calculation,  and 
one  observes  that  levels  of  public  exposure  exceeding  6  ppm  would  be  unlikely 
during  this  "worst  case"  scenario. 

One  notes  that  the  combined  effects  of  the  phased  daily  cycles  of 
aircraft  operations  and  meteorological  parameters  related  to  dispersion,  and 
wind  direction  variability  make  violation  by  aircraft  exhaust  plumes  alone 
improbable.  The  single  violation  of  the  8-hcvur  standard  of  9  ppm  observed  at 
DCA  was  accompanied  by  simultaneous  violation  over  much  of  the  Washington  area. 


OBSERVED  CO  STATION  4 

BKGD. SUBTRACTED 
INTEGRATED  PULSE 
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1.4.2  Hydrocarbons 

Hydrocarbon  emissions  are  of  concern  as  the  presence  of  reactive 
hydrocarbon  species  are  conductive  to  the  subsequent  formation  of  ozone. 
Motivation  for  control  of  aircraft  emitted  HC  results  partly  from  estimates 
that  aircraft  account  for  1-3%  of  the  total  HC  emitted  in  an  Air  Quality 
Control  Region  (AQCR)  on  an  annual  basis;  thus,  while  aircraft  are  not  a 
dominant  source,  they  reDresent  a  significant  source  for  control  as  they  are 
comparable  with  many  other  source  categories. 

Measurements  at  Dulles  and  AVAP  modeling  agree  that  total  HC  concentra¬ 
tions,  expressed  as  ppm  equivalent  methane  (CH^),  correspond  well  in  space, 
time,  and  magnitude  with  CO  levels  associated  with  these  commercial  aircraft. 
This  is  not  particularly  surprising  since  hydrocarbons  and  CO  are  both  emitted 
during  the  same,  low  power  setting,  aircraft  operational  modes.  Figure  1.6 
shows  the  peak  hourly  total  HC  (THC)  contour  resulting  from  a  "worst  case" 
modeling  of  LAX  during  the  8-9AM  period.  Though  restrictive  considerations  of 
the  reactive  (RHC)  component  and  the  three-hour  average  would  act  to  somewhat 
reduce  the  area  of  this  contour,  it  is  still  anticipated  that  the  0.25  ppm 
contour  covers  an  area  several  times  the  airport  size. 

1.4.3  Oxides  of  Nitrogen 

The  issue  of  oxides  of  nitrogen  (N0X)  impacts  created  by  aircraft  is, 
as  with  CO,  a  localized  "hot-spot"  problem  related  to  existing  and  possible 
additional  NAAQS  for  nitrogen  dioxide  (N02),  and  has  been  addressed  primarily 
through  the  DCA  monitoring  program.  Unfortunately  the  issue  is  further 
complicated  by  the  fact  that 

•  present  and  possible  future  NAAQS  standards  pertain  to 
NO2  levels  and  not  N0X  levels.  (N0X  ^  NO  +  NO2) 

•  there  is  presently  only  an  annual  average  NAAQS  of  0.05 
ppm  NO2  though  a  peak  hourly  standard  in  the  range 
0.2-0. 5  ppm  is  currently  oeing  reviewed  by  the  EPA 

•  plume  dispersion,  while  reducing  the  concentration  of 
inert  species,  will  entrain  mere  amibient  oxidant  result¬ 
ing  in  further  conversion  of  engine  emitted  NO  to  NO2 ; 
thus,  NO2  levels  will  peak  at  some  distance  downwind  of 
the  aircraft 

•  the  peak  NO2  attributable  to  aircraft  is  a  function  of 
existing  ambient  levels  of  NO,  NO2,  O3,  and  sunlight. 
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In  order  to  determine  the  impact  of  aircraft  emissions  on  annual 
average  air  quality  for  comparison  with  the  annual  NO2  NAAQS,  it  was  necessary 
to  regress  measured  NO2  concentration  levels  against  airplane  activity. 

Figure  1.7  shows  this  relationship  for  the  pulse  integrated  NO2  levels  at 
Station  1.  The  statistical  significance  of  the  slope  of  the  regression 
coupled  with  a  regression  !,y"  intercept  consistent  with  zero  enables  or  to 
confidently  estimate  the  annual  average  aircraft  impact.  The  projected  annual 
average  aircraft  impact  of  0.005  ppm  (5  ppb)  is  small  compared  with  the  0.05 
ppm  NAAQS . 

This  paucity  of  data  results  in  greater  uncertainties  for  estimation  of 
maximum  hourly  average  N0X  levels.  Figure  1.8  indicates  that,  assuming  log- 
normality  of  the  hourly  NOx  cumulative  frequency  distribution  (CFD),  maximum 
hourly  per  annum  N0X  levels  of  0.2  to  0.5  ppm  due  to  aircraft  operations  alone 
may  be  expected  several  hundred  meters  downwind  of  the  location  where  aircraft 
begin  their  takeoff  roll.  Depending  on  the  oxidation  rate  of  the  aircraft 
emitted  NO  into  NO2,  NO2  levels  is  excess  of  0.2  ppm  may  materialize. 

AVAP  modeling  of  a  typical  busy  commercial  airport  under  worst  case 
activity  and  dispersion  conditions  indicates  (Figure  1.9)  N0X  levels  exceeding 
0.5  ppm  more  than  one-half  mile  from  the  end  of  the  runway  complex,  but  the 
key  question  is  how  these  N0X  levels  translate  into  NO2  levels.  The  N02/N0X 
ratio  is  a  function  of  plume  dispersion  rate  and  transport  time,  sunlight 
intensity,  and  background  levels  of  NO,  N02,  and  O3  and  a  reactive  plume 
cal  culation  is  required  to  obtain  a  more  definitive  prediction;  however,  using 
simple  assumptions  regarding  the  amount  of  NO2  emitted  directly  by  the  aircraft, 
the  rate  of  NO  oxidation,  and  the  ambient  O3  level,  it  is  reasonable  to  expect 
several  tenths  of  pom  of  NO2  at  distances  of  possible  public  exposure.  This 
is  within  the  range  of  levels  under  consideration  by  the  EPA  as  a  possible 
short  term  NAAQS  for  NO2 . 

1.5  CONCLUSIONS  AND  RECOMMENDATIONS 

Recent  airport  air  quality  monitoring  studies  at  four  airports  suggest 
maximum  hourly  average  CO  concentrations  of  5  ppm  in  areas  of  expected  public 
exposure.  These  measurements  and  estimates  based  on  extrapolation  of  measured 
results  to  probabilities  corresponding  to  one  hour  per  year  suggest  small 
liklihood  of  violating  the  35  ppm  NAAQS. 


SLOPE 


Station  1  Average  Pulse  Integrated  NO2  Concentration  Versus  Departure  Rate  at  DCA. 
The  slope  of  the  regression  line  (solid  line)  times  the  average  hourly  departure 
rate  of  10.1  aircraft  indicates  an  aircraft  impact  of  5  ppb  on  the  annua?,  average 
NO?  nearby  the  runway.  The  annual  average  NAAQS  for  NO2  is  50  ppb. 
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Updated  AVAP  modeling  studies  of  the  major  airports  under  assumed  worst 
case  meteorological  and  aircraft  operations  conditions  also  suggest  maximum 
hourly  CO  concentrations  of  ~5  ppm  at  distances  where  public  exposure  might 
occur . 

In  light  of  anticipated  CO  reductions  that  will  accompany  and  hydro¬ 
carbon  control,  the  air  quality  benefits  of  a  separate  emission  standard  for 
CO  are  highly  questionable. 

The  one  potentially  problematic  area  with  respect  to  CO,  and  not  dealt 
with  in  this  report,  is  at  and  immediately  around  the  terminal,  where  the 
combined  effect  of  other  CO  sources  and  the  effect  of  the  building  wake 
could  conspire  to  cause  violations;  however,  reduced  engine  running  time  in 
the  terminal  area  (resulting  primarily  from  a  desire  to  save  fuel)  has  greatly 
reduced  the  aircraft  role  in  such  an  event . 

Hydrocarbon  emissions  from  aircraft  engines  contribute  to  the  pervasive 
ozone  air  quality  problem  in  many  urban  areas.  While  the  overall  contribution 
by  aircraft  is  on  the  order  of  1%  to  3%  of  the  total  hydrocarbons,  ciu_  need  to 
control  these  emissions  can  be  justified  in  that  it  will  be  necessary  to 
control  many  other  small  and  diversely  located  sources  to  achieve  the  ozone 
NAAQS  in  areas  around  most  major  airports.  Further  control  impetus  is  gener¬ 
ated  by  the  desire  to  reduce  hydrocarbon  odors  around  airpots. 

Analyses  of  the  oxides  of  nitrogen  data  from  these  airport  studies 
suggest  that  aircraft  impact  on  local,  annual  average  NC>2  levels  will  be 
confined  to  10-20%  of  the  0.05  ppm  annual  average  NAAQS.  However,  the  possi¬ 
bility  of  aircraft  alone  causing  peak  hourly  NO2  concentrations  in  excess  of 
0.2  ppm  exists.  Similarly,  AVAP  modeling  suggests  possible  high  NO2  near 
the  ends  of  runways.  At  the  present  time  there  is  no  short-term  (e.g.,  hourly 
average)  NO2  standard;  however,  the  EPA  is  reviewing  the  need  to  establish 
such  a  standard.  The  suggested  range  of  such  a  standard  is  0.2  ppm  to  0.5  ppm 
hourly  average.  However,  because  the  actual  NO2  level  is  strongly  dependent 
upon  a  number  of  atmospheric  parameters,  incorporation  of  a  reactive  plume 
model  with  realistic  background  levels  of  NO,  NO2 ,  and  O3  is  needed  before 
modeling  can  yield  a  more  definitive  prediction  of  peak  hourly  NO2  concentra- 
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Just  as  the  issue  of  NO2  impact  assessment  is  more  complicated  than  for 
CO  or  HC,  so  too  is  the  issue  of  N0X  control.  Plagued  by  poor  control  tech¬ 
nology  and  high  control  costs,  NOx,  which  originates  primarily  from  the  high 
thrust  takeoff  mode,  cannot  be  "managed"  as  effectively  through  minimization 
of  engine  idle  time  as  can  be  CO  and  HC. 
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2  HISTORICAL  SURVEY  OF  AIRPORT  AIR  QUALITY  STUDIES 

As  consequences  of  the  1967  Air  Quality  Act  (U.S.  Congress,  1967) 
and  the  1970  Clean  Air  Act  Amendments  (U.S.  Congress,  1970)  a  number  of  studies 
have  been  conducted  to  determine  the  contributions  of  aircraft  emissions  to 
the  air  quality  in  the  vicinity  of  airports.  This  chapter  summarizes  the 
history  of  these  studies,  their  purposes,  and  the  conclusions  which  have  been 
previously  drawn  from  them.  Also  pointed  out  is  a  number  of  study  difficulties 
which  have  brought  those  study  conclusions  into  question,  and  have  led  to 
additional  studies  aimed  at  resolving  remaining  questions.  The  final  section 
of  this  chapter  discusses  the  meteorological  aspects  of  "worst  case"  air 
quality  conditions. 

2.1  HISTORICAL  OVERVIEW 

The  contributions  of  aircraft  as  sources  of  air  pollution  were  not 
seriously  considered  until  the  introduction  of  turbojet  aircraft  into  air 
carrier  service  in  the  late  1950's.  Even  though  the  particulate  emissions  of 
those  earlier  engines  were  highly  visible,  the  first  two  reviews  of  their 
potential  contributions  by  the  Los  Angeles  Air  Pollution  Control  District 
(George  and  Burlin,  1960)  and  the  Coordinating  Research  Council  (1960)  did  not 
consider  the  total  emissions  significant  enough  to  warrant  further  investiga¬ 
tion. 

It  was  not  until  a  second  study  by  the  LAAPCD  (Lemke  et  al,  1965) 
and  the  Report  of  the  Secretary  of  HEW  on  the  "Nature  and  Control  of  Aircraft 
Engine  Exhaust  Emissions"  (U.S.  DHEW ,  1968)  that  the  subject  of  control  was 
brought  into  serious  discussion  as  "feasible  and  desirable."  Of  principal 
concern  were  CO  and  organic  particulates.  These  findings  by  the  HEW  led  to 
further  quantification  of  emissions  inventories  by  Northern  Research  and 
Engineering  Corporation  (NREC)  under  contract  to  the  Public  Health  Service  of 
HEW  and  later  to  the  U.S.  EPA  (e.g.,  Bastress  et  al ,  1971).  This  was  followed 
by  the  development  of  an  air  quality  impact  assessment  technique  by  the 
Northern  Research  and  Engineering  Corporation  (NREC)  with  assistance  from 
Environmental  Research  &  Technology,  Inc.  (ERT).  The  result  was  the  first  air 
quality  model  specifically  designed  for  airports. 
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A  major  study  of  the  potential  impact  of  aircraft  emissions  upon 
air  quality  at  six  airports,  including  five  of  the  busiest,  was  then  launched 
by  the  EPA  using  the  NREC  methodology  (Platt  et  al,  1971).  It  required 
compilation  of  emission  inventories  for  Washington  National  (DCA),  Chicago's 
O'Hare  (ORD) ,  New  York's  John  F.  Kennedy  (JFK),  and  Los  Angeles  International 
(LAX),  as  well  as  for  the  new  Tamiami,  Florida  and  Van  Nuys  California  airports 
Since  aircraft  emissions  data  were  lacking,  the  result  was  a  good  initial  set 
of  relative  estimates  but  with  uncertain  overall  accuracy.  The  dispersion 
model  was  relatively  simple  in  its  derivation,  but  the  validity  for  its  use 
with  aircraft  sources  hadn't  been  fully  tested. 

At  about  this  time  (1971-72)  air  quality  measurement  data  were  being 
collected  in  the  vicinity  of  LAX  to  provide  a  comparative  data  base  in 
addition  to  a  direct  assessment  of  air  quality  in  the  airport  vicinity  (I 
1971).  The  results  of  these  LAAPCD  measurements  and  the  NREC  modeling  st.  y 
served  as  the  basis  tor  the  initial  assessment  by  the  EPA  (EPA,  1972)  of  the 
air  quality  impact  of  aircraft  engine  emissions.  Although  the  EPA  used 
updated  emissions  data,  the  comparisons  of  model  calculations  with  the  few 
measurements  available  yielded  several  large  discrepancies  including  some 
severe  underpredict  ions .  This  led  to  a  model  verification  study  at  DCA  by 
Geomet .  After  several  modifications  the  model  was  compared  with  a  new  set  of 
monitoring  data  from  DCA,  but  it  still  had  a  tendency  to  underpredict  changes 
in  concentrations  as  identified  by  Argonne  National  Laboratory  staff  (Wangen 
and  Conley ,  1975 )  . 

With  FAA  support  Argonne  developed  the  more  complex  Airport  Vicinity 
Air  Pollution  (AVAP)  model.  AVAP,  also  a  Gaussian  dispersion  and  transport 
model,  had  several  refinements  for  treating  accelerating  line  sources  as  well 
as  other  types  of  point,  area,  and  line  sources.  AVAP's  emission  inventory 
routine  was  also  more  sophisticated  and  tracked  patterns  of  aircraft  activities 
as  a  function  of  wind  direction  and  functional  mode  (Rote  et  al,  1973). 
Simultaneous  development  of  a  military  version,  the  Air  Quality  Assessment 
Model  (AQAM) ,  was  undertaken  by  Argonne  for  the  U.S.  Air  Force. 

The  comparison  between  AQAM  (adopted  for  civilian  aoplicafion)  and  the 
Geomet  model  indicated  that  AQAM  was  relatively  successful  in  describing  the 
measurements  made  at  DCA.  An  initial  validation  exercise  for  AVAP  was  perform¬ 
ed  at  O'Hare  Airport  where  especially  detailed  emissions  data  had  been  obtained 
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Due  to  problems  resulting  from  inadequate  determination  of  off-airport  sources 
of  CO,  the  ORD  validation  effort  (Rote  et  al,  1973)  was  more  successful 
when  N0X  measurements  were  compared  with  AVAP  predictions.  A  subsequent 
effort  to  compare  AVAP  results  with  monitoring  data  from  Hartsville  Internation¬ 
al  Airport  in  Atlanta  was  complicated  by  troubles  with  the  measurement  program 
(Cirillo,  et  al ,  1975).  Though  absolute  concentrations  of  CO,  N0X ,  and 
HC  remained  in  doubt,  the  Atlanta  study  results  have  been  useful  in  exploring 
the  relative  sensitivity  of  concentration  patterns  to  various  aircraft  opera¬ 
tional  procedures. 


2.2  SUMMARY  OF  PREVIOUS  AIRPORT  MODELING  AND  MONITORING  RESULTS 

There  is  continuing  concern  about  impacts  of  aircraft  related  emissions 
upon  air  quality  in  public  areas  both  inside  and  outside  of  airport  boundaries. 
On  a  regional  basis  HC  and  N0X  emissions  must  be  considered  because  of  their 
role  in  photochemical  oxidant  formation,  and,  not  surprisingly,  many  of  the 
cities  with  oxidant  standard  attainment  problems  have  large  airports  associated 
with  them.  An  airport's  contribution  in  these  major  urban  areas  usually 
constitutes  1-3%  of  the  emissions  burden  from  all  sources  of  HC  and  N0X . 

An  airport  at  a  rural  location  may  represent  the  largest  single  contribu¬ 
tor  on  an  annual  basis  to  the  inventories  of  CO,  HC ,  and  N0X  emissions  in  its 
area  of  air  quality  influence;  but,  it  is  local  effects  of  these  pollutants  in 
comparison  with  shorter  term  standards  which  are  currently  of  greatest  concern. 

The  continuing  questions  about  the  adequacy  of  airport  modeling  methods 
can  often  be  traced  to  one  of  the  following: 

1.  There  are  several  airport  models  which  give  widely 
vary  ing  resul ts . 

2.  The  modeling  assumptions  are  not  always  clearly 
defined;  and  their  applicability  to  specialized 
airport  source  geometries  (e.g.,  jet  engine  ex¬ 
hausts  in  terminal  areas)  has  been  questioned. 

3.  The  validation  experiments  for  airport  models  are 
few,  and  those  validations  that  do  exist  have  not 
been  particularly  successful  even  for  the  relatively 
sophisticated  models. 

A.  The  scales  of  interest  to  the  user  may  not  coincide  with 
those  for  which  the  mathematical  model  was  developed. 

A  user  may  often  expect  finer  resolution  of  concentration 
patterns  than  is  reasonable  or  better  agreement  between 
a  few  short-term  measurements  and  predicted  ensemble 
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average  concentrations  than  is  reasonable.  Similarly, 
computations  performed  using  an  excessively  crude  grid 
may  exaggerate  the  spatial  extent  of  highly  localized 
pollution  problems. 

5.  Analysis  of  a  "hot  spot"  problem  may  be  attempted 

with  a  model  that  does  not  contain  the  right  physical 
model  to  adequately  simulate  or  describe  it. 

Many  of  these  problems  result  as  often  from  misunderstanding  of  model  capabili¬ 
ties  or  user  intents  as  from  model  design  or  documentation  inadequacies. 

Often  the  relevant  question  to  ask  is  whether  the  problem  of  interest 
fits  the  assumptions  and  scale  of  resolution  of  the  "standard"  model  to  be 
selected.  If  not,  a  custom-designed  model,  or  more  often,  a  hybrid  f 
several  sub-models  or  complimentary  models  is  called  for.  This  new  model 
may  also  require  some  validation  via  comparison  with  measurements,  but  at 
least  it  should  be  constructed  to  represent  each  of  the  physical  processes  by 
a  sub-model  that  has  been  separately  validated  to  be  an  adequate  description 
of  that  process,  consistent  with  the  state  of  current  knowledge. 

There  are  relatively  few  major  airports  at  which  violations  of  air 
quality  standards  and  guidelines  have  been  measured  and  documented.  Most  of 
these  have  been  described  previously  in  the  review  of  airport  impacts  prepared 
by  the  EPA  (1972)  and  updated  recently  by  Lorang  (1978).  There  are  also  a 
larger  number  of  airports  at  which  violations  have  been  predicted  as  a  result 
of  air  quality  modeling  studies  though  supporting  measurements  are  unavailable. 
Controversy  often  arises  over  relative  contributions  of  emissions  from 
sources  related  to  aircraft,  nonaircraft  sources  at  airports,  and  those  due  to 
nonaircraft  sources  in  the  airport  environs.  A  number  of  studies  (Cirillo  et 
al  ,  1975;  Wang,  1975;  Lorang,  1978)  have  clearly  recognized  the  importance  of 
environmental  sources  of  CO,  N0X,  and  HC  in  that  their  dominance  often  inter¬ 
feres  with  the  model  validation  or  calibration  goals  of  the  study.  As  shown 
in  the  tables  given  in  the  subsections  below,  repeated  assessments  have  given 
widely  varying  results  at  the  same  airports.  The  result  has  been  that  the 
federal  EPA  and  several  regulatory  agencies  continue  to  raise  some  of  the 
same  questions  that  were  at  issue  before  these  modeling  studies  were  undertaken. 

In  each  of  the  study  results  reviewed  in  the  following  three  subsections, 
the  reliability  of  conclusions  relating  to  the  need  for  additional  controls 
was  thrown  into  some  doubt  because  there  was  not  convincing  agreement  between 
model  predictions  and  measured  concentrations.  Although  inventory  quality  or 
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measurement  data  quality  can  often  be  challenged,  the  lack  of  verification  for 
the  simpler  sub-components  of  the  model  calculations  has  made  it  difficult  to 
identify  with  any  uncertainty  the  main  reasons  for  poor  overall  comparisons. 

As  just  recently  pointed  out  by  Turner  (1979),  it  is  only  through  the 
repeated  verification  of  each  of  the  sub-components  of  an  air  quality  disper¬ 
sion  model,  that  the  validity  of  the  overall  model  results  will  be  eventually 
demonstrated.  The  last  section  considers  the  question  of  the  identification 
of  the  meteorological  conditions  that  are  associated  with  the  highest  ground 
level  concentrations  and  that  are  presumably  simulated  by  air  quality  models. 

2.2.1  Carbon  Monoxide  Studies 

The  results  of  the  principal  investigations  of  CO  concentrations  at 
major  airports  are  given  in  Table  2.1a.  The  results  of  these  studies  have  been 
reported  over  the  1971  to  1979  period  during  which  EPA  development  of  aircraft 
engine  control  strategies  has  been  continuing.  The  initial  modeling  estimates 
from  the  1971  NREC  study  (Platt  et  al,  1971)  for  LAX,  ORD,  DCA  and  JFK  airports 
are  presented  because  they  represent  comparable  modeling  assumptions  applied 
to  a  variety  of  airports.  As  can  be  seen  in  the  table  the  majority  of  the 
remaining  modeling  studies  utilized  AQAM  or  AVAP .  Careful  examination  of 
these  results  reveals  that  the  majority  of  maximum  measurements  observed  are 
within  a  factor  of  two  of  the  predicted  maximum  for  the  same  receptor  areas. 

The  model  estimates  of  maxima  are  also  generally  higher  even  though  the  model 
may  underpredict  the  majority  of  the  cases  which  result  in  moderate  concentra¬ 
tions.  Interpretation  of  CO  monitoring  and  modeling  data  is  not  complicated 
by  significant  chemical  reactivity  or  measurement  uncertainty.  However, 
precise  modeling  of  hot  spot  concentrations  adjacent  to  obstructions  requires 
modeling  the  flow  around  obstacles  which  is  beyond  the  capability  of  present 
Gaussian-type  airport  models.  Moreover,  the  locations  where  CO  violations  are 
often  suspected  have  a  large  contribution  of  CO  emissions  by  automobiles  and 
access/service  vehicles  as  well  as  aircraft. 

When  violations  of  the  8-hour  standard  are  encountered,  the  aircraft 
may  generally  be  identified  as  a  major  source  only  near  the  end  of  runways 
with  heavy  queuing  activity  during  worst-case  me t eor log ic a  1  conditions.  Since 
that  is  not  an  area  in  which  the  public  is  generally  exposed,  a  modeled 
violation  of  ambient  standards  downwind  of  a  queueing  area  probably  does  not 
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Table  2.1a.  CO  Measurement  and  Modeling  Study  Summary 
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present  a  health  hazard.  Lorang's  review  (1978)  claims  that  monitoring  has 
shown  that  CO  violations  have  occurred  in  terminal  areas  and  that  suggests 
that  aircraft  emissions  were  important  contributors  at  both  LAX  and  DCA.  It 
may  also  be  pointed  out,  however,  that  the  CO  levels  due  to  the  airport  (as 
modeled  in  Atlanta)  usually  drops  off  very  rapidly  with  distance.  Therefore 
potential  problems  are  localized  within  the  airport  property. 

The  recent  measurement  program  at  Boston's  Logan  Airport  (Smith  and 
Heinold,  1980)  illustrated  that  measured  concentrations  were  much  below 
standards  during  periods  of  high  airport  activity  in  areas  near  the  ends  of 
runways  with  long  queues  of  taxiing  aircraft.  The  highest  concentrations 
occurred  instead,  when  winds  from  nearby  urban  centers  coincided  with  a  strong 
nocturnal  inversion.  These  CO  concentrations  tend  to  be  overpredicted  by  most 
present  airport  dispersion  models.  Therefore,  it  appears  that  proper  modeling 
techniques  must  consider  these  situations  as  well  as  the  microscale  CO  problem 
in  terminal  areas.  The  latter  modeling  must,  however,  account  for  building 
wake  effects  and  local  sources,  such  as  vehicular  traffic,  if  the  true  relative 
impacts  of  aircraft  sources  are  to  be  realistically  portraved  at  the  terminal. 

2.2.2  Hydrocarbon  Studies 

The  results  summarized  in  Table  2.1b  are  for  the  same  airports  as  those 
given  in  Table  2.1a,  except  that  there  were  no  studies  of  HC  performed  at  the 
Seattle  or  the  Van  Nuys  airports.  It  is  immediately  apparent  that  both 
nonraethane  hydrocarbon  (NMHC)  and  total  hydrocarbon  (THC)  measurements  as  well 
as  predictions  are  well  above  the  160  yg/m^  (6  AM-9  AM  average)  established  as 
an  EPA  guideline  for  management  of  photochemical  pollutants.  Comparisons  of 
the  wide  ranging  concentrations  among  receptor  points  on  and  near  airports 
reveals  that  aircraft  do  indeed  contribute  to  the  elevated  values  in  the 
vicinity  of  airport  boundaries.  The  maximum  on-airport  concentrations  occur 
in  idling  and  taxiing  areas,  and  particularly  in  queues  awaiting  takeoff.  The 
studies  of  pollutant  control  strategies  at  Atlanta  (Cirillo,  et  al,  1975)  and 
the  recent  Boston  study  (Smith  and  Heinold,  1980)  both  indicated  that  regula¬ 
tion  of  queuing  and  taxiing  times  may  serve  as  effective  measures  for  diminish¬ 
ing  hydrocarbons  and  organic  particulates  (and  the  odors  associated  with 
these)  with  current  aircraft  engine  designs. 
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A  report  on  the  air  quality  associated  with  Air  Force  bases  (Daley  and 
Naugle,  1978)  and  (Naugie  et  al ,  1978)  suggests  that  HC  and  N0X  emissions  from 
aircraft  at  airports  present  the  greatest  potential  harm  according  to  the 
EPA's  pollution  standards  index  (PSI).  Since  present  and  projected  jet  engine 
designs  are  able  to  effectively  decrease  hydrocarbon  (and  CO)  emissions  by 
increasing  combustion  efficiency,  control  of  HC  and  CO  are  expected  to  be  less 
difficult  than  N0X.  Because  of  its  rural  location,  the  study  at  Williams  AFB 
avoided  the  problem  of  high  urban  background  pollution  conditions  for  model- 
measurement  comparisons.  Using  the  AQAM  model,  THC ' s  displayed  the  highest 
PSI  levels  at  distances  beyond  the  airport  boundary.  However,  this  PSI 
approach  for  HC  analysis  suffers  from  the  problems  inherent  in  using  simple 
guideline  HC  levels  as  measures  of  O3  production  and  oxidant  health  effects. 

Unfortunately,  most  studies  make  no  distinction  between  total  and 
reactive  hydrocarbons.  Even  when  conservative  assumptions  are  invoked,  a 
distinction  should  be  made  between  representing  N0X  (NO,  NO2)  and  NO2  and  THC 
and  NMHC.  In  an  oxidizing  atmosphere,  NO  is  converted  to  NO2 ,  whereas  CH4  is 
nonreactive  at  ambient  temperatures  and  ozone  concentrations.  To  acknowledge 
the  inconsistency  but  ignore  it  in  the  interpretation  of  monitoring  and 
modeling  studies  [as  in  Lorang  (1978)],  leads  to  excessively  pessimistic 
predictions  about  the  role  of  airports  in  violations  of  the  Air  Quality 
Guideline. 

2.2.3  Oxides  of  Niiroger  Studies 

The  results  of  N0X  aid  NO2  measurements  and  N0X  model  predictions  are 
presented  in  Table  2.1c  for  che  same  airports  (and  studies)  for  which  HC 
results  were  given  in  Table  2.1b.  It  should  be  noted  that  the  values  given  in 
the  NREC  modeling  study  relate  to  an  annual  average  standard.  For  most  of  the 
studies  involving  both  measurements  and  modeling,  hourly  values  are  given  in 
both  instances.  In  addition,  most  measurement  studies  report  both  NO2  and 
total  N0X ,  even  though  modeling  generally  assumes  that  N0X  is  the  more 
reliable  parameter  to  predict.  (This  is  especially  true  for  long  term  average 
predict  ions ) . 

Considering  an  annual  standard  of  0.05  ppm  (100  pg/m^)  for  NO2 ,  the 
conservative  assumption  that  NO2  =  N0X  concentrations  leads  to  the  conclusion 
that  most  of  the  airports  modeled  by  NREC  might  have  a  problem  meeting  the 
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Table  2. Id.  Reference  Key  for  Measurement 
and  Modeling  Study  Summary 


Symbol  Reference 


a  Platt  et  al  . ,  1971 

b  Cirilloet  al.,  1975 

c  Rote  etal.,  1973 

d  Wangen  and  Conley,  1975 

e  Smith  et  al .  ,  1977 

f  Shelar,  1978 

g  Smith  and  Heinold,  1980 

h  Yamartino  et  al . ,  1980 

i  Greenberg,  1978 

j  Schewe  et  al . ,  1978 


annual  standard;  however,  such  an  assumption  may  be  grossly  incorrect  as 
aircraft  emissions  are  predominantly  NO. 

Modeling  of  the  nitrogen  dioxide  (NO2)  air  quality  impact  resulting 
from  airport  operations  has  traditionally  been  based  on  the  use  of  the  standard 
Gaussian  dispersion  equation,  assuming  that  100%  of  all  N0X  emissions  are 
converted  to  NO2  (e.g.,  EPA  1972,  Cirillo  et  al ,  1975).  The  '00%  NO2 
conversion  assumption  is  used  because  the  Gaussian  dispersion  equation 
cannot  accurately  simulate  the  complex  photochemical  reactions  that  govern 
actual  NO2  concentrations.  This  assumption  results  in  an  overpt ed ict ion  of 
NO2  concentrations;  however,  it  is  useful  in  identifying  potential  problem 
areas  deserving  of  more  detailed  modeling  and/or  comprehens ive  monitoring 
ef  forts . 

It  has  been  pointed  out  by  Jordan  (1977a)  that  although  this  conservative 
convention  for  interpreting  long-term  averges  may  often  be  acceptable,  more 
serious  questions  arise  when  the  issue  of  a  short-term  (1  to  3  hr)  NO2  standard 
is  considered.  Jordan  has  examined  both  commercial  airports  (1977a)  and 
general  aviation  aircraft  (1977b)  and  presents  evidence  that  the  latter  are 
not  likely  to  contribute  significantly  to  either  long-term  or  short-term  NO2 
concentrations.  The  magnitude  of  the  maximum  1-hour  NO2  measurements  at  the 
ai.rports  listed  in  Table  2.1c,  however,  suggests  that  some  commercial  airports 
could  have  a  problem  in  the  near  future  if  a  stringent  hourly  standard  is 
established.  Of  principal  concern  is  the  total  NO2  value  at  distances  of  1 
to  5  km  from  the  runway  areas  associated  with  the  maximum  emission  rates  for 
NO.  The  rate  of  conversion  of  NO  to  NO2  in  these  areas  seems  to  be  a 
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critical  factor  in  determining  the  need  for  further  action  in  controlling 
aircraft  emissions  of  N0X . 

The  actual  conversion  of  aircraft  NOx  emissions  to  NO2  is  a  complex 
function  of  meteorology,  atmospheric  photochemistry,  and  ambient  concentra¬ 
tions  of  N0X,  ozone,  and  hydrocarbons.  N0X  emissions  from  aircraft  mainly 
consist  of  nitric  oxide  (NO).  For  example,  emission  measurements  from 
Pratt  and  Whitney  JT3D,  JT8D,  and  JT9D  jet  engines  have  shown  a  typical  N02/N0X 
emissions  ratio  of  4  to  8%  by  volume  (Pratt  and  Whitney,  1972).  This  is 
reflected  in  ambient  air  monitoring  measurements  at  airports,  where  the 
N02/N0x  ambient  ratio  was  found  to  be  lower  on  the  airport  grounds  than  in 
areas  surrounding  the  airport  (Lorang,  1978). 

A  qualitative  assessment  has  been  made  of  the  influence  of  aircraft  N0X 
and  hydrocarbon  emissions  on  ozone  formation  downwind  (Whitten  and  Hogo, 

1976)  .  The  conclusion  was  that  the  mixing  of  aircraft  jet  exhaust  with 
automobile  exhaust  can  cause  a  more  favorable  hydrocarbon/NOx  ratio  for  ozone 
formation  than  automobile  exhaust  alone.  A  simple  semi-quant i t at ive  treatment 
of  NO  to  NO2  conversion  at  airports  considered  only  one  main  chemical  reaction 
(Jordan  and  Broderick,  1978,  Jordan  and  Broderick,  1979);  this  treatment  is 
valid  only  over  short  transport  time  scales  where  the  presence  of  hydrocarbons 
can  be  neglected. 

In  order  to  quantitatively  predict  the  NO2  conversion  of  aircraft 
N0X  emissions  and  the  effect  of  aircraft  N0X  and  hydrocarbon  emissions  on 
downwind  ozone  concentrations,  a  sophisticated  photochemical  air  quality 
simulation  model  may  be  necessary.  A  number  of  photochemical  models  have  been 
developed  which  can  simulate  chemistry,  emissions,  and  atmospheric  transport 
processes  with  detailed  spatial  and  temporal  resolution.*  A  methodology  was 
developed  to  integrate  an  early  photochemical  model  (NEXUS/P)  with  airport 
land  use  development  (Norco  et  al.,  1973);  however,  this  photochemical 
model  used  a  chemical  kinetics  mechanism  that  is  now  obsolete. 

Only  one  study  has  been  conducted  that  has  used  a  detailed  photochemic¬ 
al  air  quality  simulation  model  to  examine  the  effect  of  N0X  and  hydrocarbon 

*The  discussion  here  is  limited  to  photochemical  models  applicable  in  the 
ui ban  troposphere.  The  impact  of  aircraft  emissions  aloft  on  the  strato¬ 
sphere  ozone  layer  requires  the  use  of  very  different  photochemical  model¬ 
ing  techniques  (Oliver  et.al.,  197/). 
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emissions  from  airport  operations  on  air  quality  in  the  vicinity.  This 
study  (Duewer  and  Walton,  1978)  was  done  in  the  San  Francisco  Bay  area  using 
the  LIRAQ-2  grid-based  photochemical  model.  The  modeling  showed  that  doubling 
airport  emissions  reduced  ozone  concentrations  slightly  at  San  Francisco 
Airport,  but  increased  ozone  downwind  by  approximately  0.003  ppm. 

However,  a  grid-based  photochemical  model  such  as  LIRAQ-2  is  very 
expensive  to  run,  both  in  terms  of  manpower  and  computer  time,  because  concen¬ 
trations  must  be  calculated  at  a  large  number  of  grid  cell  points  covering  the 
entire  urban  region.  A  more  useful  modeling  tool  for  studying  the  impact  of 
airport  emissions  would  be  a  trajectory-based  photochemical  model,  such  as  the 
new  ELSTAR  model  (Lloyd  et  al.,  1979),  which  calculates  concentrations  along  a 
specific  path  or  trajectory  of  an  air  parcel.  A  trajectory-based  photochemic¬ 
al  model  could  economically  study  the  effects  of  various  airport  emissions 
control  strategies  with  a  detailed  cons idera t  ior  of  both  NO2  and  ozone 
formation  in  the  vicinity  of  the  airport. 

2.3  METEOROLOGICAL  ASPECTS  OF  AIRPORT  AIR  POILUTION  WORST  CASE  ANALYSIS 

The  ambient  levels  of  air  pollutants  depend  not  only  on  the  amount  of 
pollutant  emitted  into  the  atmosphere  buc  also  upon  the  prevailing  meteorolo¬ 
gical  conditions.  The  dispersive  capability  of  the  atmosphere  depends  upon 
such  meteorological  parameters  as  the  wind  speed  and  the  vertical  temperature 
profile.  Of  course,  the  wind  direction  also  plays  an  important  role  when 
considering  any  particular  source-receptor  pair.  These  parameters  vary 
hourly,  diurnal ly,  and  seasonally  as  both  small-  and  large-scale  weather 
patterns  change. 

The  air  qualitv  effects  of  the  prevailing  meteorological  conditions  are 
not  the  same  for  all  sources.  Elevated  sources  have  their  greatest  impact 
during  unstable  or  neutral  atmospheric  conditions.  High  wind  speeds,  which 
may  occur  during  periods  of  neutral  stability  will  also  reduce  plume  rise 
and  bring  plumes  to  the  ground  closer  to  the  source  than  under  lighter  wind 
cases,  thus  diminishing  the  effect  of  greater  initial  dilution.  Under  stable 
atmospheric  conditions  or  during  a  temperature  inversion  (ambient  temperature 
increasing  with  height)  the  plume  from  an  elevated  point  source  may  remain 
aloft  and  intact  for  many  kilometers  downwind.  If  the  inversion  layer  exists 
above  the  elevation  of  the  source,  wh i 1 e  the  layer  below  is  unstable,  maximal 
concentrations  may  occur  at  short  distances  from  the  base  of  the  source. 


36 


Conversely,  the  contribution  of  ground-level  sources  to  ambient  air 
pollutant  levels  is  usually  maximized  under  stable  or  inversion  conditio  is 
provided  wind  direction  is  steady.  At  airports,  it  is  the  ground  level 
emissions  which  predominantly  impact  nearby  receptors.  When  stable  conditions 
occur  at  the  surface,  the  thermally  turbulent  surface  mixing  layer  has 
effectively  zero  thickness  and  pollutants  emitted  at  low  levels  are  trapped 
near  the  surface.  The  dilution  of  pollutants  by  mechanically  induced  atmos¬ 
pheric  turbulence  is  also  inhibited.  In  these  cases  the  concept  of  mixing 
depth  is  not  useful  since  the  vertical  dimension  of  the  plume  usually  remains 
smaller  than  the  thickness  of  the  stable  layer.  It  should  be  noted,  however, 
that  light  and  directionally  variable  winds  are  also  often  associated  with 
such  stable  atmospheric  conditions.  These  wind  conditions  serve  to  spread  the 
impact  zone  for  any  aggregate  of  low  level  sources.  This  spreading  causes 
contributions  from  many  spatially  distributed  sources  to  partially  impact  any 
particular  receptor  area  over  an  hour's  period,  but  also  mitigates  the  impact 
of  any  single  source.  Thus,  the  wind  direction  meander  acts  like  a  spatial 
lowpass  f  Iter:  reducing  the  highest  observed  concentration  by  distribu¬ 
ting  "hot  spot"  impacts  to  other  nearby  cross  wind  receptors. 

A  study  of  the  meteorological  factors  involved  in  pollutant  levels 
monitored  around  Los  Angeles  International  Airport  (Hallanger,  1974)  showed 
that  the  potential  for  airport  related  air  pollution  at  any  given  moment  was 
greatest  during  stable  atmospheric  conditions  accompanied  by  light  and  direc¬ 
tionally  variable  winds. 

To  assess  the  frequency  of  high  pollution  potential  a  statistical 
compilation  of  the  weather  conditions  observed  at  each  of  the  five  study 
airports  was  examined.  Table  2.2  lists,  by  airport,  the  time  periods  of 
meteorological  observations  examined  and  the  frequency  of  stable,  low-moderate 
wind  speed  conditions  having  high  pollution  potential.  Table  2.2  shows  that 
low  wind  speed  stable  cases  are  indeed  frequent  at  all  of  these  airports; 
however,  it  should  be  mentioned  that  a  large  percentage  of  these  cases  are 
associated  with  low  airport  activity,  nighttime  conditions. 

In  order  to  assess  the  "worst-case"  frequency  for  maximum  concentra¬ 
tions  at  individual  airports  the  orientation  of  "critical  receptors"  with 
respect  to  areas  of  high  emissions  density  must  be  considered.  A  critical 
receptor  is  defined  as  an  area  which  is  accessible  to  the  public  for  periods 
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Table  2.2.  Frequency  of  Poor  Atmospheric  Dispersion 
Conditions3  at  Five  Major  Airports 


Airport 

Time  Period 
Examined 

Frequency*5 

DCA 

1968-1972 

30.7 

LAX 

1955-1964 

35.1 

LGA  (NY) 

1965-1970 

19,7 

ORD 

1960-1964 

29.0 

Dulles 

1966-1970 

37.0 

aDefined  here  as  stable  (E  of  F)  s t ra t i f ic at  ion  with 
a  wind  of  1-5  m/sec  from  any  direction. 
bF 

requencies  are  based  on  a  five  year  annual  average 
period . 

exceeding  an  hour,  and  is  relatively  near  an  area  of  high  emission  density. 
This  definition,  in  itself,  requires  some  knowledge  of  the  pollutant  concen¬ 
tration  patterns  associated  with  each  wind  direction  under  a  range  of  wind 
speeds  and  atmospheric  stabilities.  Critical  receptors  are  most  precisely 
defined  by  evaluating  a  series  of  dispersion  model  analyses  which  cover 
the  range  of  potentially  critical  cases.  Temporal  variation  of  emission 
patterns  must,  of  course,  also  be  considered.  Preliminary  estimates  may  be 
made,  however,  based  upon  the  knowledge  of  the  receptor  map,  the  emissions 
map,  and  wind  frequency  tables.  The  worst  case  frequencies  identified  here 
are  based  upon  those  considerations  and  model  analyses  carried  out  for  LAX, 
JFK,  and  ORD  airports  as  reported  in  Volme  II  of  this  report.  The  stability 
classification  scheme  used  here  is  the  well-known  method  of  Turner,  1964. 

While  this  method  may  be  less  precise  than  one  which  uses  actual  onsite 
measurement  of  turbulence  intensity,  it  is  generally  the  only  method  available 
for  prospective  studies  of  impacts  based  on  historical  meteorological  records. 

For  LAX  airport  the^e  are  two  potentially  critical  receptors:  (1) 
the  terminal  area  and  (2)  the  restaurant  and  golf  course  to  the  East  of 
runway  24L.  The  worst  case  wind  directions  for  the  terminal  as  a  receptor  are 
N  and  E  to  ESE.  For  the  restaurant  and  goif  course  receptors  W  to  WNW  winds 
are  most  important.  North  winds  result  in  higher  concentrations  but  are 
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less  frequent  than  E  or  ESE  winds.  For  LAX  a  previous  study  made  available 
separate  frequency  distribution  data  for  the  6  to  9  AM  hour  interval,  which  is 
the  period  of  maximum  aircraft  takeoff  activity  coinciding  with  stable  condi¬ 
tions  . 

Table  2.3  shows  LAX  wind  frequencies  for  each  of  these  critical  receptor 
directions  for  the  worst  case  1-3  knot  and  4-6  knot  wind  speed  classes  during 
stable  (E  and  F)  conditions.  (For  1-3  knots,  only  F  stability  exists;  for  4-6 
knots,  60  to  70%  of  the  cases  are  F  stability;  for  7-10  knots  only  E  stability 
exists).  Although  many  air  quality  studies  indicate  that  the  highest  concen¬ 
trations  arising  from  a  combination  of  ground  level  sources  are  associa¬ 
ted  with  1-3  knot  winds  and  stable  (F)  turbulence  levels:  consideration  of  the 
buoyant  plume  rise  of  hot  aircraft  exhaust  plumes  leads  to  the  conclusion 
that,  for  aircraft,  wind  speeds  of  4-6  ( o.  even  7-10)  knots  are  necessary  to 
bring  the  plume  to  the  ground  at  receptors  close  to  these  sources  (see  Section 
4.5).  Frequencies  of  these  stable  cases  for  the  6  to  9  AM  period  are  also 
given  in  Table  2.3. 

The  important  wind  directions  for  DCA  airport  are  defined  by  the  air 
quality  monitoring  experiments  of  1972  and  1979.  The  N  to  ENE  sectors  are  of 
principal  interest.  Critical  airport  receptors  include  the  southern  terminals 
and  the  parking  areas  SW  to  WSW  of  the  ramp  areas.  SSE  winds  would  also  be 
important  for  aircraft  emissions  impacts  at  the  public  terminal  areas.  For 
that  reason  Table  2.4  lists  frequencies  of  low  wind  speed  stable  conditions 
for  each  of  these  directions.  The  N  and  SSE  directions  each  account  for 
more  than  100  hours  per  year  of  stable  conditions  with  less  than  a  6  knot 
wind  speed.  However,  in  determining  the  frequency  of  worst-case  occurrences 
it  is  important  to  consider  the  number  of  hours  that  these  meteorological 
conditions  coincide  with  periods  of  high  a i r por t /a i rc r a f t  activity.  At 
PCA  there  is  virtually  no  activity  between  midnight  and  7:00  AM  (per  observa¬ 
tions,  Feb.  19  79  arid  OAG)  .  Since  stable  conditions  occur  predominantly  at 
night,  there  are  likely  to  be  fewer  than  50  hours  per  year  of  "airport  worst 
case"  conditions  at  UCA . 

For  JFK  a  280  degree  direction  has  been  assumed  for  the  modeling  in 
Volume  II,  where  impacts  of  2  knot  stable  winds  upon  the  terminal  building  were 
evaluated.  The  important  directions  for  off-airport  receptors  are  more 
likely  SE  to  S  winds.  Table  2.5  compares  the  frequencies  for  low  speed 
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Table  2.3.  Annual  Percentage  Frequencies3  of 
Stable  Stratification  at  LAX 


All  Hours  Hours  06-09 


Wind  Direction 

Wind 

1-3 

Speed 

4-6 

(Knots) 

7-10 

Wind 

1-3 

Speed 

4-6 

(Knots) 

7-10 

N 

0.54 

0.67 

0.34 

0.30 

0.47 

0.25 

E 

1  .61 

2.36 

0.26 

1.10 

2 .75 

0.23 

ESE 

0.95 

1  .31 

0.13 

0.85 

1  .10 

0 .10 

W 

1  .03 

2.30 

1.75 

0 .17 

0 .27 

0.06 

WNW 

0  .48 

0.81 

0.28 

0.19 

0.16 

0.06 

Table  2.4.  Annual  Percentage  Frequencies3  of  Stable 
Stratification  at  DCA.  All  Hours 


Wind 

Direction 

Wind 

1-3 

Speed 

4-6 

(Knots) 

7-10 

N 

0.34 

0 .96 

0.41 

NNE 

0.15 

0.50 

0.30 

NE 

0.24 

0.68 

0.31 

ENE 

0 . 10 

0.45 

0.12 

SSE 

0.50 

1 .30 

0 .12 

aNote  that  0.02%  is  approximately  two  hours  per  year. 
Any  calculated  concentration  exceeding  a  federal 
ambient  air  quality  standard  and  associated  with  a 
greater  frequency  could  result  in  a  violation. 
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Table  2.5.  Annual  Percentage  Frequencies3  of  Stable 
Stratification  at  JFK.  All  Hours 


Wind  Speed  (Knots) 
Wind  Direction  1-3  4-6  7-10 


ESE 

0.07 

0.13 

0.05 

SE 

0 .06 

0.25 

0.09 

SSE 

0.06 

0.25 

0.09 

S 

0.02 

0.11 

0.08 

W 

0.11 

1  .44 

0.76 

stable  flow  in  these  three  sectors  with  the  much  higher  frequency  occuring  for 
W  winds . 

For  ORD  the  worst  case  directions  considered  in  Volume  II  are  NNE  and 
SE.  These  directions  were  also  selected  to  evaluate  impacts  upon  the  terminal 
area.  However,  the  communities  to  the  North  of  the  airport  may  also  be 
defined  as  important  off-airport  receptors.  Thus  the  high  frequency  of  low 
speed  stable  winds  from  the  S  are  also  shown  in  Table  2.6  for  comparison  with 
the  NE  and  SE  frequencies. 

Further  consideration  of  the  results  of  the  DCA  and  Dulles  model 
validation  experiments  has  shown  that  any  final  determination  of  the  critical 
wird  speed  for  a  worst  case  condition  at  a  particular  airport  is  highly 
sensitive  to  the  specific  geometric  relationship  of  critical  receptors  to 
areas  of  high  emission  density.  It  is  therefore  expected  that  some  model 
sensitivity  analyses  will  be  necessary  before  conclusions  about  a  single 
critical  wind  speed  and  stability  can  be  extended  to  any  other  airport. 
However,  on  the  basis  of  the  frequency  tables  presented  here,  it  can  be 
concluded  that  stable  conditions  accompanied  by  a  wide  range  of  wind  speeds 
and  directions  occur  at  all  of  the  listed  airports. 

On  the  basis  of  the  mon'toring  experience  at  Dulles  International  and 
Washington  National  Airports,  the  definition  of  worst  case  conditions  for 
on-airport  monitoring  stations  located  within  500  m  of  taxiways  and  runways 
was  a  somewhat  different  problem  from  that  outlined  above.  For  detection  of 
impact  of  CO  pollutant  emissions  from  individual  taxiing  aircraft  at  Dulles, 
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Table  2.6.  Annual  Percentage  Frequencies3  of  Stable 
Stratification  at  ORD .  All  Hours 


Wind  Speed  (Knots) 
Wind  Direction  1-3  4-6  7-10 


NE 

0. 

14 

0 

.68 

0 

.25 

SE 

0  . 

32 

0 

.75 

0 , 

.08 

SSE 

0  . 

35 

0 

.82 

0, 

.16 

S 

0. 

79 

2 

.28 

0, 

.62 

SSW 

0  . 

60 

1 

.25 

0 

.55 

sw 

0. 

55 

1 

.30 

0 

.00 

aNote  that  0.02%  is  approximately  two  hours  per  year. 
Any  calculated  concentration  exceeding  a  federal 
ambient  air  quality  standard  and  associated  with  a 
greater  frequency  could  result  in  a  violation. 


the  critical  wind  speed  was  about  10  mph  (Smith,  et  al . ,  1977).  This  speed 
was  sufficient  to  delay  plume  rise  until  after  the  exhaust  plume  was  advected 
past  the  monitoring  sites.  In  addition,  the  higher  wind  speeds  reduced  total 
final  plume  rise  and  therefore  increased  relative  ground  level  concentrations. 
However,  at  DCA  the  highest  N0X  and  CO  total  concentrations  were  measured 
during  episodes  of  light  and  directionally  variable  winds  accompanied  by 
moderately  stable  temperature  stratification,  though  it  should  be  noted  that 
high  concentrations  were  simultaneously  observed  throughout  the  D.C.  area, 
thus  obscuring  the  role  of  aircraft  emissions. 

thus,  one  arrives  at  a  somewhat  complex  and  uncertain  characterization 
of  the  meteorological  conditions  leading  to  worst  case,  ground-level  concentra¬ 
tions  at  airports.  Considerations  include: 

a)  the  predominant  ground-level  nature  of  the  sources 
suggests  that  F  stability  and  very  low  wind  speeds  (e.g., 

2  mph)  leads  to  the  highest  concentrations; 

b)  the  plume  rise  dynamics  associated  with  the  buoyancy  of 
the  aircraft  exhaust  suggests  that  peak  ground-level 
concentrations  will  occur  at  higher  wind  speeds  (e.g., 

7-10  mph ) ; 

c)  the  strong  dynamic  peaking  of  aircraft  operations  at  most 
airports  biases  worst  case  overall  conditions  toward  less 
stable  atmospheric  conditions; 
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d)  the  proximity  of  large  bodies  of  water  or  urban  areas 
adjacent  to  many  airports  influences  the  range  of  possible 
stabilities  and  is  generally  not  considered  in  worst  case 
air  duality  assessments;  and 

e)  the  large,  engine  generated  turbulence  causes  engine 
emitted  pollutants  to  undergo  an  initial  mixing  that  is 
somewhat  stability  class  independent. 

These  considerations  together  with  modeling  approximations  and  limita¬ 
tions  discussed  in  Volume  II  suggest  tha,'.  the  compromise  choice  of  E  stability 
and  a  wind  speed  of  ^  2  mph  might  best  characterize  worst  case  conditions. 
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3  THE  WASHINGTON  NATIONAL  AIRPORT  (DCA)  STUDY 


3 . 1  INTRODUCTION 

A  major  constituent  of  the  United  States  program  to  determine  the 
impact  of  aircraft  emissions  on  airport  air  quality  consists  of  the  recently 
completed  measurement  of  pollutants  at  Washington  National  Airport  (DCA).  The 
air  quality  monitoring  program  at  DCA,  sponsored  by  the  EPA  and  FAA,  was 
designed  to  measure  NO,  N0X,  and  CO  within  0.3  km  of  queueing  and  takeoff 
operations  to  permit  an  unambiguous  assessment  of  the  aircraft  contribution  to 
nearby  pollutant  levels. 

In  this  section  the  DCA  monitoring  program  is  briefly  described, 
followed  by  statistical  analyses  of  the  CO  and  oxides  of  nitrogen  (i.e.,  NO, 

NO2 ,  and  N0X  =  NO  +  NO2)  concentration  data.  Two  methods,  background  subtrac¬ 
tion  and  pulse  integration,  are  described  and  applied  to  the  data  to  estimate 
the  aircraft  portion  of  observed  pollutant  levels.  These  hourly  data,  together 
with  extrapolations  of  the  cumulative  frequency  distributions  of  these  data, 
are  then  compared  to  relevant  National  Ambient  Air  Quality  Standards  (NAAQS) . 

In  addition,  regression  analysis  is  employed  to  estimate  the  impact  of  air¬ 
craft  operations  on  annual  average  NO2  levels  adjacent  to  the  runway. 

Further  analyses  of  the  N0X  data  include  consideration  of  the  N0/N0x 
ratio  as  a  function  of  plume  travel  time  and  a  comparison  of  the  predictions 
of  four  dispersion  models  with  measured  N0X  levels.  In  addition,  inidividual 
N0X  takeoff  plumes  passing  across  the  array  of  receptors  are  analyzed  to 
extract  initial  takeoff-mode  plume  dimensions  for  subsequent  input  to  air 
quality  dispersion  models.  Finally  the  results  of  an  AVAP  model  simulation, 
using  these  latest  plume  dimensions,  is  presented  to  indicate  the  quantitative 
agreement  between  model  predictions  and  monitoring  data  in  suggesting  the 
potentially  serious  impact  of  aircraft  relative  to  a  possible  NO2  hourly 
standard  in  the  0.2-0. 5  ppm  range. 

3.2  THE  DCA  MONITORING  PROGRAM 

Aerometric  data  was  recorded  from  rapid  response  pollutant  monitors 
(TECO  and  Bendix  chemiluminescent  monitors  for  N0/N0x  and  Bendix  NDIR  and 
Ecolyzer  electrochemical  monitors  for  CO)  at  the  six  sites  shown  in  Fig.  3.1 


NORTH 
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Fig.  3.1.  Monitoring  Site  Locations  at  DCA. 

Sites  were  chosen  to  focus  on  the 
pollution  clouds  from  the  takeoff 
and  queueing  modes. 


during  the  period  January  15-Feburary  27,  1979.  With  the  exception  of  station 
2,  which  was  sited  specifically  to  obtain  background  pollutant  levels  under 
northerly  to  easterly  wind  conditions,  Table  3.1  indicates  the  air  quality 
parameters  measured  at  each  station.  In  addition,  airport  meteorological  data 
were  supplemented  by  measurement  of  wind  speed  and  direction,  wind  azimuth  and 
elevation  angles,  temperature,  vertical  temperature  gradient,  and  dew  point 
temperature  at  a  site  r  20  m  north  of  station  5.  During  the  latter  half  of 
the  program  a  decibel  meter  nearby  station  4  measured  aircraft  engine  noise 
and  provided  a  convenient  time  reference  for  takeoff  and  landing  operations. 

Pollutant  concentrations,  noise  level,  and  meteorological  parameters 
were  recorded  on  three  independent  systems:  individual  strip  chart  recorders 

for  each  instrument,  a  set  of  multi-pen  strip  chart  recorders  synchronized  by 
an  external  time  reference,  and  on  magnetic  tape  via  a  13-channel  data  acqui¬ 
sition  system  (DAS) .  While  hourly  average  concentrations,  as  extracted  from 
the  strip  charts,  have  been  previously  analyzed  and  compared  to  PAL  predicrlons 
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Table  3.1a.  Air  Quality  Parameters  Measured  at  Each  Station 


Parameter 

Station 

1 

3 

4 

5 

6 

N0X 

X 

X 

X 

X 

- 

NO 

X 

X 

- 

- 

- 

CO 

X 

X 

X 

- 

X 

CO  Ecolyzer 

- 

•• 

X 

- 

X 

X:  measured 
- :  not  measured 


(EPA,  1980).  thiu  report  utilizes  the  data  provided  via  the  DAS;  thus  a 
brief  description  of  its  operational  characteristics  and  subsequent  data 
analysis  assumptions  are  essential. 

The  DAS,  operating  basically  as  a  recording  voltmeter,  interrogates 
each  of  its  15  inputs  once  every  3  sec,  converts  these  analog  voltage  signals 
to  digital  form  (0-255),  and  then  writes  these  numbers  onto  magnetic  tape. 
Assumptions  concerning  the  conversion  from  engineering  units  (e.g.,  ppm)  to 
voltage  and  subsequently  back  to  engineering  units  are  based  or.  manufacturer 
specifications  and  EPA  calibrations.  In  addition,  an  independent  EPA  audit 
(Arey  et  al ,  1979)  confirmed  that  all  of  the  CO  and  N0/N0x  instruments 
exhibited  satisfactory  or  excellent  response  to  the  audit  concentrations. 
Hence,  the  air  quality  data  in  this  report  should  be  considered  accurate  to 
the  o<10 %  level. 

Some  of  the  program  objectives  can  be  visualized  with  the  aid  of  Figs. 
3.2a  and  3.2b  which  show  sample  sets  of  "strip  charts"  produced  using  the  DAS 
acquired  magnetic  tapes.  Numerous  comparisons  between  these  synthesized  strip 
charts  of  pollutant  and  noise  level  and  actual  strip  charts  of  pollutant 
levels  and  log  sheets  of  aircraft  operations  permitted  confident  matching  of 
NO  and  NCX  pollutant  pulses  with  specific  aircraft  operations.  Figure  33 
illustrates  the  analysis  possibilities  created  through  this  high  sampling  rate 
approach.  Takeoff  plumes,  easily  identified  by  the  N0X  rise  above  background 
levels  and  the  presence  of  a  noise  pulse  accompanying  takeoff,  are  seen  to 
disperse  (i.e.  reduce  in  peak  concentration  while  spreading  out  in  time)  as 
they  move  from  the  nearby  station  1  to  the  more  distant  station  3.  The  acutal 
impact  of  these  plumes  may  be  estimated  via  a  background  subtraction  (i.e.  as 
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Table  3.1b.  List  of  Monitoring  Equipment 


Equipment 

Shelter 

Parameter 

Instrument 

Voltage 

Site  No. 

Shelter  No. 

Dimens  ions 

Measured 

Manufacturer 

Output 

I 

Self  Propelled 

27 1 x8 1 x  14'* 

NO,  N0X 

Thermo  Electron 

10V 

EPA  #313 

Company  (TECO) 

Bendix 

10V 

IV 

Single  pen 

strip  chart  recorders 

(SCR) 

for  each  parameter. 

2 

Self  Propelled 

27'x8’x30'** 

CO 

Bendix 

10V 

Background 

EPA  #376 

03 

Das ib i 

HC 

Bendix 

IV 

N0X 

Bendix 

IV 

Wind  Direc- 
t  ion  & 

Ve loc i ty 

Climatronics 

Single  pen 

strip  chart  recorders 

(SCR) 

for  each  parameter  which  is  also 
into  data  processing  computer. 

input 

3 

Trai ler 

8 ' xl4 1 xl4 ' * 

CO 

Bendix 

10V 

EPA  #577 

IV 

N0X  NO 

Wind  Direc- 

TECO 

10V 

tion  & 
Velocity 

Climel 

Wind  Direc- 
t  ion  & 

Ve loci ty 
(2  Dimen¬ 
sions  ) 

MRI  Vector  Vane 

Temperature 

and 

Temperature 

Gradient 

Cl ime  t 

Single  pen 

SCR  for  each  paramter. 

4 

Self  Propelled 

Same  as 

N0X 

Bendix 

IV 

EPA  #315 

Site  1 

CO 

Bendi x 

10V 

IV 

Single  pen 

SCR  for  each  parameter 

4  multi  pe n 

SCR  coordinated  to  common 

time  reference  to  simultaneously 

record 

concentrations  at  Sites  1,  3,  4, 

5,  6. 

Data  logger 

computer  to  record  15 

chan- 

nels  of  data  from  Sites  1,  3,  4, 

5,  6. 

5 

Trailer 

Sane  as 

N0X 

Bend  ix 

IV 

EPA  #5  75 

Site  3 

Single  pen 

SCR  for  each  parameter 

6 

Trailer 

Sane  as 

CO 

Bendix 

10V 

EPA  #5  76 

Site  3 

CO 

Energet ic 

Sc iences  Co .  ( 2 ) , 

mobile 

HC 

Beckman  400 

Single  pen 

SCR  for  each  parameter 

•Includes  Air  Intake  Probe. 


••Includes  22  foot  high  wind  set. 


NOX  -  TRAILER  1 


"WlC 


"Tr»  in; 

Tint  iMMnrii  -  irrCBRuMi  1979 

WD  -  TRRILER  5 _ 


Tm:' 


TrTo  T7i  Tnc  Vmo 

unc  rnmni  -  is  rcBRu«Ri  1979 
HS  -  TRAILER  5 


TtTo  inc  irx 

Ti»c  ttmii  *  ts  rtam^ri  ir9 


Fig.  3.2a.  Synthesized  Strip  Chart  from  DAS  Tapes  for  a  Typical  2-hr 
Period  of  Northeast  Winds  Showing  Noise  and  N0X  Signals 


OEGREES 

190.0  360. 
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seen  by  the  dotted  line)  or  through  estimation  of  the  area  under  the  major, 
aircraft  related  pulses  Ci.e.  as  indicated  by  the  hatched  area)  via  a  tech¬ 
nique  requiring  a  threshold  concentration  above  background  to  be  reached 
before  integration  and  thus  inclusion  of  a  particular  pulse.*  Analysis  of 
these  individual,  aircraft  generated  N0/N0x  signals,  with  the  aid  of  an 
integrated  Gaussian -puf f  type  model  that  simulates  the  high  thrust,  high 
NOx~eraitting  takeoff  roll  and  parameterizes  the  N0X  plume  behavior  in  terras  of 
initial  plume  dimensions  and  subsequent  plume  transport  and  dilution  rates  is 
reported  in  Section  3.6.  These  new  plume  parameterizations  permit  more 
realistic  prediction  of  peak,  short-term,  N0X  levels  near  runways.  Analysis 
of  the  CO  signals,  which  are  more  difficult  to  associate  with  single  aircraft 
operations  due  to  the  complexity  of  aircraft  queueing,  will  now  be  discussed. 

3.3  HOURLY  AVERAGE  MONITORING  DATA  FOR  CO 

While  a  number  of  agencies  and  groups  concerned  with  the  air  quality 
impact  of  major  airports  have  undertaken  monitoring  programs  as  well  as 
theoretical  studies  based  on  the  use  of  atmospheric  dispersion  algorithms,  a 
recent  review  of  these  efforts  by  Lorang  (1978)  suggests  that  the  issue  is 
particularly  confusing  with  respect  to  carbon  monoxide.  Ambient  measurements 
conducted  at  Los  Angeles  International  by  Thayer  et .  al.  (1974)  and  Washington 
National  Airports  by  Platt  et.  al.  (1971)  were  ambiguous  as  to  their  attribu¬ 
tion  of  measured  ievels  to  either  aircraft  or  non-aircraft  sources.  Similarly, 
initial  modeling  predictions  using  the  NREC  model  [Platt  et  .  al  .  (1971)] 
indicated  the  likelihood  of  violations  of  both  the  1-hr  (40  -ig/nr  .  35  ppm)  and 
8-hr  (10  ]ig/m^=  9  ppm)  standards  for  CO,  while  a  more  recent  modeling 
exercise  [Yamartino  end  Rote,  (1978)]  for  LAX  suggests  "worst  case"  hourly  CO 
concentrations,  attributable  to  aircraft  alone,  of  less  than  5  ppm  beyond  1000 
ft  from  the  aircraft  queueing  area  and  2  ppm  or  less  at  the  passenger  terminals. 
Given  the  uncertainties  generated  by  these  monitoring  and  modeling  experiences, 
the  EPA  and  FAA  chose  to  monitor  CO  near  an  aircraft  queueing  area  at  DCA . 

*These  techniques  will  subsequently  be  referred  to  as  "background  subtracted" 
and  "pulse  integrated"  concentrations. 
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In  order  to  facilitate  comparison  with  air  quality  standards,  hourly 
average  concentrations  were  computed  from  the  high  repetition  rate  data 
(i.e.,  1200,  3-sec  scans  *  1  hr).  Figure  3.4  shows  synthesized  "strip  charts" 
for  hourly  averaged  CO  concentrations  for  the  entire  monitoring  period.  In 
addition,  a  statistical  summary  of  the  hourly  average  CO  concentrations  is 
given  in  Table  3.2  for  the  708  hours  of  nearly  continuous  DAS  operation. 
Reductions  in  sample  size  result  from  a  combination  of  equipment  calibrations 
and  failures  and/or  deliberate  changes  in  the  allocation  of  the  15  DAS  channels 
among  the  20  channels  of  incoming  data.  DAS  logging  of  station  3  CO  was 
suspended  after  10  days  by  just  such  a  DAS  channel  pt ior it ization  and  because 
station  3  levels  were  not  markedly  different  in  magnitude  or  time  dependence 
from  those  at  station  4. 

Figure  3.5  shows  a  histogram  of  hourly  CO  concentrations  observed  at 
station  6.  This  nearly  lognormally-shaped  distribution  peaks  near  1  ppm, 
which  is  not  surprising  for  a  receptor  surrounded  by  a  large  metropolitan 
area.  The  more  interesting  higher  concentration  values  and  interstation 
comparisons  are  more  conveniently  presented  in  the  cumulative  frequency 
distributions  of  Fig.  3.6.  One  observes  that  concentrations  at  monitors  1  and 
6,  located  within  500  ft  of  the  aircraft,  are  less  than  3  ppm  90%  of  the  time. 
Concentrations  at  station  4,  located  1000  ft  from  the  queueing  aircraft  and 
at  a  distance  comparable  to  that  where  the  closest  public  exposure  might 
occur,  are  not  observed  is  exceed  3.2  ppm. 

It  should  be  pointed  out  that  all  station  1  and  6  concentrations  in 
excess  of  10  ppm  correspond  to  a  single  episode  that  persisted  from  2100  hrs 
February  22  through  0500  hrs  February  23,  under  near  calm,  stable  conditions. 

CO  levels,  wind  speed  and  direction,  and  noise  level  (indicting  departures) 
for  this  period  are  shown  in  Figure  3.7.  Though  the  1-hr  standard  was  never 
approached,  the  8-hr  standard  was  exceeded  by  38%.  More  careful  analysis  of 
this  rather  unique  period  suggests  that  airport  snowplows,  operating  intensive¬ 
ly  in  the  vicinity  of  station  1  and  6,  may  have  been  significant  contributors 
to  these  maximum  observed  CO  values.  In  addition,  the  generally  westerly  nature 
of  the  wind  would  not  have  been  conducive  to  transporting  pollutants  from  thp 
queueing  area.  However,  because  of  the  very  low  wind  speeds  2  mph  and  the 
highly  variable  wind  direction  (  -  60°),  significant  contribution  of  aircraft 

cannot  be  precluded.  Further,  linear  extrapolation  (appropriate  here  assuming 
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Synthesized  Strip  Charts  of  Hourly  Average  CO  Concentrations 
Wing  Speed,  and  Wind  Direction  for  the  Period  Jan.  24-Fob.  2 
1979.  These  averaged  data  are  computed  from  3-sec  interval 
DAS  data. 


Cumulative  Frequency  Distributions  of  Hourly  CO  Concentrations  at 
Stations  1,  4,  and  6. 
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Fig.  3.7.  Synthesized  Strip  Chart  from  DAS  Tapes  for  the  8  hour 
Period  Beginning  2100  February  22  and  Including  the 
Entire  CO  Episode 
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lognormal  concentration  distributions)  of  the  cumulative  frequency  distribution 
upward  from  below  10  ppm  suggests  that  these  greater-than-1 0  ppm  values  could 
occur.  However,  there  is  some  question  whether  the  observation  of  high 
concentrations  at  locations  in  such  close  proximity  to  the  aircraft  is  relevant 
to  the  question  of  NAAQS  violations.  Lastly,  we  note  that  similarly  high  CO 
values  were  observed  during  this  period  by  the  local  pollution  control  agency 
monitors  throughout  the  entire  Washington  area. 

Extrapolation  of  the  curves  in  Fig.  3.6  out  to  the  99.99 %  probability 
level  (i.e.,  1  hr  in  10,000  or  approximately  once  per  year),  while  informative, 
should  be  viewed  with  caution,  as  not  only  is  such  extrapolation  based  on  only 
1/12  of  a  year's  data,  but  all  these  data  come  from  a  single  contiguous  set  of 
hours  rather  than  from  a  random  selection  of  hours  throughout  the  year. 

Figures  3.8-3.10  represent  an  attempt  at  setting  bounds  on  the  aircraft 
contribution  to  the  three  curves  in  Fig.  3.6.  In  each  of  these  figures  the 
uppermost  curve  represents  the  distribution  of  total  hourly  average  CO  concen¬ 
trations  (same  as  Fig.  3.6),  and  thus  represents  the  maximum  possible  impact 
of  aircraft.  The  next  lower  lying  curves  (labeled  2)  represent  the  "background' 
subtracted  concentrations,  where,  in  lieu  of  station  2  observations,  'pack- 
ground"  is  defined  as  the  average  of  the  12  minimum  concentrations  observed 
during  the  consecutive  5-min  periods  making  up  the  hour.  Coincidence  of 
curves  1  and  2  at  the  lower  concentrations  arises  when  actual  background  is 
below  instrument  threshold  and  thus  yields  a  "zero  background  upon  subtraction 
The  lowest  lying  curves  (labeled  3)  represent  the  average  concentrations  con¬ 
tributed  by  pollution  pulses  rising  at  least  0.35  ppm*  above  a  15-min  average 
"background"  and  subsequently  corrected  upward  by  the  factor  1/erf  [)/ 1  n  ( Cp/  Cy ) 
where  Cp  is  the  peak  pulse  concentration  above  background,  to  compensate  for 
this  0.35  ppm  "barrier"  CT.  Curve  3  thus  isolates  the  contribution  of  nearby 
transient  pollution  sources.  The  fixed  size  of  the  "barrier"  accounts  for 
convergence  of  curves  2  and  3  at  high  concentrations.  Thus,  the  actual 
aircraft  (or  more  properly,  local  source)  contribution  to  observed  concentra¬ 
tions  probably  lies  somewhere  within  the  band  defined  by  curves  1  and  3. 


*Th is  threshold  barrier  and  the  0.035  ppm  barrier  for  the  N0X  analysis  were 
chosen  by  searching  for  the  plateau  region  which  is  observed  when  the 
pulse  integrated  concentration  is  plotted  as  a  function  of  threshold 

barrier  level. 
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Keeping  in  mind  previous  cautions  about  extrapolation  of  these  curves, 
one  notes  that  at  station  4,  located  =  1000  ft  from  the  queueing  area,  a 
maximum  hourly  CO  concentration  of  -  5  ppm  may  be  expected  about  once  per 
year.  This  result  is  consistent  with  "worst  case"  predictions  for  LAX, 
seen  in  Fig.  3.11,  where  a  fleet  mix  and  queueing  emissions  comparable  to  DCA 
are  assumed.  Further  extrapolations  or  generalizations  from  the  DCA  results 
to  other  airports  should  be  tempered  by  the  following  considerations: 

•  DCA  is  closed  between  10:00  p.m.  and  6:00  a.m. 

These  nighttime  hours  are  associated  with  stable 
atmospheric  conditions  and  thus  potentially  poor 
pollutant  dispersion  conditions. 

•  Runway  36  at  DCA  is  shared  between  arrivals  and 
departures.  Though  not  an  unusual  situation, 
airports  having  dedicated  departure  runways  should, 
for  the  same  departure  rate,  have  shorter  queueing 
times  and  correspondingly  reduced  CO  emissions  and 
concentrations . 

•  Nearly  all  operations  at  DCA  are  by  medium  range 
jets  (e.g.,  727,  737,  DC9)  primarily  using  the  JT8D- 
17  engine.  This  engine  has  a  relatively  low  CO 
emission  rate  at  idle  (~40  lbs/hr)®  compared  to  some 
other  engines  (e.g.,  88  lbs/hr  for  the  CF6-5oC  and 

=  140  lbs/hr  for  the  JT9D-7  and  RB-2 1 1- 22B)  .8 


3.4  HOURLY  AVERAGE  MONITORING  DATA  FOR  OXIDES  OF  NIRTOGEN 

The  issue  of  oxides  of  nitrogen  (N0X)  impacts  created  by  aircraft  is, 
as  with  CO,  a  localized  "hot-spot"  problem  related  to  existing  and  potential 
NAAQS .  Unfortunately  the  issue  is  further  complicated  by  the  fact  that 

•  present  and  possible  future  NAAQS  standards  pertain  to 
NO2  levels  and  not  NOx  levels.  (N0X  «  NO  +  NO2) 

•  there  is  presently  only  an  annual  average  NAAQS  of  0.05 
ppm  NO2  through  a  one  hour  average  standard  in  the  range 
0.2-0. 5  ppm  is  presently  being  considered  by  the  EPA 

•  plume  dispersion,  while  reducing  the  concentration  of 
inert  species,  will  entrain  more  amibient  oxidant  result¬ 
ing  in  further  conversion  of  engine  emitted  NO  to  NO2 ; 
thus,  NO2  levels  will  peak  at  some  distance  downwind  of 
the  aircraft 

•  the  peak  NO2  attributable  to  aircraft  is  a  function  of 
existing  ambient  levels  of  NO,  NO2,  O3,  and  sunlight. 
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The  potential  for  violation  of  possible  NO2  peak  hourly  standards  at 
airports  has  recently  been  reveiwed  by  Jordan  and  Broderick  (1978).  Finding 
that  both  worst  case  modeling  predictions  and  previous  monitoring  results  were 
in  the  same  0.25-0.5  ppm  range  as  the  potential  NO2  standards,  was  one  of 
the  principal  motivating  factors  for  the  DCA  experiment.  Rather  than  simply 
accumulating  more  N0/N0x  data,  the  placement  of  monitors  and  data  recording 
rate  at  DCA  were  chosen  so  as  to  enable  separation  of  the  aircraft  contribution 
(i.e.  in  the  form  of  short  pulses  associated  with  takeof f/landing)  from 
continuous  source  and  background  contributions.  Such  a  resolution  of  aircraft 
from  non-aircraft,  sources  was  considered  vital  to  the  assessment  of  the 
aircraft  impact  on  the  NO2  standard  since  previous  monitoring  [Lorang, 

(1978)]  identified  NO2  levels  of  0.3  ppm  without  such  a  separation  while  AVAP 
modeling,  unable  to  separately  predict  NO2  levels,  indicates  that  under  worst 
case  conditions,  aircraft  contribute  N0X  concentrations  of  the  order  of  1  ppm. 

Statistical  summaries  of  the  hourly  average  concentrations,  as  computed 
from  the  high  sampling  rate  data,  are  given  in  Tables  3. 3-3. 5  for  N0X,  NO,  and 
NO2  (**  NOX~NO)  respectively.  The  fact  that  the  highest  observed  values  of  NO 
and  N0X  saturate  the  recording  equipment  and  are  outside  the  calibrated  range 
of  the  N0/N0x  instruments  is  indeed  unfortunate  and  casts  some  doubt  upon  the 
validity  of  the  NO2  data  computed  by  subtraction  of  the  NO  from  the  NOx  con¬ 
centrations  . 

Figures  3.12-3.14  show  the  cumulative  frequency  distributions  of  concen¬ 
trations  for  N0X,  NO,  and  NO2  respectively.  Aside  from  slightly  lower  NO  and 
N0X  values  at  station  1,  one  notes  a  striking  similarity  between  the  distribu¬ 
tions  from  the  different  stations.  Examination  of  these  plots  indicates  that 
95%  of  the  time  concentrations  of  N0X,  NO,  and  NO2  are  less  than  0.2  ppm, 

0.1,  and  0.07  ppm  respectively. 

Interestingly,  the  hours  corresponding  to  saturations  of  the  recorders 
for  NO  and  N0X  are  the  same  hours  of  the  CO  episode.  The  fact  that  NO  and 
N0X  from  station  1  were  not  recorded  during  this  episode  period  (due  to  the 
severing  of  the  signal  lines  by  a  snowplow)  accounts  for  the  lowered  distribu¬ 
tion  in  these  highest  percentile  ranges.  Further,  the  fact  that  this  pollu¬ 
tion  episode  affected  stations  3  and  5,  in  addition  to  1  and  6,  tends  to 
further  confirm  that  the  episode  covered  a  wider  area  than  could  be  inferred 
from  the  CO  data. 


64 


Table  3.3.  DCA  Monitoring  Experiment  N0X 


Station  1 

Station  4 

Station  3 

Station  5 

Mean  (ppb) 

54.79 

62.16 

58.63 

63.29 

Standard  Deviation  (ppb) 

50.54 

74.81 

69.62 

74.64 

Minimum  (ppb) 

2.50 

5.80 

7.00 

6.40 

Maximum  (ppb) 

433.00 

559.20s 

551  .90s 

534.20s 

Number  of  Hours 

683 

708 

683 

704 

Geometric  Mean  (ppb) 

38.80 

42.58 

41.24 

43.31 

Geometric  Standard  Deviation 

2.35 

2.23 

2.15 

2.24 

Number  of  Values  <  Threshold* 

7 

0 

0 

0 

*Values  <  instrument  threshold  (5  ppb)  included  as  1/2  threshold 
sindicates  instrument  saturation 
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Table  3 .4 . 

DCA 

Monitoring  Experiment 

NO 

Station  1 

Station  3 

Mean  (ppb) 

25.06 

29.26 

Standard  Deviation  (ppb) 

38.92 

62  .70 

Minimum  (ppb) 

2.50 

2.50 

Maximum  (ppb) 

360.50 

549.20s 

Number  of  Hours 

68  3 

683 

Geometric  Mean  (ppb' 

10.85 

11.32 

Geometric  Standard  Deviation 

3.66 

3.62 

Number  of  Values  <  Threshold* 

227 

196 

*Values  <  instrument  threshold  (5 

sindicates  instrument  saturation 

ppb)  included  as  1/2  threshold 

Table  3.5. 

DCA  Monitoring  Experiment 

N02 

Station  1 

Station  3 

Mean  (ppb) 

30.22 

29.38 

Standard  Deviation  (ppb) 

17.39 

15.04 

Minimum  (ppb) 

2.5 

2.7 

Maximum  (ppb) 

91  .90 

86.70 

Number  of  Hours 

679 

683 

Geometric  Mean  (ppb) 

24.84 

25.49 

Geometric  Standard  Deviation 

1  .97 

1  .74 

Number  of  Values  <  Threshold* 

- 

- 

*Values  <  instrument  threshold  (5  ppb)  included  as  1/2  threshold 


Cumulative  Frequency  Distributions  of  Hourly  NO 
Concentrations  at  Stations  1  and  3. 

Dashed  lines  indicate  instrument  saturation. 
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In  an  analogous  manner  to  the  method  used  for  the  CO  analysis,  Figures 
7.15-7.20  represent  an  attempt  to  estimate  the  aircraft  contributions  to  the 
distributions  given  by  Figs.  3.12-3.14.  Figures  3.15-3.17  represent  the 
"background"  subtracted  concentrations,  where,  in  lieu  of  continuous  station  2 
observations,  "background"  is  defined  as  the  average  of  the  12  minimum  concen¬ 
trations  observed  during  each  5-min  period  making  up  the  hour.  Figures 
7.18-7.20  represent  the  average  concentrations  contributed  by  pollution 
pulses  rising  at  least  0.035  ppm  above  a  15-min  average  "background"  and 
subsequently  corrected  upward  to  compensate  for  this  0.035  ppm  "barrier." 

These  curves  thus  isolate  only  the  nearby  transient  sources  of  pollution. 

As  with  CO,  extrapolation  of  the  curves  in  Figs.  3.12-3.20  out  to  the 
99.99%  probabi’ity  level  (i.e.,  1  hr  in  10,000  or  approximately  once  per 
year),  as  is  presented  in  Table  3.6,  should  be  viewed  with  some  caution  as 
such  an  extrapolation  is  based  on  only  1/12  of  a  year's  data  and  all  these 
data  come  from  a  single  contiguous  set  of  hours  rather  than  from  a  random 
selection  of  hours  throughout  the  year.  That  the  sensitivity  of  such  an 
extrapolation  to  extreme  measured  values  is  very  high  is  evidenced  by  the  fact 
that  elimination  of  the  one  severe  episode  period  lowers  the  estimated  "worst 
case"  values  of  NO  and  N0X  from  ^4  ppm  to  ^  1  ppm. 

Though  observed  values  of  background  subtracted  and  integrated  pulse 
concentrations  do  not  exceed  0.2  ppm  for  NO  or  N0X  and  0.05  ppm  for  NO2 ,  the 
extrapolations  to  99.9?%  probability,  present  in  Table  3.6,  while  highly 
speculative  in  nature,  generate  some  cause  for  concern  as  NO2  concentrations 
of  order  0.3  ppm  would  be  very  close  to  potential  hourly  average  NO2  standards 
being  investigated  by  the  EPA . 

As  with  CO,  these  N0X  results  are  consistent  with  current  "worst  case" 
model  predictions  for  LAX,  where  a  fleet  mix  and  takeoff  emissions  comparable 
to,  though  perhaps  as  much  as  25%  higher  than,  DCA  are  assumed.  AVAP  modeling 
of  LAX  under  worst  case  activity  and  dispersion  conditions  indicates  (Figure 
3.21)  N0X  levels  exceeding  0.5  ppm  more  than  one  half  mile  from  the  end  of  the 
southern*  runway  complex  but  the  key  question  is  how  these  N0X  levels  translate 
into  NO2  levels.  The  N02/N0x  ratio  is  a  function  of  plume  dispersion  rate  and 


*even  more  widespread  high  N0X  values  are  associated  with  the  jumbo  jet 
departures  on  the  northern  runway  complex 


70 


Table  3.6.  Estimated  Highest  Hourly  per  Annum* 


Concentrations  for 
gen  ( in  ppm) 

Oxides  of 

Nitro- 

NOx 

NO 

NO  2 

Total  Concentration 

1. 0-4.0 

0. 8-4.0 

0.1-0. 3 

Background  Subtracted 

0.2-0. 4 

0. 1-0.3 

0 .05®-0 . 1 1 

Pulse  Integrated 

0.2-0. 5 

0.2-0. 3 

O.^-O.S1 

*Based  on  visual  linear  extrapolation  of  cumulative 
frequency  distributions  to  99.99%  probability.  Care 
has  been  taken  to  avoid  underestimates  caused  by  NO/ 
N0X  saturation  at  0.5  ppm.  A  range  is  given  where 
linear  extrapolation  is  not  unambiguous. 

^Station  1 

^Station  3 


transport  time,  sunlight  intensity,  and  background  levels  of  NO,  NO2,  and  O3 
and  a  reactive  plume  calculation  is  required  to  obtain  a  more  definitive 
prediction;  however,  using  simple  assumptions  regarding  the  amount  of  NO2 
emitted  directly  by  the  aircraft,  the  rate  of  NO  oxidation,  and  the  ambient  O3 
level,  it  is  reasonable  to  expect  several  tenths  of  ppm  of  NO2  at  distances 
of  possible  public  exposure.  Further  extrapolations  or  generalizations  to 
other  airports  should  be  tempered  by  the  following  considerations: 

•  DCA  is  closed  between  1 0 : 00  p.m.  and  6:00  a.m. 

These  nighttime  hours  are  associated  with  stable 
atmospheric  conditions  and  thus  potentially  poor 
pollutant  dispersion  conditions. 

•  Nearly  all  operations  at  DCA  are  by  medium  range 
jets  (e.g.,  727,  737,  DC9)  primarily  using  the  JT8D- 
17  engine.  This  engine  has  a  relatively  low  N0X 
emission  rate  at  takeoff  (~200  lbs/hr)®  compared  to 
some  other  engines  (e.g.,  670  lbs/hr  for  the  CF6-50C 
and  =  47 5-500  lbs/hr  for  the  JT9D-7  and  RE-2 1 1 -22B) ,8 

•  This  experiment  was  conducted  during  winter  months. 

A  similar  experiment  during  summer  months  could  be 
accompanied  by  higher  oxidant  levels  with  resulting 
higher  NO2  levels. 
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In  order  to  relate  observed  NO2  and  N0X  levels  with  modeled  N0X  concen¬ 
trations  it  is  essential  to  note  the  difference  in  behavior  of  the  [NO]/[NOX] 
ratios  for  ambient  air  and  for  aircraft  plumes,  before  such  plumes  have 
completely  mixed  with  the  ambient  air.  During  periods  of  good  ventilation 
(sufficient  mixing  depth  and  moderate  windspeed),  the  value  of  [NO]/[NOx]  is 
generally  near  0.4  during  the  morning  decreasing  to  0.10  or  less  during  the 
afternoon,  remaining  low  throughout  the  night.  This  temporal  decrease  in  the 
ratio  is  due  to  photochemical  processes  during  daylight  hours.  It  has  been 
observed  that  at  takeoff  thrust  more  than  95%  of  the  NOx  emitted  by  jet 
aircraft  engines  is  in  the  form  of  NO.  The  concentrations  above  background 
for  aircraft  induced  peaks  measured  at  the  sites  generally  have  a  l NO ] / [ N0X ] 
ratio  exceeding  0.8,  indicating  that  some  transformation  of  NO  to  NO2  is 
taking  place  in  the  near  field  under  the  winter  conditions  observed. 

The  actual  NO  to  NO2  plume  oxidation  rate  is  a  complex  function  of 
plume  dispersion  rate  and  transport  time,  sunlight  intensity,  and  background 
levels  of  NO,  NO2,  and  O3,  and  a  reactive  plume  calculation  is  required  to 
predict  the  [ NO ] / [ N0X ]  ratio.  Clear  evidence  for  this  oxidation  process,  on 
the  short  transport  time  scales  of  the  DCA  experiment,  is  seen  in  Figure  3.22 
Plotted  is  the  ratio  of  the  hourly  integrated-pulse  concentrations  of  NO  fc 
N0X  as  a  function  of  estimated  plume  travel  time  from  the  departing  aircraft 
to  the  receptor.  The  linear  regression  line  is  indicated  for  comparison 
purposes  only  and  has  no  theoretical  basis. 

During  periods  of  light,  variable  winds,  the  hourly  N0X  concentrations 
surpassed  0.1  ppm  on  more  than  ten  (10)  separate  occasions.  For  such  periods 
hourly  background  accounts  for  greater  than  70%  of  the  mean  total  N0X  implying 
that  most  of  the  important  sources  are  nonlocal.  The  [ NO ] / l N0X ]  ratio  of  this 
background  component  generally  is  between  0.5  and  0.8.  Several  of  these  high 
concentration  episodes  are  coincident  with  low  airport  activity  and/or  non¬ 
airport  wind  directions. 

Considering  now  the  issue  of  aircraft  impact  on  the  anneal  NO2  standard 
of  0.05  ppm,  one  notes  that  the  regression  of  observed  nollutant  levels  against 
departure  rate  should  provide  some  insight.  Figure  3.2)  shows  the  oi str ibut ions 
and  regression  lines  for  the  hourly  total,  background  subtracted,  and  pulse 
integrated  NO2  concentrations  at  stations  1  and  3  versus  aircraft  departure 
rate.  The  regression  parameters  are  then  summarized  in  Table  3.7  along  with 
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similarly  obtained  parameters  for  NO  and  N0X .  Aircraft  impacts,  as  indicated 
by  the  slope  parameter,  are  seen  to  be  substantially  smaller  at  station  3  than 
at  station  1,  but  some  of  this  may  be  due  to  the  greater  difficulty  of  separat¬ 
ing  background  and  aircraft  pulse  contributions  at  station  3.  The  fact  that 
the  intercept  paramter  is  generally  consistent  with  zero  (i.e.,  that  in  the 
limit  of  zero  aircraft  there  is  no  remaining  background  pollution)  for  pulse 
integrated  concentration  regressions  strongly  suggests  that  the  pulse  inte¬ 
grated  concentration  represents  a  good  estimate  of  the  actual  aircraft  impact. 
Also,  the  fact  that  the  three  slope  estimates  for  a  given  station  and  species 
are  consistent  with  one  another  indicates  that  the  regression  estimated  impact 
is,  as  expected,  independent  of  the  assumed  background.  Thus,  using  the 
station  1,  NOj  slope  parameter  and  the  observed  monthly  averaged  departure 
rate  of  10.1  aircraft  per  hour,  one  projects  an  annual  average  aircraft  impact 
of  about  0,005  ppm  (5  ppb) .  Even  allowing  an  additional  factor  of  two  for 
preferred,  seasonal  wind  direction  effects,  these  levels  are.  small  compared 
with  the  0.05  ppm  NAAQS .  In  addition  the  observed  0.005  ppm  impact  of  aircraft 
is  in  approximate  agreement  with  the  difference  in  NO2  seen  between  DCA  and 
other  Washington  area  locations.  Table  3.8  further  suggests  that  use  of  the 
DCA  data,  acquired  during  February,  will  lead  to,  if  anything,  a  slight 
overestimate  of  annual  NO2  at  DCA.  Hence,  one  concludes  that  at  DCA,  aircraft 
account  for  only  about  one-sixth  of  the  estimated  NO2  annual  average  of  0.030 
ppm. 


Table  3.8.  Annual  NO2  Levels  at  Various  Washington 
Area  Locations 


Annua  1 
Average 

Location  (ppm)  February/Annual 


Lewi nsvi lie*  0.028  1.2 
Seven  Corners*  0.023  1  .1 
Massey  Building*  0.020  2.6 
DCA  (Feb.  only)  0.030 


*data  of  1977  from  local  pollution  control  agencies 
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3.5  SUMMARY  AND  CONCLUSIONS  BASED  ON  THE  HOURLY  AVERAGE  CONCENTRATIONS 

Results  of  CO  monitoring  near  an  aircraft  queueing  area  at  DCA  suggest 
that  observed  worst-case,  aircraft  contributions  to  hourly  CO  concentrations 
on  the  order  of  10  ppm  are  not  unexpected  at  distances  of  500  ft  from  the 
aircraft.  These  highest  observed  concentrations  decrease  to  about  5  ppm  at  a 
distance  of  1000  ft.,  a  minimum  distance  where  public  exposure  might  normally 
be  anticipated  and  in  good  agreement  with  worst-case  modeling  results  for  LAX. 
No  violations  of  the  35  ppm  hourly  standard  were  observed  even  as  close  as  500 
ft  from  the  aircraft  and  the  single  observed  violation  of  the  8  hr  standard  is 
thought  to  be  primarily  related  to  high  observed  CO  values  throughout  the  D.C. 
area  and  augmented  by  intensive  operations  of  airport  snowplows  very  near  the 
monitoring  stations. 

Preliminary  results  of  CO  monitoring  at  DCA  suggest  that  violation  of 
the  hourly  NAAQS  CO  standard,  in  areas  accessible  to  the  general  public  and  by 
aircraft  alone,  is  highly  improbable. 

Results  of  NO,  NOx  monitoring  indicate  that  NO,  NOx,  and  NOy  concentra¬ 
tion  distributions  are  nearly  independent  of  station  location  (i.e.,  within 
the  limited  spatial  regime  of  monitoring).  Worst  case  N0X  concentrations  of 
•v  1  ppm  are  consistent  with  modeling  predictions  for  LAX. 

No  NO2  concentrations  in  excess  of  0.1  ppm  were  observed  though  a 
conservative  extrapolation  to  once  a  year  probability  yields  a  concentration 
(i.e.,  '1.0.3  ppm)  in  the  same  range  as  nossible  short-term  standards.  Regres¬ 
sions  of  the  NO2  hourly  average  data  against  aircraft  departure  rate  suggests 
that  aircraft  are  responsible  for  only  about  0.005  ppm  of  the  estimated  annual 
average  NO2  of  0.03  ppm  seen  near  the  runway  at  DCA.  Th i s  projected  annual 
average  airctaft  impact  of  0.005  ppm  is  small  compared  to  the  0.05  ppm  NAAQS 
and  is  in  agreement  with  the  concentration  differential  observed  between  DCA 
and  other  Washington  area  monitors. 

3.6  ANALYSIS  OF  SINGLE  EVENT  DATA 

The  locations  of  the  monitors  with  respect  to  takoffs  and  landings  on 
Runway  36  and  the  terminal  area  at  DCA  allow  the  impact  of  airport  operations 
to  be  measured  in  several  different  ways.  In  addition  to  the  hourly  average 
analysis,  pollution  from  take-off  events  may  be  evaluated  by  subjective 
analysis  of  single  events  or  by  use  of  an  objective  single  event  evaluation 
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method.  The  complexity  of  the  concentration  patterns  at  the  various  sites 
makes  the  formulation  of  a  completely  general  event  finding  program  most 
difficult.  The  results  presented  here  are  based  on  subjective  analysis  of 
N0X  peak  concentrations  and  dosages  due  to  single  events. 

The  monitors  are  positioned  such  that  for  a  northerly  to  easterly 
wind  direction,  the  highest  possible  instantaneous  ambient  ground  level 
concentrations  of  N0X  due  to  single  aircraft  sources  should  be  experienced. 
Within  the  measurement  range,  the  exhaust  plume  is  likely  to  behave  initally 
like  a  momentum  jet  in  a  cross  flow  as  discussed  in  Section  4.5.  After  this 
initial  phase,  the  jet  momentum  is  depleted  sufficiently  for  buoyancy  forces 
to  take  effect  and  yield  an  initial  plume  rise.  The  ambient  surface  layer 
turbulent  flow  then  becomes  the  dominant  plume  transport  and  dispersion 
mechanism  during  this  jet  to  plume  transition.  The  nature  and  magnitude  of 
the  ambient  turbulence  in  the  near  field  due  to  wakes  from  nearby  aircraft  in 
various  modes  of  operation  is  expected  to  produce  enhanced  mixing.  The 
detailed  analysis  of  this  unique  data  set  should  improve  our  knowledge  of  the 
"effective  dispersion"  in  the  initial  mixing  zone,  though  the  absence  of 
monitors  in  the  vertical  plane  imposes  clear  limitations  and  renders  some 
plume  rise/dilution  ambiguities  unreso Ivab le  . 

3.6.1  Subjectively  Analyzed  Single-Event  Data 

As  was  seen  qualitatively  from  Figure  3.3,  it  is  possible  to  extract 
information  about  the  dispersion  of  the  takeoff  plume  through  examination  of 
the  N0X  concentration  time  histories  recorded  at  the  four  N0X  instrumented 
stations.  However,  the  interpretation  of  the  single  event  peak  concentrations 
and  dosages  in  terms  of  typical  Gaussian  dispersion  parameters  is  not  completely 
straightforward.  The  parametric  Gaussian  formulation  requires  that  the 
turbulent  elements  have  temporal  and  spatial  scales  much  smaller  than  the 
observed  events  which  in  this  case  are  jet  exhaust  plumes  generated  on  takeoff. 
At  airports,  surface  layer  turbulence  is  enhanced  by  aircraft  wakes  and 
trailing  vortices.  Therefore,  especially  when  conditions  are  far  from  neutral, 
the  scale  separation  is  insufficient  to  analyze  each  event  with  a  simple 
Gaussian  model.  However,  a  large  ensemble  of  events  may  more  appropriately  be 
approached  with  this  methodology. 
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With  this  objective  in  mind,  approximately  120  individual  aircraft 
departure  events  were  measured  through  digitization  of  the  concentration  time 
histories.  Peak  concentration,  time  from  peak  noise  to  peak  concentration, 
full  width  at  half  maximum  above  background,  and  background  concentration  were 
extracted  for  subsequent  input  into  the  ensemble  research  model  for  the 
takeoff  mode,  described  in  Section  3.6.3.  The  data  extraction  technique  is 
referred  to  as  the  "subjective"  approach  primarily  because  of  the  element  of 
judgment  involved  in  the  estimation  of  background  levels. 


3.6.2  DCA  Single  Event  Finding  Program 

An  attempt  has  been  made  to  develop  an  objective  algorithm  to  extract 
peak  concentrations  and  dosages  from  the  month  of  high  sampling  rate  data 
gathered  at  DCA.  After  the  measurements  had  been  transferred  to  a  master 
archive  tape,  calibrated  in  a  preliminary  fashion,  and  edited  to  minimize 
inclusion  of  periods  of  uncertain  data,  the  following  procedure  was  adopted  to 
isolate  and  quantify  events  associated  with  individual  aircraft  departures. 


1.  Period  of  interest  is  specified  such  that  periods 
of  missing  data,  zero  airport  activity,  non-optimal 

wind  direction,  etc.  may  be  avoided.  Only  wind  directions 
between  10  and  80°  were  selected  since,  for  other  direc¬ 
tions,  the  plume  is  transported  away  from  the  monitoring 
sites  . 

2.  The  search  for  a  usable  event  occurs  as  follows:  a 
noise  pulse  is  searched  for  that  is  sufficiently  separated 
from  other  pulses  to  allow  transport  to  the  furthest 
receptor  before  the  next  aircraft's  plume  impacts  the 
closest  receptor.  An  event  pulse  is  defined  by  an  80  dB 
noise  threshold.  When  runway  18  is  in  use  the  noise 
spikes  are  similar  but  the  wind  is  from  the  south!  The 
screening  in  step  !  then  becomes  important  to  avoid  the 
possibility  of  erroneous  results.  The  transport  time  is 
defined  as  D/ u  where  D  is  the  along  wind  distance  from 
station  3  to  the  runway.  Five  minute  average  wind  speed 
(u)  and  direction  values  are  utilized.  If  the  transport 
time  exceeds  the  noise  event  separation  che  event  is 

sk  i  pped  , 

3.  Once  a  noise  event  is  identified  as  giving  rise  to  an 
analyzable  "dose  event",  3-minute  averages,  background 
values,  <’d,  and  o  are  evaluated  (1  min  prior  to  noise 
pulse,  2  min  after  noise  pulse).  Missing  data  or  data 
exceeding  0.5  ppm  will  cause  the  event  to  abort. 
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4.  The  peak  concentrations,  triangular  doses,  and 
numerically  integrated  doses  are  now  found  by 
subtracting  out  the  background  values  of  the  3 
minute  period. 

a.  For  each  receptor  the  search  for  the  event 
begins  after  a  delay  time  equal  to  the  time  of 
transport  from  the  runway.  A  window  the  size 
of  the  minimum  pulse  separation  used  in  step  2 
is  searched  for  a  maximum. 

b.  The  times  of  half-maximum  (above  background) 
are  determined,  to  calculate  the  triangular 
dose  . 

c.  The  endpoints  of  the  numerical  dose  integration 
are  defined  as  follows: 

The  starting  point  is  the  time  where  the  search 
for  the  maximum  begins.  The  end  point  the  time 
at  which  the  concentration  becomes  smaller  than 
(1.1  x  background)  or  (background  +  10  ppb) ,  wich- 
ever  is  greater,  or  the  end  of  the  window  in  (a), 
whichever  comes  first.  For  a  weak  pulse  above  a  high 
background  at  station  3  the  size  of  the  integration 
interval  may  be  underestimated.  Early  versions  use 
Simpson's  rule  but  later  attempts  employ  a  Gaussian 
quadrature  technique. 

Application  and  limitations  of  the  above  described  technique  include: 

1.  Care  must  be  taken  to  specify  periods  with  the 
proper  wind  direction  or  for  each  5  minute  period  the 
average  wind  direction  must  be  screened  to  insure 
acceptable  events  ti.e.,  10°  £  0  <_  80°  only). 

2.  Priority  should  be  given  to  the  subset  of  hours 
when  onsite  observations  are  available  because 
noise  network  printouts  of  runway  activity  are  not 
entirely  reliable. 

3.  The  dose  of  NO  sometimes  exceeds  that  of  N0„:  this 
is  most  likely  due  to  a  calibration  problem  and  not 
inherent  in  the  method.  The  same  problem  is 
evident  in  the  subjectively  analyzed  single  events. 

4.  The  specification  of  the  "window"  to  search  for  an 
event  is  subject  to  further  experimentation  and 
refinement.  A  better  way  to  find  the  endpoints  for 
the  integration  may  be  required. 

5.  The  temporal  spacing  requirements  will  cause  many  legiti¬ 
mate  events  to  be  skipped.  However,  with  the  data  base 
provided,  the  present  program  should  provide  an  adequate 
cross-section  of  "well  defined"  single  events. 
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Though  the  above  described  approach  had  yielded  a  data  base  consisting 
of  several  hundred  N0X  plumes  from  departing  aircraft,  subsequent  analysis  of 
these  objectively  obtained  events  has  indicated  the  model's  inability  to  deal 
with  numerous  pathological  c-.uuat  ions .  Research  model  results  using  these 
single  events  as  input  thus  remain  high  suspect. 

3.6.3  DCA  Research  Model  Results 

The  DCA  research  model  basically  considers  the  aircraft  to  undergo 
uniform  linear  acceleration  down  the  runway  while  emitting  pollutant  at  a 
uniform  rate  from  a  tail  of  length  L  attached  directly  behind  the  aircraft. 

The  initial  transverse  and  vertical  extent  of  the  plume  is  specified  by 
parameters  ox(0)  and  az(0)  respectively  and  subsequent  plume  growth  is 
assumed  to  be  governed  by  the  equations 

ax=oy=ox(°)  +  bxayX(x) 

and 

az  =  az(0)  +  bzCTzT<*> 

where  oyT(x)  and  azT(x)  are  the  dispersion  coefficients  taken  from  Turner's 
Workbook  and  bx  =  bz  =  0.7  is  a  correction  for  averaging  time.  Unfortunately, 
the  fact  that  all  receptors  are  near  ground-level  prevents  independent  determina¬ 
tion  of  plume  center  line  height  (assumed  to  be  at  H0  E  1.2  az( 0)),  so  that 
effects  of  plume  rise  are  erroneously  incorporated  into  the  available  parame  i_ers . 

Minimization  of  the  chisquare  between  predicted  and  observed  peak 
concentrations  and  pulse  durations  yields  the  parameters 

3X(0)  =  21.1  m, 

cz( 0)  =  25.7  m, 

and 

L  =  52  m 

One  notes  that  the  tail  length  parameter  L  is  unreasonably  small  and  probably 
is  related  to  simplifying  assumptions  made  in  this  research  model.  Never¬ 
theless,  the  research  model  algorithm  matcl  es  the  AVAP  line  source  model  very 
well  so  that  reasonable  concentration  predictions  should  be  expected  from  AVAP 
despite  the  small  tail  length. 
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Figure  3.24  shows  a  scatter  plot  of  observed  and  research  model  predict¬ 
ed  peak  concentrations.  Because  of  an  approximation  which  assumes  that  the 
peak  concentration  is  primarily  determined  by  pollution  coming  from  the  region 
of  the  intersection  of  the  wind  dirction  and  the  runway,  the  full  range  of 
peak  concentrations  is  not  reproduced.  An  approach  is  currently  under  develop¬ 
ment  that  will  avoid  this  problem  while  maintaining  maximal  overlap  with  the 
AVAP  method  so  that  research  model  determined  paramters  may  confidently  be 
used  in  AVAP.  However,  it  should  be  pointed  out  that  the  present  model  does 
indicate  the  gross  error  associated  with  neglect  of  the  initial  turbulent 
mixing:  assumption  of  an  initial  accelerating  point  source  leads  to  a  theory/ 

observation  regression  slope  of  ~30. 
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4  ANALYSIS  OF  EXPERIMENTS  AT  DULLES  INTERNATIONAL  AIRPORT 

4 . 1  INTRODUCTION 

The  measurement  program  at  Dulles  International  Airport  was  initiated 
in  1976  in  response  to  an  order  by  the  Secretary  of  Transportation  to  monitor 
pollutant  emissions  and  noise  levels  associated  with  Concorde  aircraft  opera¬ 
tions.  Three  pollutants  for  which  there  are  engine  emission  standards  (carbon 
monoxides  (CO) ,  nitrogen  oxides  (N0X) ,  and  hydrocarbons  (HC) )  were  monitored. 
Measurements  were  obtained  at  nearby  regional  monitoring  stations,  and  these 
data  were  analyzed  using  statistical  inference  techniques  as  well  as  by  means 
of  a  source  finding  algorithm;  a  method  designed  to  locate  sources  and  assess 
source  culpability  based  on  observed  concentrations. 

Principal  effort,  both  during  the  experimental  program  and  in  subse¬ 
quent  analysis  efforts,  was  devoted  to  arrays  of  sensors  placed  at  the  airport 
in  the  vicinity  of  aircraft  operations.  The  locations  and  periods  of  use  for 
these  arrays  are  detailed  elsewhere  by  Smith  et  al .  (1977) 

Early  measurements  indicated  low  concentrations  that  would  not  be 
explicable  with  many  conventional  airport  models,  apparently  because  previous 
models  developed  specifically  for  airports  have  generally  ignored  plume  rise 
and  initial  plume  dilution.  Although  many  model  applications  are  not  severe¬ 
ly  limited  by  this  omission,  when  conservative  estimates  of  airport  impact  at 
distances  of  2  to  10  miles  are  at  issue,  it  is  essential  to  consider  both  the 
direct  and  indirect  (augmentation  of  oz)  effects  of  plume  rise  for  validating 
model  predictions  at  closer  distances.  Thus,  the  early  results  led  to  a 
design  of  a  progressively  more  sophisticated  experiment,  using  first  one,  then 
two,  and  finally  three  towers  instrumented  at  three  to  five  levels  with  CO 
sensor  probes.  Meteorological  data  included  two  levels  of  wind  direction  and 
speed  and  temperature  and  its  gradient.  These  data  have  been  analyzed  in 
detail  to  provide  information  on  jet  plume  rise,  actual  atmospheric  dispersion 
parameters,  and  vertical  and  horizontal  "profiles"  of  exhaust-plume  pollutant 
concentrations  for  individual  aircraft  in  actual  service. 
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4.2  REGIONAL  MEASUREMENTS 

As  previously  mentioned,  the  monitoring  program  at  Dulles  initially 
consisted  of  a  network  of  stations  located  in  and  around  the  airport.  These 
stations  recorded  the  concentrations  of  various  pollutants  and  relevant 
meteorological  data  during  the  time  period  from  May  to  September,  1976. 

The  locations  of  these  stations  are  shown  in  Fig.  4.1,  and  the  parame¬ 
ters  measured  at  each  site  are  tabulated  in  Fig.  4.2.  Two  stations,  South 
Ramp  and  Sterling  Park,  were  especially  important  to  the  monitoring  program 
because  of  their  proximity  to  the  airport  and  the  amount  of  data  recorded 
at  each  site.  The  South  Ramp  site  was  located  approximately  200  ft  south  of 
the  jet  taxi  ramp,  and  the  Sterling  Park  station  was  located  3  miles  north- 
northeast  of  Dulles  in  the  community  of  Sterling  Park. 

The  data  at  these  two  sites  was  recorded  on  strip  charts  which  were 
later  digitized  and  processed  by  computer  to  yield  hourly  averaged  concentra¬ 
tion  and  meteorological  values.  Most  of  the  data  in  the  Dulles  regional  data 
base  consists  of  measurements  made  at  those  two  sites. 

4.3  SOURCE  FINDING  APPLICATION  AT  DULLES 

4.3.1  Introduction 

The  increasing  number  and  complexity  of  air  quality  monitoring  networks 
in  operation  throughout  the  world  coupled  with  the  greatly  increased  quantity 
of  data  such  networks  generate,,  suggests  the  need  for  innovative  data  analy¬ 
sis  techniques. 

Many  analysis  tools  (such  as  the  use  of  frequency  and  cumulative 
frequency  distributions,  time  series  analysis,  and  factor  analysis),  while 
well  suited  for  understanding  the  significance  of  data  from  a  single  receptor, 
become  somewhat  unwieldy  when  applied  to  data  from  a  network  of  receptors. 
Ideally,  one  might  hope  that  the  data  from  an  n-station  array  would  result  in 
more,  and  not  less,  information  than  obtainable  from  n-stations  considered 
ind ividval ly . 

The  capability  of  utilizing  observed  meteorological  and  pollutant 
concentration  data  to  determine  the  location  and  strength  of  those  sources 
having  the  greatest  influence  on  measured  concentrations  is  su^h  a  tool, 
having  potential  applications  which  include: 
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a)  an  independent  means  of  checking  source  inventories 
and  determining  source  culpability, 

b)  the  ability  to  locate  inadvertent  leaks  of  hazard¬ 
ous  effluents  into  the  atmosphere,  and 

c)  as  ar,  enforcement  tool  for  air  pollution  regulatory 
agencies . 

One  possible  approach  to  this  problem  is  to  develop  a  technique  for 
relating  the  aerometric  data  obtained  from  a  monitoring  network  under  a  wide 
range  of  meteorological  conditions  to  the  spatial  distribution  and  strengths 
of  emissions  sources  via  the  use  of  existing  air  quality  dispersion  algo¬ 
rithms.  The  approach  may  be  envisioned,  quite  literally,  as  the  running  of  an 
air  quality  dispersion  model  "in  reverse";  where  an  effective  emission  density 
map  is  determined  from  knowledge  of  the  concentrations  at  the  receptors  of  a 
network  and  the  values  of  a  few  relevant  meteorological  parameters  over  the 
area  of  interest. 

The  regional  data  base  available  from  the  Dulles  monitoring  program 
can,  therefore,  be  used  as  input  to  this  source  finding  algorithm  in  an 
attempt  to  determine  the  significance  of  local  or  on-site  sources  on  the 
overall  measured  concentrations, 

4.3.2  Model  Development 

Air  quality  dispersion  models  are  most  often  concerned  with  the  deter¬ 
mination  of  pollutant  concentrations  at  a  receptor  given  known  source  streng¬ 
ths  and  locations.  Assuming  steady  state  conditions  have  been  achieved,  the 
concentration,  ,  at  the  k^  receptor  during  the  t^  time  interval  can  be 
expressed  as 

C^t  ~  ^  Rjkt  Qj  » 

j 

where  Qj  is  the  time  independent  strength  of  the  j1-*1  source  and  Rjkt  is 
the  transport  coupling  coefficient  between  the  source  and  kc^  receptor  for 

the  meteorological  conditions  existent  during  the  t1-^  time  interval.  The 
solution  of  the  simple  inverse  problem 

Qj  =  I  ^jkt  ckt 

is  not  of  particular  interest  as  a  unique  solution  might  exist  only  if  the 
number  of  receptors,  K,  equalled  the  number  of  scarce  candidates,  J.  A 
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more  reasonable  problem  would  consist  of  determining  the  set  of  Qj  values 
leading  to  the  minimization  of  the  quantity  y^- ,  defined  as 


where  ACkt;  is  the  uncertainty  in  Ckt ,  i  is  a  dummy  index,  and  a  summation 
over  both  t  and  k  is  required.  This  problem  will  generally  lead  to  a  unique 
set  of  J  source  strengths  provided  the  number  of  measurements  M  (nominally  KT, 
where  T  is  the  number  of  time  periods  of  data  available)  exceeds  J,  and 
secondly,  that  the  sources  Qj  of  interest  actually  couple  to  the  data  in  hand 
(i.e.,  Rjkt  ?  0  for  all  j  and  some  k,  t).  Given  these  conditions,  one  may 
write  down  the  set  of  J  equations  generated  by  the  relations 


3Qj 

These  equations  are  of  the  form: 

l  Aij  Qj  *  Bi  (2) 

j 

where 

V  Rikt  Rjkt 

Ajj  =  L  - 

t,k  (A  t ^ 

and 


r,  ckt  Rikt 

Bi  =  l - r~ 

t.k  (ACkt)z 

and  upon  obtaining  the  inverse  A-i  of  the  positive  definite  matrix  A,  one 
arrives  at  the  solution 

Qj  -  l  Aj-  Bi  (3) 


This  approach,  with  no  constraint  on  the  Qj ,  is  quite  acceptable 
provided  the  number  of  candidate  source  locations  J  is  small  (e.g.,  J<100) 
and  source  candidate  locations  are  well  separated;  however,  if  nothing  is 
known  about  the  source  locations,  and  instead,  an  array  of  point  or  area 
sources  is  conjectured  to  exist  on  an  X-Y  grid  of  size  n  x  in,  then  one  must 
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contend  with  a  matrix  A  of  order  (n  x  m)  by  (n  x  m)  ,  which,  in  turn  demands 
a  minimum  of  n  *  ra  concentration  measurements.  In  addition,  the  solution  Qj 
may  be  quite  unstable,  with  adjacent  grid  point  Qj  values  oscillating  wildly 
positive  and  negative  in  a  region  known  to  be  source  free.  Fortunately  the 
availability  of  a  non-negative  least  squares  algorithm,  developed  by  Lawson 
and  Hanson  (1974),  enables  one  to  solve  equation  (2)  subject  to  the  physically 
reasonable  requirement  that  all  source  strengths  Qj  be  positive  (i.e., 

Qj  >  0).  Further  theoretical  developments  are  presented  by  Yamartino  and 
Latnich  (1979). 


4.3.3  Model  Application 

In  order  to  apply  the  source  finding  algorithm  to  the  regional  Dulles 
data,  the  coupling  coefficients,  and  a  source  grid  must  be  defined. 

Computation  of  Rjkt*  The  coupling  coefficients,  Rjkt,  were  computed  via  the 
Gaussian  plume  formalism;  however,  due  to  the  apriori  unknown  nature  of  the 
sources,  the  calculations  were  performed  as  if  each  source  were  a  ground- 
level,  square,  area  source  of  length  AX  (AX  =  320  meters  -  0.2  miles)  so  that 
the  array  of  source  candidates  is  contiguous.  Plume  rise,  which  also  must  be 
assumed  to  be  unknown  in  any  "source-finding"  program,  has  been  neglected, 
and  it  should  be  noted  that  this  highly  source  oriented,  yet  meteorologically 
dependent,  variable  represents  the  major  stumbling  block  in  the  otherwise 
trivial  extension  of  this  two-dimensional  "source-finding"  theory  to  three 
dimensions.  This  neglect  of  plume  rise  limits  the  usefulness  of  this  approach 
with  the  Dulles  data  since  the  South  Ramp  station  is  close  enough  to  the 
aircraft  for  concentrations  there  to  be  strongly  affected  by  jet  plume  rise. 

Computational  Grid  and  Receptors.  In  Fig.  4.3  the  computational  grid  used  by 
the  source  finding  algorithm  is  shown  superimposed  on  a  map  of  the  airport. 
Also  shown  are  the  locations  of  the  two  receptors  which  gathered  the  concen¬ 
tration  data. 

Only  two  of  the  receptors  from  the  regional  monitoring  network  were 
used  initially  because  of  the  steady-state  requirement  in  the  original  source 
finding  algorithm.  The  assumption  of  steady-state  conditions  in  the  computa¬ 
tion  of  the  coupling  coefficients  (Rj^t)  means  that  the  sources  must  couple  to 
each  receptor  during  the  time  interval  that  the  computations  are  being  done. 


Hence,  there  is  a  cutoff  distance  around  the  airport;  any  receptors  lying 
outside  this  area  cannot  be  used  because  t'.,£  pollutants  emitted  from  air¬ 
port  sources  during  some  time  interval  will  not  arrive  at  those  receptors 
until  a  later  time  interval.  The  only  two  receptors  that  were  close  enough  to 
Dulles  to  be  suitable  for  such  source  finding  application  were  the  Sterling 
Park  and  South  Ramp  stations.  Taking  into  fccount  network  sensitivity 
and  resolving  power,  it  was  expected  that  using  only  two  receptors  would  not 
result  in  a  very  accurate  picture  of  the  significant  sources.  A  subsequent 
version  of  the  source  finding  model  circumvented  this  limitation  through  the 
use  of  backward  trajectories,  which  allowed  for  up  to  six  hours  of  pollutant 
transport . 

Each  mesh  point  of  the  computational  grid  is  treated  by  the  algorithm 
as  a  potential  area  source.  Minimization  of  the  z  value  produces  an  optimal 
set  of  relative  source  strengths,  where  each  member  of  this  set  corresponds  to 
an  area  source  located  on  the  grid.  When  each  relative  source  strength  is 
plotted  on  a  map  in  its  proper  location,  the  resulting  picture  ideally  will 
show  the  locations  and  strengths  of  the  actual  sources  inside  the  area 
covered  by  the  grid . 

4.3.4  Model  Results 

Contours  proportional  to  logjg  of  CO  emission  density  are  plotted  in 
Fig.  4.4.  Despite  the  use  of  additional  data  through  application  of  the 
trajectory  method,  the  solution  is  nearly  identical  with  that  obtained  (not 
shown)  using  only  the  Sterling  Park  and  South  Ramp  data.  In  essence  these  are 
the  only  stations  close  enough  to  the  airport  to  be  influenced  by  it. 

The  presence  of  strong  sources  at  the  edges  of  the  computational  grid 
is  an  indication  that  much  of  the  observed  CO  actually  originates  from  points 
outside  of  the  grid.  One  notes  further  that  the  strength  of  these  outside 
sources  is  as  much  as  two  orders  of  magnitude  greater  than  CO  emission 
densities  found  on  the  airport,  and,  in  particular,  in  the  vicinity  of  the 
South  Ramp  where  a  major  fraction  of  the  aircraft  CO  would  originate. 

4.3.5  Conclusions 

A  formalism  has  been  developed  for  utilizing  data  obtained  from  an  air 
quality  monitoring  network  to  determine  the  locations  and  strengths  of  those 
sources  contributing  most  significantly  to  the  observed  concentrations.  The 
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problem  is  best  treated  exactly  for  the  case  of  J  source  candidates  via 
solution  of  the  coupled  equations 

l  Aij  Qj  =  Bi 

i 

subject  to  the  condition  Qj  >  0  for  all  j. 

Application  of  this  source  finding  algorithm  to  the  CO  data  obtained 
from  the  regional  network  at  and  around  Dulles  International  Airport,  correct¬ 
ly  locates  the  leading  aircraft  CO  emission  zone  but  further  indicates  the 
presence  of  substantially  stronger  or f-airport  sources. 

4.4  THREE  TOWER  MEASUREMENTS  OF  CO 

The  Concorde  sir  quality  monitoring  and  analysis  program  conducted 
at  Dulles  International  Airport  during  1 9  7  6  —  7  7  provided  a  unique  opportunity 
to  measure  CO  plumes  from  taxiing  aircraft.  The  transport  and  dispersion  of 
these  CO  plumes  was  monitored  at  13  points  on  the  three  tower  array  shown  in 
Fig.  4.5,  measured  with  Ecolyzers,  and  recorded  on  high-speed  strip  chart 
recorders.  CO  values  were  extracted  from  measurements  of  strip  chart  records. 
A  sample  set  of  CO  traces  is  also  seen  in  Fig.  4.5.  Wind  speed  and  direction 
were  measured  at  the  80'  and  14'  levels  on  the  first  tower.  Temperature 
gradient  was  measured  between  67'  and  14'  on  the  same  tower.  Several  hundred 
plumes  were  observed  under  neutral/ uns tab le  daytime  conditions  during  the 
one-,  two-,  and  three-tower  phases  of  the  experiment.  Commercial  aircraft 
types  monitored  included  the  Concorde,  707,  727,  737,  747,  DC8,  DC9,  DC10, 
and  L1011.  Though  peak  instantaneous  CO  levels  reached  10  ppm  at  the  first 
tower  (only  215  ft  from  the  taxiway  centerline),  maximum  aircraft  contribution 
to  the  hourly  average,  ground  level,  CO  concentration  remained  below  0.06  ppm 
per  aircraft.  Extrapolations  to  1000  ft  from  the  taxiway  indicate  a  maximum 
hourly  average  concentration  of  0.03  ppm  CO  per  aircraft.  Thus,  hourly 
concentrations  in  excess  of  several  ppm,  adjacent  to  a  busy  taxiway,  would 
be  unlikely  to  occur. 

Though  maximum  CO  impacts  are  expected  from  queueing  operations  rather 
than  taxi,  the  data  from  the  taxi  mode  provide  interesting  information  on 
initial  plume  dimensions  and  buoyant  plume  rise.  These  plume  parameters  may 
then  be  used  in  airport  air  quality  models  to  increase  their  accuracy  and 
predictive  power. 
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As  indicated  in  Fig.  4.3,  CO  was  measured  by  pumping  air  samples 
continuously  trough  identical  volume  sampling  line3  into  individual  Ecolyzer 
units.  These  nits  were  housed  in  an  air  conditioned  shelter.  They  were 
periodically  calibrated  by  sequential  switching  of  the  intakes  to  the  same  18 
ppm  concentration.  The  calibration  system  was  designed  to  allow  precise 
timing  of  sensor  exposure  to  calibration  gases  of  different  concentrations  so 
that  response  time  constraints  and  linearity  of  signal  amplitude  could  be 
determined.  Since  an  aircraft  passage  "event"  was  expected  to  produce  a  pulse 
representing  concentration  versus  time  (as  shown  if  Fig.  4.6),  the  measurement 
of  sensor  system  time  characteristics  was  deemed  important.  The  time  constant 
of  the  Ecolyzers  averaged  12  seconds,  and  their  threshold  sensitivity  averaged 
0 . 25  ppm . 

The  concentration  shown  in  Fig.  4.5  represent  the  instantaneous  peak 
values  from  the  relatively  high  speed  chart  records.  Figure  4.6  is  idealized 
in  the  sense  that  the  skew  (to  the  right)  observed  as  a  result  of  the  time 
response,  and  potentially  the  pollutant  distribution,  is  not  shown.  When 
the  data  was  reduced,  both  the  time-to-peak  and  the  t ime-to-half-peak  were 
recorded  in  addition  to  the  peak  CO  value,  the  f ul 1-width-at-half-maximum 
time,  and  the  background  CO  so  that  skewness  could  be  accounted  for  in  future 
modeling.  Details  of  the  method  for  correcting  the  peak  concentrations  when 
one  uses  an  event  modeling  technique  arc  given  by  Smith  (1977). 

In  addition  to  concentration  measurements,  documentation  for  each  event 
included  event  time,  direction  of  aircraft  travel,  departure  or  arrival  made, 
time  to  travel  50  m,  and  meteorological  conditions.  Wind  direction  and  speed 
were  averaged  over  three  minutes.  The  value  of  for  the  same  averaging  time 
was  found  from  30  six-second  samples  for  each  event,  commencing  at  the  record¬ 
ed  event  time.  The  specific  ranges  of  meteorological  conditions  are  given 
with  the  results  b-iow.  All  selected  tests  were  conducted  in  the  daytime  and 
had  winds  between  290°  and  70°,  The  taxiing  activity  pattern,  the  orientation 
of  the  towers  at  Dulles,  and  simplicity  for  modeling  were  factors  leading  to 
this  selection. 
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4.5  PLUME  RISE  FROM  JET  AIRCRAFT  DURING  THE  TAXI  MODE 

4.5.1  Introduct ion 

This  section  discusses  the  results  of  an  investigation  of  the  behavior 
of  buoyant  jet  engine  exhaust  plumes  in  a  crosswind;  it  attempts  to  identify 
the  degree  to  which  the  plume  rise  can  be  described  by  relationships  developed 
for  other  types  of  sources. 

At  least  four  factors  affecting  the  rate  of  dilution  of  jet  exhaust 
before  it  reaches  receptor  adjacent  to  taxi-ways  or  runways  have  been  previ¬ 
ously  identified  [Heywood  et  al .  (1971)]: 

1.  turbulent  mixing  of  the  jet  exhaust  at  the  engine 
exit 

2.  buoyant  plume  rise 

3.  advective  dilution 

4.  dispersion  by  ambient  turbulence 

Observational  studies  of  plumes  generally  allow  only  one  or  two  sepa¬ 
rate  processes  to  be  measured  (plume  rise  and  total  dispersion  rate).  Al¬ 
though  the  bending  of  an  exhaust  plume  from  its  original  release  axis  until  it 
is  aligned  with  the  prevailing  wind  direction  is  also  observable  (particularly 
from  above) .  This  change  in  orientation  or  "bending"  may  also  be  viewed  as 
the  transition  from  plume  dilution  dominated  by  the  first  mechanism  to  dilu¬ 
tion  controlled  by  the  latter  three.  For  this  reason,  the  maximum  length  of 
the  highly  turbulent  jet  trail  as  a  function  of  wind  speed  is  of  inter- 
es  t . 

The  assumption  is  made  that  the  two  phenomena,  plume  bending  and  plume 
rise,  can  be  treated  independently  as  a  first  approach.  Both  theoretical  and 
empirical  models  are  available  to  describe  plume  bending  in  a  perpendicular 
wind  [e.g.  Abramovich  (1963)].  Estimates  of  the  maximum  distances  of  domi¬ 
nance  of  jet  exhaust  mechanical  turbulence  are  made  for  taxiing  aircraft. 

The  estimates  here  are  restricted  to  perpendicular  winds  for  simplicity. 

Analysis  of  the  experimental  data  revealed  that  the  precision  of 
measurement  of  the  initiation  time  of  each  aircraft  passage  "event"  was  not 
adequate  for  analysis  of  differences  between  expected  arrival  times  for  CO  at 
the  first  tower  under  alternate  plume  bending  hypothesis.  Thus,  although 
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these  alternative  descriptions  of  aircraft  plume  bending  are  given,  the 
present  comparisons  with  experimental  evidence  are  restricted  to  the  phenome¬ 
non  of  plume  rise.  It  is  this  mechanism  for  aircraft  plume  dilution  that  was 
of  primary  concern  in  the  Dulles  experiments  [Smith  (1977)],  although  the 
other  three  mechanisms  listed  were  also  considered. 

4.5.2  Modeling  Turbulent  Jet  Exhausts  without  Plume  Rise 

Aircraft  jet  exhausts  discharged  horizontally  into  a  uniform  crosswind 
may  be  described  •’  i  two  stages:  the  momentum-dominant  stage  and  the  buoyancy- 
dominant.  stage.  In  the  momentum-dominant  stage,  the  horizontal  velocity  of 
the  jet  plume  decays  through  turbulent  mixing  with  ambient  air,  and  plume  rise 
is  suppressed.  1'  the  buoyancy-dominant  stage,  the  plume  rises  and  is 
entrained  by  the  vertical  motion.  If  it  is  assumed  that  there  are  no  interac¬ 
tions  between  adjacent  engine  plumes  and  plume  rise  ignored,  the  benching  path 
of  a  nonbuoyant  momentum  jet  in  a  crosswind  may  be  estimated  from  Eq .  (1): 

y  =  1.5  Vr2/3  x  1/3  (1) 

where  the  coordinate  system  is  that  shown  in  Fig.  4.7,  with  B  (angle  between  the 
y-axis  and  aircraft  path)  equal  to  zero  [see  Briggs  (1969)]. 

Vr  =  Ve/u 

Ve  »  exhaust  velocity 
u  =  windspeed 

y  =  y/° 

x  x/D 

D  =  exit  diameter 

the  maximum  penetration  length  of  the  exhaust  plume  behind  the  aircraft  can  be 
estimated  from: 

ymax  ”  3Vr  and  xmax  =  8Vr  (2) 

At  this  ymax  distance  behind  the  aircraft,  the  angle  between  the  plume  center- 
line  and  the  y-axis  may  be  found  from  Eq .  (3): 


0  =  tan  (dy/dx)  =  tan~l  1/2  ^ ^2/3 


therefore,  r-. 


Fig,  4.7.  Coordinate  System  for  Model 
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To  evaluate  ymax  and  x^x  for  engines  of  the  specific  aircraft,  Table  4.1 
should  be  consulted.  Also  presented  in  Table  4.1  are  appropriate  values  of 
the  exit  velocities  and  temperatures  for  calculation  of  effective  velocity 
ratio: 


V'  3 


'  V  2 
P  e 


pu' 


1/2 


where  p  =  density  of  the  ambient  air 
P!  =  density  of  the  jet  exhaust 


(4) 


It  is  expected  that  substitution  of  in  the  Eq .  (1)  through  (3)  will  yield 
more  accurate  estimates  for  jet  exhausts.  In  Table  4.1,  the  exhaust  diameter, 
exit  velocity,  and  exit  temperature  is  given  for  the  JT-3  and  JT-8  engine 
during  taxi/idle  mode  operation.  Thrust  values  and  mass  emission  rates  are 
also  given  for  comparison.  For  average  surface  winds  of  3  m/sec,  Vr  would 
range  from  15  for  the  JT3s  to  23  for  the  JT-8s  during  taxiing  operations. 

For  the  range  of  8  <  V^.  <  54  and  34,  experimental  evidence  [Patrick 

(1967)]  indicates: 

y  =  Wr)  (x)  and  ymax  -  2.3  (Vr)  (5) 


These  relationships  yield  similar  results  to  those  obtained  from  Eqs. 
(1)  and  (2).  Thus,  for  taxiing  B70/s,  ymax  -  40  m  and  for  B7.?7s  ymax  53  m, 
and  corresponding  xmax  values  of  108  m  and  138  m  (with  wind  speeds  of  5 
m/sec).  For  sensors  near  the  edge  of  the  taxiway,  the  value  of  xmax  would 
determine  whether  dilution  of  the  plume  reaching  those  sensors  was  dominated 
by  jet  trail  turbulence  or  ambient  turbulence. 


Table  4.1.  Aircraft  Engine  Emission  Parameters 


Aircraft  Type 

B707 

B727 

Engine  Type 

JT-3 

JT-8 

Diameter,  Exit  (m) 

0.9 

0.75 

Ve  (m/sec) 

76 

1  14 

Te  (°K) 

38  6 

440 

Mass  Rate  (Kg/sec) 

45 

40 

Thrust  (Nt) 

270 

250 

Source:  Goldberg  (1978) 
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4.5.3  Plume  Rise  Modeling 


To  obtain  a  simple  plume  rise  equation  for  a  buoyant  plume  it  is 
necessary  to  make  some  basic  assumptions: 

1.  The  flow  is  fully  turbulent,  thus  the  effect  of 
molecular  viscosity  or  Reynolds  number  is 
negligible. 

2.  Boussinesq  approximation  is  valid,  i.e.,  local  den¬ 
sity  variations  are  neglected  except  when  multi¬ 
plying  by  gravity. 

3.  The  buoyancy  is  assumed  to  be  conserved. 

4.  The  fluids  are  quasi-incompressible. 

This  theory  or  the  2/3  power  relation  was  obtained  by  Slawson  and 
Csanady  (1967)  and  substantiated  by  Briggs  (1969)  and  Koult,  Fay  and  Forney 
(1969)  . 


Using  the  entrainment  hypothesis  given  by  Morton,  Taylor,  and  Turner 
(1956),  one  may  express  the  rise  of  the  buoyant  plume  from  jet  aircraft 


as : 


z 


1/3 


F  1/3  u'1  (x 


where 


(6) 


z,.,  -■  initial  height 

x0  =  initial  downwind  distance 

F  =  bouyancy  flux 

Q  =  entrainment  constant 

u  =  windspeed 

Although  this  2/3  power  law  relation  was  developed  originally  for 
stationary  sources  with  lower  exit  plumes  than  jet  exhaust,  its  use  as  a  first 
estimate  for  the  present  application  is  encouraged  two  factors: 


1.  The  heat  flux  from  a  jet  engine  is  similar  in  mag¬ 
nitude  to  a  small  stationary  source  (F  -=*  10^m^/sec^). 


2. 


This  same  power  law 
cribing  plumes  from 
stacks  located  in  a 
lence  [Hoult  (1975) 


has  been  successful  in  des- 
high  temperature  has  turbine 
region  of  high  ambient  turbu- 
and  Egan  (  1975) ]  . 
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4.5.4  Event  Modeling 

Several  distinct  approaches  were  attempted  in  the  analysis  of  the 
Dulles  three-tower  data.  The  most  straight-forward  involved  fitting  a 
gaussian  vertical  profile,  plus  ground  reflection  term  to  the  concentration 
measurements  at  each  tower  to  obLain  the  plume's  centerline  location  and 
vertical  spread  for  each  event.  After  this  individual  fitting,  the  dynamical 
pi.ume  rise  and  growth  equations  were  fitted  to  the  earlier  obtained  values  for 
the  plume  centerline  and  oz.  The  advantage  of  this  approach  was  that  plume 
successful  for  the  first  tower  (i.e.,  closest  to  the  taxiway),  but  it  proved 
unreliable  at  the  more  distant  towers  where  plume  centerlines  were  often  above 
the  highest  receptor  and/  or  where  rapid  vertical  dispersion  produced  nearly 
uniform  vertical  concentration  profiles. 

At  the  other  extreme  lies  the  ensemble-fit  method,  where  the  entire  set 
of  observations  of  a  single  aircraft  type  under  the  full  range  of  meteorologi¬ 
cal  conditions  is  applied  to  a  single  comprehensive  theory  containing  a  number 
of  adjustible  parameters.  This  method  provides  a  starting  point  for  investi¬ 
gating  single  event  deviations  from  the  ensemble  predictions  but  may  obscure 
interesting  dynamical  effects  not  built  explicitly  into  the  model.  This  method 
was  chosen  above  the  single  event  method  because  of  the  fact  that  many  events 
had  a  "non-ideal"  distribution  where  a  centerline  maximum  was  not  observable 
at  even  the  first  tower.  Figure  illustrates  the  plume  rise  at  tower  1. 

Other  methods  like  alternate  multiparameter  schemes  for  assessing 
individual  events  were  considered  but  abandoned  as  their  numerous  parameters 
could  not  be  adequately  determined  by  the  data  accompanying  each  event. 


4.5.5  The  Ensemble  Model 

Consider  the  case  of  a  source  with  emission  rate  q  moving  at  velocity  V 
along  an  infinite  line  orientated  at  an  angle  P1  with  respect  to  the  positive 
y  direction  (see  Fig.  4.7).  If  the  wind,  u,  defines  the  positive  x  direction, 
the  receptor  is  located  at  (x,  o,  z)  and  t  =  0  corresponds  to  the  source 
position  (0,0),  then  the  instantaneous  concentration  at  the  receptor  is  given 
by  Eq  .  7  . 


q  exp 


:(t)  = 


j  1 

/  \  2  / 

/  x  -  ut\  / 

1  ~  2 

\  *  )  V 

2™”z  V 

V  cos 


n-\2 


V. 


S(z,H,r z) 


(7) 


Ill 


where 


S(z,H,o  )  =  exp  -  -r- 


1  (z  -  H 


2  V  a. 


1  /  z  +  H 


-H 


ve  =  |Ut  - 


and  assuming  cx  =  = o  does  not  vary  significantly  over  the  time  interval: 

-X-^-  <  t  <  ^-a- 
u  u 

While  several  assumptions  involved  may  be  avoided  by  a  more  tedious 
calculation,  they  do  not  appear  in  conflict  with  the  experimental  situation  at 
hand.  This  simple  formulation  then  leads  directiy  to  the  expressions: 

qS(z,H,oz) 

CPeak  "  27raaz~Ve_  (10) 


T  =  / 81n2 


u  \ V  cos0 ’ 


which,  with  minor  corrections  for  Ecolyzer  response  (Smith,  1977)  correspond 
to  the  measured  peak  concentration  above  background  and  the  pulse  duration  at 
502  of  peak  concentration,  respectively. 

The  dynamic  behavior  of  the  plume  was  assumed  to  be  described  by  the 
equations : 

a  “  ox( 0)  +  bx  ut  sin  (12) 


°z  =  t o - C 0 )  +  b,  ut  sin  aJ  + 


2  <H  -  H0)‘ 


H  =  H„  +  h(x  -  utn)P/uq 


H0  -  1.2  az(0) 


(15) 
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where 


Og  =  3  min  average  measurement  of  the  standard  deviation  of 
wind  direction, 

ox(0)  =  initial  along-wing  plume  spread, 
o2( 0)  *  initial  vertical  plume  spread, 

bx,bz  =  plume  growth  parameters.  They  describe  the  growth  of  the 

plume  relative  to  o  , 

6 

H  =  plume  centerline  height  at  distance  x, 

H0  =  initial  plume  centerline  height. 


The  addition  of  the  term 


(H  -  H0)2 


10 


is  suggested  by  Pasquill. 


The  Eq .  (14)  used  for  plume  rise  is  somewhat  more  general  than  the 
equation  suggested  theoretically  in  that  the  powers  p  and  q  are  free  parame¬ 
ters.  Fits  were  done  with  p  and  q  free  and  with  these  parameters  fixed  at 
p  =  2/3  and  q  ■  1  as  given  by  the  2/3  power  law  relation. 

Equation  15,  K0  “  1.2  oz(0)  is  dictated  by  the  assumption  of  zero 
vertical  concentration  gradient  which  causes  an  uniform  concentration  profile 
near  the  ground. 

With  these  Eq .  (9-14)  we  can  define  a  measure  of  "goodness-of-f it"  ^2 
and  via  minimization  of  the  eight  free  parameters  ox(0),  oz(0) ,  bx,  bz,  h, 
tQJ  p,  and  q  are  to  be  determined.  The  equation  for  y2  £s  given  by: 


where  TC  and  6F  denote  the  approximate  measurement  errors  of  0.25  ppm  and  10 
seconds.  The  superscripts  T  and  M  denote  theory  and  measurement  respectively 
The  indicated  summation  is  over  all  measurements  for  a  single  aircraft  type. 

The  preceding  expression  for  y2  should  actually  be  normalized  by  the 
total  expected  variances  in  Cpea^  and  (i.e.,  the  statistical  plus  the 
measurement  component)  and  not  merely  by  the  measurement  variances.  However 
the  statistical  variances  are  not  determinable  from  these  data  alone. 

This  shortcoming  preclude  determination  of  overall  model  confidence  level 
and  parameter  errors.-  Some  additional  insight  into  potential  model  improve¬ 
ments  is  provided  by  alternate  consideration  of  the  linear  correlation  coef- 


1 1  1 

ficients,  for  predicted  versus  observed  Cpeak  and  F,  and  the  associated 
confidence  bounds  on  these  correlations. 

4.5.6  Results 

Table  4.2  shows  the  number  of  events  for  the  different  aircraft  types. 
The  results  are  based  on  a  somewhat  smaller  selection  because  of  the  con¬ 
straint  that  the  wind  direction  was  within  70°  of  being  perpendicular  to  the 
taxiway.  All  events  were  observed  under  near  neutral  to  unstable  atmos¬ 
pheric  conditions,  with  bulk  Richardson  numbers  ranging  from  -] .0  *  10“^  to 
-0.02,  windspeeds  in  the  range  from  1.2  to  13.4  m/sec,  and  3-minute  Oq  from 
4.7  to  37.0  degrees.  Average  taxi  speeds  ranged  from  7.6  to  11.8  m/sec. 

Fixing  the  values  of  p  and  q  as  2/3  and  1  respectively,  the  ensemble  fit 
yields  the  parameters  given  in  Table  4.3.  The  correlation  values  for  the 
concentration  values  and  their  95%  confidence  limits  are  also  given. 

The  correlation  values  for  the  pulse  duration  is  not  given.  The  theory 
is  quite  poor  in  predicting  T,  the  pulse  duration.  This  is  partially  attribu¬ 
table  to  the  fact  that  F  is  predicted  to  be  independent  of  height  but  con¬ 
siderable  fluctuation  is  observed  experimentally.  Another  factor  contributing 
to  this  poor  correlation  is  that,  contrary  to  theoretical  expectations,  the 
observed  pulse  duration  increases  very  little  between  the  first  and  third 
tower  and  is,  in  fact,  consistent  with  zero  along  wind  plume  growth  (i.e., 
bx  =  0).  Thus,  the  differences  in  horizontal  growth  rate  factor  bx  between 
aircraft  types  cannot  be  considered  significant.  The  value  of  bz  is  found  to 
be  highly  correlated  with  the  H  parameter,  which  establishes  the  rate  of  plume 
rise  and  the  associated  vertical  dispersion  due  to  entrainment.  In  addition, 
despite  significant  differences  in  engine  placement  on  the  B707,  B727,  and  DC8 
aircraft,  the  initial  along  wind  and  vertical  plume  dimensions  are  nearly 
identical . 

The  rather  poor  value  of  the  correlation  coefficient  r  is  thought  to  be 
due  to  large  variations  in  ambient  air  stratification  and  its  effect  on 
buoyancy.  Evidence  for  this  is  obtained  when  each  event  is  optimized  sepa¬ 
rately  with  h,  the  parameter  which  describes  the  buoyancy,  as  the  only  free 
parameter.  The  other  parameters  were  fixed  to  their  values,  found  in  the 
ensemble  fit.  Figure  4.8  shows  the  wide  v  iety  of  h  values  for  the  B707  when 
calculated  as  described  above.  Significant  variations  about  the  ensemble 


Table  4.2.  Number  of  Events  pe-  Aircraft  Type 


A/C  TYPE 

EVENTS 

B707 

48 

B727 

25 

B737 

4 

B747 

7 

DC  8 

12 

DC9 

3 

DC  10 

17 

A37 

7 

L1011 

15 

SST 

10 

CESSNA 

3 

TURBO 

1 

Table  4.3.  Ensemble  Fit  Parameters 

AIRCRAFT  TYPE  ox(0)  oz(0)(ft)  bx  bz  h  tQ(sec)  r(95%  C.L.) 


1  16 

value  of  16  a re  seen.  The  correlation  coefficient  of  jumps  from  0  63  to 
0.81  when  this  event-by-event  freedom  is  permitted. 

Another  insight  into  event-by-event  deviation  from  the  prelicted  x2/^ 
plume  rise  behavior  can  be  seen  in  Fig.  4.9.  It  shows  the  wide  variety  of 
p-values  (the  plume  rise  trajectories)  one  obtains  when  p  is  determined  on  an 
event-by-event  basis.  Ensemble  optimizations  in  which  p  and  q  were  allowed 
to  be  free  indicated  somewhat  higher  values  for  p  and  q  of  1.25  and  1.75, 
respectively;  however,  the  resulting  improvements  in  r  were  insignificant. 

The  plume  rise  equation  used  in  the  above  analysis  does  not  predict 
final  plume  height  and  thus  is  of  limited  usefulness  in  terms  of  airport  air 
quality  models.  A  simpler  model,  described  briefly  below  has  been  applied  to 
the  3-tower  data  and  the  results  indicate  the  importance  of  i  itial  plume 
dilution  and  rise  on  observed  concentrations. 

Assume  that  plume  growth  is  governed  by  the  equations 

ox  =  oy  =  ox( 0)  +  bx  oyT(x)  (17) 

and 

°z  =  (°z(0)  +  bz  azT(x))2  +  °.l  <H  -H0)2]  l/2  ,  (18) 

where  ox(0)  and  oz( 0)  are  the  initial  horizontal  and  vertical  plume  dimensions, 
Oy'j’(x)  and  az'p(x)  are  the  dispersion  coefficients  taken  from  Turner  (1970) 
and  bx  =  bz  =  0.7  is  a  correction  for  averaging  time.  Further  assume  that 
the  final  plume  rise  is  given  by  the  equation 

H  =  H0  +  h/u,  where  H0  h  1.2  oz(0)  (19) 

u  is  the  wind  speed,  and  h  is  the  plume  rise  factor  to  be  determined.  Taking 
the  dynamical  behavior  of  the  plume  rise  into  account  and  applying  these 
equations  to  a  sample  of  121  cases,  encompassing  all  the  aforementioned 
aircraft  types  except  Concorde,  one  obtains  the  optimized  parameters 

ox(0)  “  60  ft 
oz(0)  =  26  ft 

h  =  386  ft2/ sec 

For  ground  level  data  points  (z  =  6  ft)  the  regression  equation 


cOBSERVED  =  m  ^THEORY  +b 
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fit  to  the  peak  concentrations  yields  a  correlation  coefficient  of  0.56,  a 
slope  of  m  =  0.96,  and  an  intercept  of  b  =  0.11  ppm.  A  scatter  plot  of 
observed  versus  peak  concentrations  is  shown  in  Fig.  4.10. 

If  the  parameters  describing  the  initial  plume  dimensions  and  plume 
rise  are  instead  taken  as  ox(0)  =  oz(G)  =  h  =  0,  the  resulting  regression 
parameters  are  m  =  0.11  and  b  =  0.52,  which  implies  that  for  the  highest 
observed  concentrations,  the  theory  is  overpredicting  by  a  factor  of  9. 

Hence  the  assumptions  about  initial  plume  size  and  rise  can  have  serious  air 
quality  modeling  consequences.  For  example,  the  initial  vertical  dispersion 
of  26  ft  is  equivalent  to  the  amount  of  dispersion  realized  by  1500  ft  of 
downwind  transport  under  F  stability  conditions. 


4.5.7  Conclusions 

Preliminary  analysis  of  measurements  of  the  CO  exhaust  plume  from 
taxiing  aircraft  suggest  that  the  rise  of  these  horizontally  injected, 
buoyant  plumes  is  not  inconsistent  with  the  2/3  power  law  relation  over  the 
distance  range  of  65  m  to  165  m,  but  large  event-by-event  fluctuations  from 
this  average  behavior  lead  to  rather  mediocre  correlations  between  theoreti¬ 
cally  predicted  and  observed  pollutant  concentrations. 

There  is  an  evidence  for  about  a  10  second  delay  in  plume  rise.  The 
present  analysis  suggests  that  the  plume  rise  delay  time,  t0  is  significant 
and  at  the  highest  wind  speed  has  the  effect  of  suppressing  plume  rise  at  the 
first  tower. 


Recalling  the  comparison  between  the  B707,  B727,  and  DC8 ,  the  initial 
plume  rise  during  taxi  was  found  to  be  characteristic  of  aircraft  dimensions. 
The  significance  of  engine  geometry  one  might  expect  was  not  observed. 
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buoyancy , 


determined  here 


where  the  total  buoyancy  flux  F  for  the  B707  and  B727  aircraft 
taxiing  is  about  150  nA/sec^  [Goldberg  (1978)].  A  value  of  0. 
used  in  this  calculation.  In  the  case  of  iet  aircraft  plumes 
tor  the  entrainment  constant  might  be  expected. 
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Fig.  4.10.  Scatter  Plot  Plus  Projection  Histograms  for 
Observed  Versus  Calculated  Peak  CO  at  Ground 
Level  Receptors  ( z  =*  6  ft.).  The  calculation 
uses  equations  17-19  and  the  optimal  parame¬ 
ters  given  in  the  text  . 
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Dulles  measurements  indicate  that  aircraft  contributions  to  hourly  CO 
concentrations  do  not  exceed  0.06  ppm  per  aircraft  at  a  distance  of  200  ft 
from  the  taxiway  centerline.  Thus,  concentrations  adjacent  to  a  busy  taxiway 
accommodating  100  aircraft  per  hour  should  not  exceed  6  ppm.  The  Dulles 
data  also  permits  determination  of  initial  plume  dimensions  and  plume  rise  for 
subsequent  use  in  airport  air  quality  models,  and  further  indicate  that 
failure  to  consider  these  initial  plume  parameters  can  lead  to  order-of- 
magnitude  overpredictions  of  local  concentrations. 

4.6  RE-EXAMINATION  OF  DULLES  SINGLE  EVENTS  (THC  VS  CO) 

Although  hydrocarbons  (HC)  are  not  the  main  impetus  of  this  program, 
they  are  perceived  to  be  a  major  component  of  the  air  quality  deterioration 
experienced  at  airports.  In  this  regard,  simultaneous  measurements  of  HC  and 
CO  from  taxiing  aircraft  at  Dulles  have  shown  a  high  degree  of  c  on. 

The  THC  measurements  at  Dulles  must  be  treated  with  the  same  reservations 
concerning  the  separation  of  reactive  and  non-reactive  species.  However,  when 
dosages  due  to  single  aircraft  are  considered,  the  methane  fraction  may 
most  likely  be  ignored.  The  peak  THC  concentrations  above  background  due 
to  the  passage  of  an  aircraft  average  1.48  ppm.  With  an  average  THC  dose 
per  aircraft  of  32  ppm-sec,  it  will  require  84  aircraft  operations  between  6 
to  9  A.M.  to  reach  0.25  ppm.  The  THC  is  measured  at  a  height  of  6  ft, 

215  ft  from  the  taxiway,  with  C  stability  as  an  average.  Adjusting  these 
measurements  to  reflect  worst  case  conditions  suggests  that  0.25  ppm  may  be 
approached  at  Dulles  with  ~ 40  operations  per  3  hours  under  such  adverse 
condi tons . 

In  order  to  evaluate  indirectly  the  decrease  in  THC  with  distance  due 
to  turbulent  dispersion,  coincidently  measured  concentrations  of  THC  and  CO 
are  compared.  According  to  published  emission  factors  for  the  taxi/idle  mode, 
the  CO/THC  emission  ratio  varies  from  1.1  to  <•'*'4  depending  on  aircraft  type.  A 
scatter  plot  (Fig.  4.1!)  of  peak  concentrations  of  CO  versus  THC  for  taxiing 
aircraft  events  exhibits  significant  correlation  (r  =  0.72)  and  a  CO/THC  ratio 
of  mean  peak  concentrations  of  1.26.  A  similar  analysis  of  the  mean  doses 
gives  a  ratio  of  2.1.  As  peak  concentrations  are  a  function  of  instrument 
response  time,  the  dose  ratio  represents  the  more  useful  figure.  The  dose 
ratios,  found  to  vary  more  between  individual  events  than  between  aircraft 


(i.e.,  engine)  types,  are  seen  to  agree  in  the  mean  with  the  engine  emission 
ratios  of  CO  to  THC.  From  these  single  station  comparisons  it  appears  that 
estimation  of  THC  through  scaling  of  CO  concentrations  is  viable  for  such 
near-field  experiments  where  the  background  may  be  separately  identified. 


5  CARBON  MONOXIDE  MEASUREMENTS  AT  A  HIGH  ACTIVITY 
FLY-IN  AT  LAKELAND  AIRPORT,  FLORIDA+ 

5.1  INTRODUCTION 

Carbon  Monoxide  (CO)  concentrations  were  measured  during  a  major  fly-in 
of  general  aviation  (GA)  and  experimental  aircraft  at  Lakeland  Airport, 

Florida  from  January  23  thru  29,  1978.  Over  3000  aircraft  participated  in 
this  fly-in,  where  in  excess  of  250  aircraft  operations  per  hour  were  experi¬ 
enced.  The  purpose  of  the  measurements  was  to  quantify  the  effect  of  emis¬ 
sions  from  GA  aircraft  on  air  quality  under  extreme  conditions  of  airport 
activity.  The  Federal  Aviation  Administration  (FAA) ,  in  conjunction  with  the 
Environmental  Protection  Agency  (EPA) ,  planned  the  measurement  activity. 

Three  FAA-owned  Energetic  Sciences  "Ecolysers"  (CO  monitors)  were  used.  EPA 
personnel  participated  in  the  field  program  and  assisted  in  data  gathering. 

Figure  5.1  is  an  aerial  view  of  the  airport  and  the  lightly  populated 
surrounding  countryside  and  Figures  5.2  and  5.3  show  the  operating  pattern  at 
the  airport  during  easterly  and  westerly  winds,  respectively.  Figure  5.4 
shows  the  location  of  the  25  monitoring  sites  used  at  one  time  or  another  in 
the  course  of  the  data  gathering.  Measurements  were  made  during  aircraft 
landing,  takeoff  and  taxi  modes.  Additionally,  an  instrument  was  set  up  in  an 
auto  which  was  periodically  driven  around  the  entire  airfield  at  the  periphery 
in  attempts  to  detect  gross  airport  contributions  to  the  local  CO  "background. 
Within  the  discriminating  capability  of  the  equipment  this  was  not  possible, 
nor  were  significant  observable  levels  of  CO  measured  during  any  of  some  50 
observed  landings. 

From  all  these  measurements  taken  under  a  variety  of  airplane  activity, 
and  meteorological  conditions,  the  maximum  projected  one-hour  average  concen¬ 
tration  measured  at  positions  where  people  might  be  expected  to  be  located  was 
less  than  2  parts  per  million  (ppm)  by  volume.  This  concentration  is  insig¬ 
nificant  (Federal  Register,  June  19,  1978)  when  compared  to  the  one-hour 
National  Ambient  Air  Quality  Standard  (NAAQS)  of  35  ppm.  These  measurements 
constituted  one  consideration  in  the  formal  recommendation  bv  the  EPA  to 
withdraw  GA  engine  emission  standards  (Federal  Register,  March  24,  1978). 


+Adapted  from  "Pollution  Dispersion  Measurements  at  High-Activitv  Fly-In  of 
General  Aviation,  Military,  and  Antique  Aircraft"  by  H.M.  Segal,  1978 
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Fig.  5.2.  Airport  Operations  During  Easterly  Winds  -  Lakeland  Airport 
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Fig.  5.3.  Airport  Operations  During  Westerly  Winds  -  Lakeland  Airport 
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Fig.  5.4.  Monitoring  Sites  -  Lakeland  Airport  January  23-30,  1978 
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Analysis  of  these  measurements  not  only  showed  the  small  influence  of 
GA  aircraft  on  air  quality  but  also  for  the  first  time  quantified  the  disper¬ 
sion  rate  from  individual  and  groups  of  GA  aircraft  that  were  landing, 
taxiing  and  taking-off.  Comparison  of  these  measurements  with  theory  indicat¬ 
ed  that  some  models  may  not  account  for  the  unique  dispersion  characteristics 
of  the  hot,  turbulent  airplane  emission  source. 

5.2  APPROACH 

To  fulfill  the  need  for  mobile  monitoring,  Ecolysers  and  associated 
strip  chart  recorders  were  placed  in  the  two  automobiles  used  in  this  monitor¬ 
ing  effort.  By  reinforcing  the  radio  antenna  with  a  dowel,  an  effective 
raobil  air  sampling  system  was  devised.  Figures  5.5  and  5.6  show  this  instal¬ 
lation  with  the  instruments  placed  inside  and  outside  the  automobile.  Where  a 
third  monitoring  site  was  required  the  installation  shown  in  Figure  5.7  was 
used.  Instruments  were  calibrated  with  20  ppm  calibration  gas  twice  per  day 
and  before  and  after  each  major  series  of  measurements. 

The  monitoring  procedure  was  similar  to  the  one  used  at  Dulles  Inter¬ 
national  Airport  to  monitor  Concorde  emissions  (Segal,  1977  and  Smith  et  al , 
1977).  Three  monitors  were  placed  in  a  line  abeam  to  the  path  of  a  taking-off 
or  taxiing  aircraft  so  that  these  instruments  could  sequentially  detect  the 
exhaust  plume  of  aircraft  emissions  as  the  plume  was  transported  by  the  wind 
over  the  three  stations.  In  addition,  as  mentioned  above,  a  car  instrumented 
with  an  Ecolyser  was  periodically  driven  around  the  entire  airport  in  attempts 
to  detect  the  gross  airport  contribution  to  the  local  vicinity  CO  background. 

Measurements  were  also  taken  during  landing,  takeoff  and  taxi.  Taxi 
measurements  were  made  at  both  (multiple  event)  and  low  (single  event) 
activity  times. 

5.3  RESULTS 

Landing  Measurements  -  On  January  27,  emissions  from  over  50  landing 
aircraft  were  measured  at  site  20.  (Figure  5.4).  Winds  were  from  the  north 
at  18  mph  and  "C"  stability  was  estimated. 

No  concentration  above  background  level  was  recorded  at  this  site, 
located  450  ft.  from  the  runway. 
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Takeoff  Measurements  -  On  January  29.  emissions  from  over  30  taking-off 
aircraft  were  measured  at  sites  20,  21,  and  22.  (Figure  5.4).  Wind  was  from 
330’  at  12  mph  and  a  "B"  stability  was  estimated. 

The  strip  chart  trace  from  a  typical  taking-off  GA  aircraft  is  shown 
in  Figure  5.8.  Emissions  from  this  modern  GA  aircraft  are  quite  low. 

However,  at  sites  20,  21,  and  22,  the  highest  pollution  levels  of  the 
entire  week  was  also  recorded.  This  occurred  when  a  World  War  II  vintage 
B-25  took  off.  The  CO  strip  chart  trace  of  this  take-off  is  shown  in  Figure 
5.9.  The  high  emission  levels  from  this  aircraft's  large  radial  engines, 
characteristic  of  both  military  and  commercial  aircraft  engines  of  that  time 
period  are  to  be  compared  with  the  almost  undetectable  pollution  produced  by 
the  Lurbine  engines  used  in  present  day  commercial  aircraft  (.Segal,  1977  and 
Smith  et  al ,  1977).  This  comparison  indicates  that  pollutant  emissions  have 
been  drastically  reduced  by  the  aircraft  industry  in  developing  the  gas 
turbine  engine  technology  of  the  present  era. 

Dispersion  measurements  permit  determination  of  a  power  law  exponent 
by  which  atmospheric  dispersion  may  be  parameterized.  This  dispersion  rate 
exponent  has  been  measured  during  airplane  taxi  and  takeoff  assuming  that  the 
relationship  between  concentration  and  downwind  distance  can  be  expressed  as: 

C  X~K  (1) 

where  C  is  the  concentration  at  downwind  distance  X.  The  rate  exponent  at 
which  the  pollutant  disperses  is  defined  as  K.  Peak  concentrations  of  those 
takeoff  events  having  adequate  signal  to  background  ratio  were  averaged  and 
were  found  to  disperse  as  in  the  power  law  expression  listed  above. 

This  exponent  which  is  derived  from  measurement  data  will  be  compared  with  the 
theoretical  value  of  this  exponent  in  Section  5.4. 

Taxi  Measurements  (Low  Activity)  -  On  January  26,  emissions  fiorn  over 
40  taxiin  a’. craft  were  recorded.  Wind  was  from  340°  at  15  mph  and  a  "C  ’ 
stability  ins  estimated.  Pollution  from  this  mode  dispersed  as  X~^-®  in  the 
previously  mentioned  power  law  relationship. 
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Taxi  Measurements  (High  Activity)  -  The  most  note-worthy  data  were 
obtained  at  station  17,  18  and  19  between  1700  and  1740  hours  on  January  28, 
when  a  continuous  queue  of  over  30  aircraft  stretched  down  the  taxiway  for  more 
than  1/2  mile.  Wind  was  from  343°  at  12.  mph  and  "D"  stability  was  estimated. 
During  peak  activity  one  airplane  taxied  by  the  monitoring  station  every  ten 
seconds.  As  they  approached  the  end  of  the  taxiway,  these  aircraft  were 
almost  continuously  dispatched  down  the  two  takeoff  runways  at  the  rate  of  278 
aircraft  per  hour. 

Because  both  taxi  and  takeoff  emissions  impacted  at  the  three  monitoring 
stations  downwind  of  the  taxiing  aircraft,  it  was  necessary  to  devise  a  method 
for  measuring  emissions  from  the  taking-off  aircraft  only.  This  was  accomplish¬ 
ed  by  moving  the  instrumented  auto  at  site  18  to  site  20  which  is  directly 
upwind  of  the  taxiing  aircraft.  Takeoff  concentrations  were  measured  at  this 
location.  This  move  was  made  after  sufficient  data  had  been  collected  at  site 
18. 

The  contribution  of  takeoff  emissions  to  concentrations  at  sites  17, 

18,  19,  and  20  was  modeled  and  calibrated  with  measurements  taken  at  site  20. 
This  takeoff  contribution  was  then  subtracted  from  the  total  concentrations 
measured  at  sites  17,  18,  and  19  to  identify  concentrations  directly  attribu¬ 
table  to  the  taxiing  aircraft.  These  data  are  plotted  in  Figure  5.10.  These 
multiple  event  taxi  emissions  are  the  found  to  disperse  as  in  the 

power  law  relationship  CaX~^. 


5.4  CONCLUSIONS 

The  following  conclusions  may  be  drawn  from  the  concentration  summary 
of  Table  5.1: 

1.  From  all  measurements  taken  under  extreme  aircraft 
activity  conditions,  the  maximum  recorded  con¬ 
centration  for  CO  at  the  closest  position  where 
people  might  be  expected  to  be  located,  was  less 
than  2  ppm  for  a  projected  one-hour  time  period. 

This  concentration  is  insignificant  when  compared 
to  the  one-hour  NAAQS  of  35  ppm. 

2.  The  highest  CO  concentration  ever  recorded  of  the 
dispersing  plumes  of  a  taking-off  airplane  (22 
ppm  at  335  ft.  from  the  runway  centerline)  was 
measured  at  Lakeland  Airport  on  January  29,  1978. 

This  measurement,  which  was  from  a  World  War  II 
vintage  B-25,  indicates  that  airplanes  have  been 
significant  sources  of  CO  pollution  in  the  past. 
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Table  5.1.  Carbon  Monoxide  Concentrations  During  Different  Operational  Modes 


MODE 

STABILITY 

CLASS 

WIND 

SPEED 

(MPH) 

CHARACTERISTICS 

PEAK  CONCENTRATIONS  PPM  ABOVE  BACKGROUND 

SINGLE  EVENT  1  HOUR  AVERAGE 

DIST. 

450  FT. 

DIST. 

335  FT. 

DISTANCE 
385  FT. 

LANDING 

(SINGLE  EVENT) 

C 

18 

<1  PPM 

TAKEOFF 

(SINGLE  EVENT) 

B 

12 

<1  PPM* 

23  PPM  (B-25) 

2  PPM  (GA) 

TAXI 

(SINGLE  EVENT) 

C 

15 

<  1  PPM 

TAXI 

(DURING  QUEUE) 

D 

12 

~2  PPM 

*  Jommercial  Jet  from  Dulles  Data,  Ref.  Smith  et  al,  1977) 


The  measured  dispersion  rate  exponents  during  taxi  and  takeoff  in  some 
cases  do  not  coincide  with  expectations  based  on  dispersion  rate  curves 
(Turner,  1970)  that  are  used  in  most  airport  models.  This  inconsistency  is 
important  to  recognize,  since  it  may  contribute  to  errors  in  receptor  concen¬ 
trations  calculated  from  airport  pollution  models.  While  this  short  term 
measurement  program  was  not  designed  to  develop  a  large  data  base  or  to 
explain  dispersion  inconsistencies  (m.b.,  that  no  measurements  of  the  vertical 
dispersion  of  the  emission  plume  were  made),  the  listing  of  the  dispersion 
parameters  in  Table  5.2  represents  an  initial  quantification  of  previously 
unmeasured  dispersion  characteristics  of  several  types  of  aircraft  exhaust 
plume  s . 

Factors  contributing  to  the  inconsistency  between  measurement  and 
theory  may  be  traced  ro  the  inability  of  the  theory  to  effectively  account 
for: 
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a.  Plume  rise 

b.  Extensive  initial  dispersion  related  to  the  turbulence 
field  created  by  the  high  velocity  fan  action  of  the 
propeller.  (The  extent  and  duration  of  this  turbulent 
field  is  unknown  at  the  present  time.) 

c.  Different  emission  densities  and  turbulence  intensity 
along  the  takeoff  path  of  an  accelerating  aircraft. 


Table  5.2.  Measured  versus  Theoretical  Aircraft  Plume  Dispersion  Rates 


Mode 

Approximate 

Aire  raft 

Speed 

Propeller 

Speed 

Measured 
Power  Law 
Exponent,  K. 

Theoretical* 
Power  Law 
Exponent,  K. 

Taxi 

(Single  Event) 

15  MPH 
(Constant) 

Low 

1.0 

1.8 

Taxi 

(During  Queue) 

5  MPH 
(Constant) 

Low 

0.4 

0.9 

Takeoff 

25 

(Accelerating) 

High 

1  .9 

1.8 
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6  THE  MEASUREMENT  OF  CO  CONCENTRATIONS  FROM  QUEUING  AIRCRAFT 
AT  LOS  ANGELES  INTERNATIONAL  AIRPORT+ 

6.1  INTRODUCTION 

Carbon  monoxide  (CO)  emissions  from  queuing  aircraft  were  monitored  at 
Los  Angeles  International  Airport  (LAX)  from  April  16  to  April  20,  1979. 

Carbon  monoxide  was  selected  as  the  pollutant  of  concern  because  it  is  stable, 
easily  measured  and  predominant  during  aircraft  queuing.  LAX  was  selected  for 
this  experiment  because: 

1.  It  is  a  busy  airport. 

2.  High  data  accumulation  is  possible  under  the  influence  of 
the  predominant  sea  breeze  which  blows  emissions  directly 
down  the  main  queueing  taxiway. 

3.  The  airport  authority  was  very  cooperative  and  permitted 
equipment  positioning  directly  on  the  taxiways. 

4.  A  National  Climatic  Weather  Station  which  records  wind 
direction,  speed  and  vertical  temperature  profiles  is 
located  within  1500  meters  of  the  monitoring  sites. 

5.  Data  from  this  program  could  be  compared  with  similar 
data  generated  in  the  early  1970s  and  which  was  used  as  a 
justification  for  the  aircraft  engine  emission  standards. 

6 . 2  APPROACH 

The  approach  of  this  program  is  to  measure  and  model  the  emissions  of 
aircraft  that  are  lined  up  (or  queued)  along  a  taxiway  just  prior  to  takeoff. 
Queuing  was  measured  au  both  the  north  and  south  runway  complexes  (Figure  6.1) 
from  April  16  to  April  20,  1979  ,  Monitoring  and  wind  measur ement  equipment 
were  positioned  directly  downwind  of  the  queuing  aircraft. 

Two  Energetic  Sciences  Model  2000  "Ecolysers"  were  employed  in  this 
program.  These  instruments  were  calibrated  with  20  parts  per  million  (ppm) 
calibration  gas  before  and  after  each  intensive  measurement  period. 

Equipment  was  placed  in  a  Federal  Aviation  Administration  vehicle  with 
the  pollution  sampling  tube  extending  outside  the  vehicle  where  it  was  at¬ 
tached  to  a  vertical  probe.  Air  intake  height  was  1.7  meters.  A  second 


+Adapted  from  "Emissions  from  Queuing  Aircraft"  by  H.M.  Segal,  1980. 
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monitoring  instrument  was  located  50  meters  downwind  from  the  first  monitoring 
location.  Its  air  intake  tube  was  attached  to  a  tripod  and  elevated  to  a 
height  of  1.7  meters.  (Results  from  this  second  monitoring  location  are  not 
reported.)  Figure  6.2  shows  instrument  layout.  Equipment  was  lined  up  in  the 
direction  of  the  prevailing  wind  which,  because  of  its  westerly  direction, 
transported  a  line  of  aircraft  emissions  directly  over  the  receptors.  This 
arrangement  provided  the  desired  worst  case  pollution  geometry.  Wind 
velocity  was  measured  every  15  minutes  at  the  monitoring  sites. 


Air  quality  was  recorded  during  162  minutes  of  aircraft  activity  during 
a  five-day  time  period.  Queue  lengths  varied  from  1  to  8  aircraft.  (Figure  3 
shows  the  configuration  of  one  7  aircraft  queue  that  was  monitored.)  Distance 
from  the  first  queuing  aircraft  to  the  nearest  receptor  was  220  meters  from  the 
south  runway  and  320  meters  for  the  north  runway.  Air  quality  was  recorded 
for  wind  speeds  of  2.8  to  8.6  meters  per  second  under  Pasqui 1 1-Gi f ford  stabil¬ 
ity  classes  of  B,  C,  D,  and  E.  Airplane  entrance  to  and  exit  from  the 
various  queue  positions  was  recorded  to  the  nearest  second.  This  precise 
recording  of  the  time  when  each  aircraft  entered  and  left  its  queue  position 
and  the  simultaneous  recording  of  pollutant  concentrations  at  the  downwind 
receptors  were  essential  portions  of  this  program.  One  person  was  assigned 
full  time  to  accomplish  these  tasks. 

Upon  completing  the  monitoring  program,  the  162  minutes  of  data  were 
stratified  according  to  wind  speed,  stability  class,  and  queue  length,  and  a 
flow  diagram  such  as  the  one  shown  in  Figure  6.4  was  prepared  for  each  of 
eight  different  queuing  conditions.  These  conditions  reflected  measurements 
taken  during  different  days,  wind  speeds,  and  stability  conditions.  Emissions 
dispersion  during  transport  to  the  monitoring  sites  was  then  modeled  for 
comparison  with  measurements.  Each  airplane  was  positioned  on  the  taxiway  in 
accordance  with  its  observed  location.  Data  from  Tank  and  Hodder  (1978)  were 
used  to  determine  the  height  of  the  plume  centeiline  and  the  initial  size  of 
the  plume.  Pollutant  transport  times  from  queue  to  the  receptor  location  were 
determined  by  dividing  each  sour ce-to-receptor  distance  by  the  measured  wind 
speed , 
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6.3  ANALYSIS 

Two  types  of  analysis  were  performed.  First,  emissions  from  a  number 
of  shor t~dui at  ion  queue  events  that  occurred  during  the  entire  monitoring 
period  were  analyzed  separately.  This  ensemble  of  measurements  was  then  used 
to  verify  the  model.  Then,  the  single  longest  queueing  event  of  19  minutes 
duration,  which  occurred  on  April  18,  was  analyzed  in  detail.  Emissions  from 
this  event  were  also  summed  with  emissions  from  all  other  aircraft  activity 
during  a  70  minute  time  period  surrounding  the  event.  This  longer  averaging 
time  calculation  permitted  model  performance  to  be  determined  for  the  longer 
averaging  time  (1  hr)  considered  in  the  National  Ambient  Air  Quality  Standard 
(NAAQS)  for  carbon  monoxide  (CO) . 

Three  types  of  data  were  considered.  First,  measured  concentrations 
were  averaged  over  the  short  duration  of  a  queue  event.  Secondly,  model 
estimates  of  concentrations  were  made  corrr spond ing  to  the  meteorological 
conditions  and  averaging  times  of  these  measurements.  Thirdly,  model  estimates 
were  made  to  reflect  worst  case  meteorology  for  comparison  with  the  NAAQS. 

The  NAAQS  limit  of  35  ppm  is  not  to  be  exceeded  more  than  once  per  year  (approxi¬ 
mately  one  time  in  10,000  or  the  99.99th  percentile)  and  the  meteorology 
accompanying  such  a  violation  usually  consists  of  stable  air  3rd  low  wind 
speeds.  (These  conditions  may  be  conservative  since  they  often  do  not 
coincide  with  peak  aircraft  activity  (i.e.,  longest  queues).)  The  receptor 
for  this  latter  calculation  was  also  relocated  to  reflect  a  characteristic 
location  where  public  exposure  might  be  expected. 

6.4  MODELING  CONSIDERATIONS 

The  FAA  Simplex  model  (Segal,  1979)  was  used  in  this  study.  The  more 
complex  airport  models  were  not  employed  because  of  their  spatial  and  temporal 
averaging  assumptions.  The  Simplex  model  employs  a  simple  point  source 
algorithm  with  provisions  for  selecting  a  particular  plume  height  and  initial 
plume  dilution  box  size  for  each  aircraft  being  analyzed.  It  employs  the 
classical  Gaussian  point  source  equation  and  the  Pasqui 1  1-Gi f ford  curves  from 
Turner  (1968).  For  consistency,  the  same  sampling  times  were  used  in  all 
cases.  A  power  law  expression  for  the  standard  deviation  of  plume  size 


^  .  i  i  '■  i  — mm 


9HRSBMBBHI! 


138 


az  =  Kz  ,  Oy  =  KyX®-^  was  assumed, 

where  Ky,  Kz  were  chosen  to  match  the  P-G  curves  at  x  =  1.0  km. 

This  functional  form  was  used  in  all  calculations  which  emcompassed 
stability  classes  B  thru  E.  Such  an  approximation  greatly  facilitated  calcula¬ 
tions  and  differed  negligibly  from  the  P-G-T  predictions  even  at  the  shortest 
distances  used  in  this  study. 

6.5  RESULTS 

The  results  are  displayed  in  Figures  6.5  and  6.6.  Figure  6.6  shows  the 
comparison  of  estimated  and  measured  concentrations  during  the  variety  of  wind 
and  stability  conditions  experienced  over  the  entire  monitoring  period.  The 
average  ratio  of  estimated  to  measured  concentrations  for  the  ensemble  of 
measurement  events  reduces  substantially  when  finite  values  for  plume  height 
and  initial  plume  size  are  used.  The  ratio  of  1.7  for  the  latter  condition 
is  within  the  factor-of-two  considered  in  determining  an  acceptable  level 
of  model  performance. 

Measurements  performed  during  the  longest  queue  were  analyzed  separate¬ 
ly.  The  results  of  this  analysis  are  shown  in  Figure  6.5.  The  horizontal  bars 
represent  the  times  during  which  each  airplane  occupied  a  particular  queue 
position.  The  queue  length  at  any  time  is  determined  by  summing  the  number  of 
horizontal  bars  crossed  at  any  particular  time.  Between  1742  and  1801  hours 
on  April  18  the  queue  length  increased  to  eight  airplanes.  During  this 
time  period  estimated  concentrations  were  3  ppm  while  measured  (average) 
concentrations  were  1.5  ppm.  A  second  comparison  was  then  made  for  the 
extended  averaging  time  period  (70  minutes).  Under  these  conditions,  the 
estimated  concentration  was  1.3  ppm  and  the  measured  concentration  was  0.9 
ppm.  Both  of  these  model  test  conditions  fall  within  the  factor-of-two 
criteria  for  determining  model  performance. 

A  final  model  calculation  was  performed  to  reflect  worst  case  meteoro¬ 
logical  conditions  for  comparison  to  the  NAAQS .  Worst  case  conditions  of  "E" 
stability  and  one  meter  per  second  wind  speed  were  assumed.  The  receptor  was 
relocated  750  meters  downwind  from  the  end  of  the  taxiway,  a  distance  which  is 
characteristic  of  where  people  might  first  experience  aircraft  emissions.  The 
uppermost  curve  in  Figure  6.5  is  a  plot  of  these  conditions.  When  this  curve 
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is  averaged  over  the  entire  70  minute  period,  an  estimated  concentration  of  4 
ppm  results. 

6.6  CONCLUSIONS 

The  model  appears  to  reflect  measured  concentrations  quite  well  under 
both  the  short  averaging  times  of  specific  queue  events  and  the  longer  averag¬ 
ing  times  associated  with  all  types  of  airplane  activity  at  the  end  of  the 
taxiway . 

It  is  interesting  to  note  the  ability  of  the  model  to  track  the  pollu¬ 
tion  peaks  and  valleys  during  the  period  from  1742  to  1801.  This  capability 
is  quite  impressive  considering  the  number  of  times  the  airplane  changes  its 
queue  position.  The  time  shift  between  the  two  curves  is  in  the  expected 
direction  and  is  probably  related  to  the  time  taken  for  the  high  velocity 
engine  exhaust  to  slow  down  to  ambient  conditions  prior  to  atmospheric  trans¬ 
port  and  to  the  actual  slowing  down  of  the  ambient  wind  field  by  a  line  of 
queuing  aircraft.  (This  latter  condition  has  been  observed  at  Washington 
National  Airport) . 

When  the  verified  Simplex  model  is  used  to  estimate  concentrations  at 
expected  populated  locations  during  the  highest  activity  hour  monitored, 
concentrations  of  only  4  ppm  from  aircraft  alone  result.  This  value  is  small 
when  compared  to  the  NAAQS  limit  value  of  35  ppm. 

This  assessment  indicates  that  a  simple  point  source  algorithm  can 
successfully  accomplish  the  "verification  by  parts"  procedure  suggested  by 
Turner  (  1979  ).  Parametenzat  ions  from  this  verified  model  may  subsequently  b 
incorporated  into  a  number  of  more  complex  models  if  validation  efforts  at 
other  airports  confirm  these  results. 

The  results  of  this  study  should  apply  to  other  engine  exhaust  gases 
and  should  be  particularly  useful  in  defining  the  queuing  concentrations  of 
engine  N0X  (NO2  under  a  high  ozone  environment).  This  can  be  accomplished  by 
merely  changing  the  model  inputs  to  reflect  N0X  rather  then  CO  emissions. 
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Fig.  6.4.  Airplane  Flow  Diagram 
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CONCENTRATION 
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1- A-  ZERO  PLUME  AND  BOX 

2- 7 -  TLUME  AND  BOX  PER'IR  MEASUREMENTS 

OF  PLUME  RISE  AND  DIMENSIONAL 
PLACEMENT  OF  ENGINES  ON  AIRPLANE 

3  -  WIND  SPEEDS  -  2.8  -  8.6  m/s 

4  -  "Q"  LENGTH  1-8  AIRCRAFT 

3  -  STABILITY  CLASSES  ■■  C,  C,  D,  E 

6  -  TOTAL  EVENT  QUEUING  TIME  -  48  MIN. 

7  -  AIRCRAFT  TYPES  -  707,  727,  737,  DC-9 

DC-8,  DC-10,  A-300,  L-101),  747 

8  -  WIND  DIRECTION  -  240  -  280° 

9  -  MONITORING  DATES  -  APR  16-20,  1979 

10  -  CLOSEST  RECEPTOR-SOURCE  DISTANCE  - 

217  METERS  SOUTH  R/W 
320  METERS  NORTH  R/W 

11  -  MODEL  -  FA/.  SIMPLEX 


MEASURED 

CONCENTRATION 

(PPM) 


Fig.  6.6.  Reconciliation  of  Measured  with  Estimated  Carbon 
Monoxide  Concentrations  During  Aircraft  Queuing 
at  Los  Angeles  International  Airport  (LAX) 
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7  ACCURACY  DEFINITION  OF  THE  AIR  QUALITY  ASSESSMENT  MODEL 
(AQAM)  AT  WILLIAMS  AIR  FORCE  BASE+ 

7.1  OBJECTIVES 

The  U.S.  Air  Force  must  have  an  analytical  tool  which  can  reliably 
predict  the  impact  of  their  airbases  on  air  quality.  Such  a  tool  is  essential 
for  assessing  the  significance  of  existing  airbase  emissions  in  relation  to 
National  Ambient  Air  Quality  Standards  (NAAQS)  and  for  estimating  the  impact  of 
proposed  facilities  and  missions  changes  as  is  required  by  federal,  state  and 
local  laws.  The  Air  Quality  Assessment  Model  (AQAM),  developed  by  ANL  for 
the  Air  Force  is  such  a  tool  that  meets  these  air  quality  impact  assessment 
requirements.  The  AQAM  model  treats  airbase  emissions  in  sufficiently  great 
detail  and  incorporates  a  state-of-the-art  treatment  of  pollutant  transport 
and  diffusion.  However,  until  now,  a  suitable  data  base  has  not  been  avail¬ 
able  for  determining  the  accuracy  limits  of  this  mathematical  computer-based 
mode)  . 

Ambient  air  quality  measurements  of  CO,  NO,  N0X  THC,  and  CH^,  and 
visibility  were  made  at  Williams  AFB,  nea^  Phoenix  from  June  1976  to  June 
1977.  The  objectives  were  (1)  to  determine  the  effects  of  local  aircraft 
operations  on  air  quality  and  (2)  to  provide  a  data  base  for  evaluation  of 
AQAM.  In  this  report  model  independent  statistical  analyses  of  the  data  and 
comparisons  of  AQAM  predicted  and  observed  concentrations  are  performed  to 
determine  the  accuracy  limits  of  AQAM. 


7.2  APPROACH 

The  AQAM  is  designed  to  predict  the  impact  of  airbase  operations, 
including  aircraft,  on  the  air  quality  in  the  vicinity  of  the  airbase.  It 
does  not  predict  background  concentrations  as  such.  Hence,  it  is  important 
from  a  model  evaluation  point  of  view  that  the  influence  of  the  local,  modeled 
sources  dominate  the  measured  pollutant  levels.  Spatial  and  temporal  depen¬ 
dences  of  these  local  sources  must  also  be  well  described.  These  considera- 


+From  "Analyses  for  the  Accuracy  Definition  of  the  Air  Quality  Assessment 
Model  (AQAM)  at  Williams  Air  Force  Base"  by  Yamartino,  et  al,  1980. 
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tions  played  a  significant  role  in  the  selection  of  Williams  AFB  at  the  site 
meeting  these  requirements  for  this  model  validation  effort.  Williams  AFB  had 
a  high  volume  aircraft  operations  ano  is  relatively  remote  from  an  urban 
area. 

Thirteen  months  of  hourly  average  concentrations  of  CO,  NMHC,  and 
N0X  monitored  at  the  five-station  network,  shown  in  Fig.  7.1,  constitutes 
the  data  base  used  for  assessing  the  predictive  capabilities  of  AQAM. 

Parallel  data  bases  of  hour-by-hour  meteorology  and  aircraft  activity 
were  utilized,  in  conjunction  with  the  standard  emissions  inventory  input 
to  AQAM,  to  compute  pollutant  concentrations  at  the  locations  of  inter¬ 
est.  To  define  the  incremental  AQAM  predictive  power  obtained  through  the 
use  of  higher  time  resolution  aircraft  data,  AQAM  predictions  were  made 
based  on  both  the  standard  AQAM  input  of  annual  total  aircraft  operations 
(referred  to  as  AQAM  I  predictions)  and  on  the  hour-by-hour  aircraft 
operations  mentioned  above  (referred  to  as  AQAM  II  predictions). 

7.3  IMPACT  OF  AIRCRAFT  ON  LOCAL  AIR  QULAITY 

As  seen  in  Table  7.1,  pollutant  levels  at  several  receptors  are  found  to 
depend  in  a  significant  way  on  aircraft  emissions  though  the  average  concentra¬ 
tion  impacts  are  small  relative  to  the  National  Ambient  Air  Quality  Standards 
(NAAQS).  In  all  cases  AQAM  overpredicts  the  percentage  role  of  aircraft  emis¬ 
sions  but  much  of  this  is  simply  due  to  the  background  levels  not  accounted 
for  in  AQAM.  The  fact  that  AQAM  overpredicts  the  absolute  role  of  aircraft  at 
most  stations  is  thought  to  be  related  to  the  model's  neglect  of  plume 
rise  and  plume  turbulence  enhanced  dispersion:  two  mechanisms  which  act 
to  reduce  concentrations  nearby  the  aircraft.  The  largest  observed  average 
daytime  impact  of  aircraft  occurs  at  station  4  where,  on  the  average,  aircraft 
account  for  36%  of  the  CO,  28%  of  the  NMHC  and  24%  of  the  N0X . 

Both  AQAM  predictions  and  measurements  agree  that  station  4,  atypical 
in  the  sense  of  its  close  proximity  to  buildings,  trees,  and  automobiles,  sees 
the  highest  concentrations:  a  factor  of  2-3  higher  than  station  1,  2,  3,  and 
3  collectively  in  the  cumulative  frequency  distribution  (CFD)  sense.  The 
failure  of  the  AQAM  to  correctly  reproduce  the  observed  rank  ordering  among 
stations  1,  2,  3,  and  5  is  also  thought  to  be  due  to  dynamical  factors  such  as 
the  neglect  of  aircraft  plume  rise  (which  clearly  leads  to  overpr ed i c t ion  of 


Table  7.1.  Williams  AFB  Aircraft  Emissions  Impact  on  Annual  Average 
Hourly  6AM-6PM  Concentrations 


Concent  ration 

Percent 

Total 

Wi thout 

Aircraft 

Aircraft 

Concentrat ion 

Aircraft* 

Contribution 

Contribution 

( ppm) 

(ppm) 

(ppm) 

CO 

Station  1 


Observed 

0.114 

0.105 

0.009 

8 

AQAM  II 

0.071 

0.040 

0.031 

44 

Station  2 

Observed 

0.134 

0.119 

0.015 

11 

AQAM  II 

0.106 

0.039 

0.067 

63 

Station  3 

Observed 

0.108 

0.093 

0.015 

14 

AQAM  II 

0.195 

0.064 

0.131 

67 

Station  4 

Observed 

0.362 

0.230 

0.132 

36 

AQAM  II 

0.345 

0.168 

0.177 

51 

Station  5 

Observed 

0.156 

0. 136 

0.020 

13 

AQAM  II 

0.099 

0.046 

0.053 

54 

NMHC 

Station  2 

consistent 

Observed 

0.123 

0.128 

-0  .  OOS^  with  zero 

0 

AQAM  II 

0.039 

0.019 

0.020 

51 

Station  4 

Observed 

0.215 

0.155 

0.060 

28 

AQAM  II 

0.209 

0.071 

0.138 

66 

NOy 

(concentrations  in  ppb) 

Station  2 

Observed 

9.44 

8.98 

0.46 

5 

AQAM  11 

3.68 

2.15 

1  .53 

42 

Station  4 

Observed 

15.1 

11.5 

3.6 

24 

AQAM  II 

8.4 

6.6 

1  .8 

41 

*Based  on  regression  of  pollutant  vs  AQAM  II  estimated  aircraft  emissions 
on  ground  level  line  sources. 
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CO  and  NMHC  at  station  3).  Finally,  using  the  computed  CFDs  for  off-base 
populated  areas  and  allowing  for  possible  underprediction  by  a  factor  of  2-3, 
one  concludes  that,  with  the  exception  of  the  6-9  AM  National  ambient  guideline 
concentration  for  reactive  hydrocarbons,  the  airbase  impact  is  negligible 
relative  to  existing  NAAQS . 

No  significant  difference  in  predictive  power  between  the  AQAM  I 
and  AQAM  II  has  been  found,  thus  extremely  detailed  time  histories  of  aircraft 
operations  do  not  have  a  significant  effect  on  the  model's  accuracy  (predictive 
power)  and  the  standard  AQAM  input  of  an  average  diurnal  distribution  of 
aircraft  operations  appears  adequate. 

7.4  PERFORMANCE  OF  THE  MODEL 

In  the  CFD  sense,  the  AQAM  predictions  for  the  upper  percentile 
concentration  range  agree  reasonably  well  in  magnitude  and  slope  with  the 
observed  concentration  distributions  (sample  case  seen  in  Fig.  7.2),  suggest¬ 
ing  that  the  model  simulation  encompasses  a  range  of  emission  and  dispersion 
conditions  comparable  with  reality.  At:  the  lower  concentration  percentile 
levels,  the  CFDs  are  often  orders-of-magnitude  different,  reflecting  the 
problem  of  absence  of  background  levels  in  the  AQAM  computations.  CFD  esti¬ 
mates  of  the  99.99  percentile  concentrations  (i.e.,  highest  hourly  average 
concentration  per  year)  of  -3  ppm  CO,  1-3  ppm  NMHC,  and  0. 1-0.3  ppm  N0X 
agree  surprisingly  well  with  observed  values  of  2-4  ppm  CO,  1-3  ppm  NMHC,  and 
0.08-0.15  ppm  NOx  if  stations  1,  2,  3,  and  5  are  considered  collectively; 
however,  such  estimates  for  any  single  station  may  underpredict  the  once  per 
year  high  by  as  much  as  a  factor  of  1.7  for  CO  and  NMHC  and  3  for  N0X .  The 
fact  that  the  CFDs  of  observed  concentrations  at  the  different  stations 
converge  at  the  upper  percentiles  while  the  individual  station  curves  diverge 
slightly  for  the  AQAM  predictions,  suggests  that  the  most  severe  pollution 
episodes  actually  exist  over  a  spatial  domain  much  larger  than  the  airbase  and 
thus  are  probably  not  solely  due  to  specific  local  sources  such  as  aircraft, 
as  suggested  by  the  model . 

In  examining  the  performance  of  AQAM  on  an  hour-by-hour  basis  one 
encounters  shortcomings  common  to  Gaussian  plume  models  in  general.  If 
no  accounting  of  background  pollutant  levels  is  made,  hour-by-heur  comparisons 
of  AQAM  with  observations  indicate  severe  underprediction  for  ail  three 
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pollutants  (a  mean  factor  of  3  for  CO  and  NMHC  and  a  factor  of  5  for  N0X). 

In  addition,  the  standard  deviations  of  these  distributions  indicate  that  the 
unadjusted  model  falls  short  of  the  "50-percent  within  a  factor-of-two" 
criteria  for  Gaussian  models.  However,  addition  of  a  modest  annual  mean 
background  (0.09  ppm  for  CO,  0.08  ppm  for  NMHC,  and  7  ppb  for  N0X)  leads  to 
a  dramatic  improvement  in  predictive  power.  The  background  adjusted  model 
yields  predictions  with  a  factor-of-two  of  observation  in  excess  of  65%  of  the 
time,  while  errors  in  excess  of  factor-of-ten  occur  at  a  tolerable  ~1%  level. 

The  reason  such  order  of  magnitude  discrepancies  exist  lies  with  the  fundamental 
limitations  of  modeling  a  stocastic  process  with  a  deterministic  model. 

7.5  DIFFICULTIES  WITH  THE  THEORY  VERSUS  OBSERVATION  COMPARISON 

At  the  time  the  experiment  was  being  planned  (circa  1975),  Williams  Ain 
had  the  highest  level  of  aircraft  operations  of  any  airbase  in  the  U.S.,  and, 
as  Williams  is  a  training  base,  it  was  expected  that  records  of  aircraft 
activity  would  be  more  accurately  maintained  than  at  other  bases.  While 
accurate  records  were  available  during  normal  training  operations  periods, 
documentation  of  off-hours  activity  (e.g.  weekends)  was  incomplete.  In 
addition,  as  most  of  the  operations  involved  small  twin-engine  aircraft  (i.e. 
T37,  T38,  and  F 5 ) ,  selection  of  the  airbase  having  the  highest  traffic  count 
was  not  necessarily  compatible  with  a  choice  based  on  highest  aircraft  pollu¬ 
tant  era i ss  ions  . 

It  was  also  thought  that  the  remoteness  of  the  base  from  other  signifi¬ 
cant  sources  would  render  the  resolution  of  airbase  and  aircraft  generated 
pollution  from  background  levels  straightforward.  Unfortunately,  Phoenix, 
though  some  50  km  to  the  Northwest,  contributed  high  background  levels 
to  the  measured  air  quality  particularly  at  night.  These  so-called  background 
levels  often  exceeded  the  local  pollutant  levels,  resulting  in  .  poor  signal- 
to-noise  ratio  and  greatly  reducing  the  effectiveness  of  the  receptor  network 
in  sensing  local  source  (i.e.,  airbase)  created  pollutant  gradients.  In 
addition,  the  entire  Valley  of  the  Sun  appears  at  times  to  exhibit  pollution 
reservoir  characteristics  which  can  not  be  predicted  by  a  short-range  Gaus¬ 
sian  plume  model  such  as  AQAM .  Even  the  several  hour  transport  and  disper¬ 
sion  of  pollutant  from  Phoenix,  though  included  in  the  AQAM  inventory  of 
environ  sources,  is  not  adequately  treated  due  to  total  reliance  on  the 
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stationary  state  assumption.  Such  multihour  transport  could  have  been  more 
realistically  modeled  using  a  backward  trajectory  technique,  which  would 
select  the  emission  rate  for  the  time  period  presently  impacting  the  receptors 
and  allow  for  varying  dispersion  rates  over  the  trajectory  of  the  plume,  but 
such  is  not  the  case  in  the  present  AQAM,  designed  for  short-range  pollutant 
transport  and  dispersion  calculations.  Thus,  it  was  necessary  to  attempt 
to  validate  the  model  under  conditions  where  a  major  portion  of  the  aerometric 
signal  was  related  to  distant,  background  sources  not  adequately  treated  by 
the  model.  The  presence  of  five  monitors  on  the  base  could  have  been  useful 
in  subtracting  out  these  unwanted  and  poorly  described  components  of  the 
observed  concentrations.  However,  two  factors  limited  the  effectiveness  of 
this  latter  approach  to  investigating  the  local  (i.e.,  airbase)  contribution 
to  the  observed  pollutant  levels.  First,  noise  in  the  form  of  spatial  inhomo¬ 
geneities  of  the  background  and,  second,  inter-instrument  random  and  systematic 
errors  which  tended  to  wash  out  many  of  the  more  subtle  effects  since  local 
signal  components  were  often  small  compared  to  the  accuracy  limits  of  the 
instruments . 

All  of  the  studies  of  model  predictive  power  versus  meteorological 
parameters  or  time  of  day  suggest  that  time  of  day  is  the  most  signifi¬ 
cant  variable  affecting  AQAM  performance  in  that  AQAM  reproduces  the  major 
trends  in  daytime  observed  concentrations  when  local  sources  dominate  but 
seriously  underpredicts  at  night  when  more  distant  sources  contribute.  This 
deficiency  is  probably  due  to  an  underestimate  of  vehicle  activity  between 
midnight  and  5  a.m.  and  to  a  breakdown  of  the  steady-state  Gaussian  plume 
assumption  used  in  the  model.  Major  revision  of  the  model  to  incorporate 
backward  trajectories  would  probably  be  required  to  rectify  this  latter 
problem.;  however,  such  a  revision  is  perhaps  of  only  academic  interest  at 
present  since  the  AQAM  is  most  successful  in  simulating  the  potential  "worst 
case"  airbase  impact  situations  associated  with  morning,  low  wind  speed, 
stable  or  low  inversion  height  conditions  coincident  with  the  commencement  of 
high  airbase  emissions. 


7.6  CONCLUSIONS 
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In  final  summary,  we  note  the  major  findings  of  this  study  as 

•  the  upper  concentration-percentile  range  of  the  CFD  of  AQAM 
predictions  agrees  well  enough  with  observation  to  make  this 
a  valuable  tool  for  predicting  rare  worst  case  situations. 

As  at  William  AFB,  the  estimate  of  worst  case  (e.g.,  highest 
hourly  per  annum)  concentrations  at  a  particular  station  may 
be  underpredicted  by  a  factor  of  two  or  three  so  that  comple¬ 
mentary  monitoring  (for  perhaps  a  <  w  .  indred  hours)  should 
be  performed  in  questionable  sit  ’  ■>  when  concentrations 

approach  50%  of  NAAQS . 

•  the  AQAM,  adjusted  for  background  levels,  provides  a  predic¬ 
tive  tool  for  hourly  average  concentrations  commensurate  in 
ability  (i.e.,  greater  than  50%  of  predictions  within  a  factor 
of  two  of  observed  concentrations)  with  well  calibrated  Gaus¬ 
sian  plume  models  applied  to  urban  areas. 

•  Without  adjustment  for  background  levels,  most  AQAM  results 
underpredict  by  more  than  a  factor  of  two  and  often  by  as 
much  as  10  to  100.  Without  background  compensation,  AQAM 
results  should  be  interpreted  with  caution.  Air  quality 
monitoring  data  sampled  in  the  Air  Quality  Control  Region 
should  be  considered  when  analyzing  AQAM  results. 

•  The  present  AQAM,  unadjusted  for  jet  plume  rise  or  other  jet 
plume  dynamical  factors,  tends  to  overpredict  the  near  field 
air  quality  impact  of  aircraft.  Studies  at  commercial  air¬ 
ports  suggest  that  this  overprediction  becomes  negligible 
beyond  a  few  hundred  meters  from  the  aircraft. 

•  The  highest  predicted,  average  hourly  aircraft  impact  of 
0.18  ppm  CO  (51%  of  total  CO)  at  station  U  only  slightly 
overprediccs  the  observed  impact  of  0.13  ppm  CO  (36%  of 
total  CO).  [Special  instances  do  exist,  however,  where 
aircraft  contribute  nearly  100%  of  observed  and  predicted 


concentr at  ions . 
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•  The  ability  of  the  AQAM  to  accurately  predict  inter-station 
concentration  differences  is  only  weakly  confirmed  because 
of  large  measurement  errors  relative  to  these  observed  con¬ 
centration  differences  and  because  of  the  unexpectedly  high 
background  concentrations. 

•  The  AQAM  could  benefit  from  minor  revisions  such  as  the  in¬ 
corporation  of  jet  plirne  rise  and  turbulence  enhanced  dis¬ 
persion  and  from  major  revisions  such  as  a  backward  trajec¬ 
tory  calculation  for  more  realistic  assessment  of  the  impact 
from  distant  sources. 

•  AQAM  is  ready  for  acceptance  under  EPA  Guidelines  on  Air 
Quality  Modeling. 
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